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I«   INTRODUCTION 

Although  many  thin-wire  computer  programs  have  been 
developed  for  the  purpose  of  analyzing  antennas  and 

E 

scatterers,  few  of  these  programs  have  been  directed  toward 
the  student  of  electro-magnetic  theory.  The  majority  of  the 
programs  are  directed  to  the  engineer  or  advanced  student 
for  the  purpose  of  analyzing  designed  structures  or 
verifying  experimental  data. 

The  purpose  of  the  study  is  to  develop  a  computer 
program  by  modifying  an  existing  computer  code  which  can  be 
utilized  as  an  educational  method  to  develop  insight  into 
radiating  structures  by  the  beginning  student  of 
electro-magnetic  theory. 

The  modified  Ohio  State  University  Antennas-Scatterers 
Analysis  Program  (OSUMOD  or  ASAP)  is  directed  toward  the 
beginning  student  who  does  not  yet  have  the  expertise 
necessary  to  manipulate  the  input  data  for  proper  execution 
of  the  larger  more  comprehensive  analysis  program.  Even 
though  ASAP  is  small  in  core  requirements  and  is  fast  in  run 
time,  it  is  capable  of  analyzing  structures  to  assist  the 
engineer  with  design  problems. 

Since  the  resulting  program,  ASAP,  is  primarily  directed 
toward  students,  the  program  has  been  limited  to  structures 
which  contain  less  than  50  monopoles  (segments) ,  no  longer 
than  one-fourth  of  a  wavelength,  and  which  have  less  than  51 
nodes  (intersections  and  endpoints) .   If  a  ground  plane, 
either  perfect  or  finite  is  present;  the  stated  limits  above 
are  halved  due  to  the  generation  of  an  image  structure. 


II*   ORIGINAL  PROGRAM 

A.  THEORY 

Reference  1  presents  the  electro-magnetic  theory  for  the 
analysis  of  antennas  and  scatterers  in  an  isotropic,  linear, 
and  homogeneous  ambient  medium.  The  analysis  is  performed  in 
the  freguency  domain  with  an  excitation  caused  by  either  a 
generator  or  an  incident  wave. 

In  the  analysis,  a  piecewise-sinusoidal  expansion  is 
used  for  the  current  distribution.  The  matrix  eguation 
Z  I  =  V  is  generated  by  enforcing  reaction  tests  with  a  set 
of  sinusoidal  dipoles  located  in  the  interior  region  of  the 
wire.  Since  the  current  distribution  has  the  same  form  as 
the  expansion  mode,  this  formulation  is  known  as  the 
"sinusoidal  reaction  technique". 

B.  COMPUTER  PROGRAM 

Reference  2  presents  the  computer  program  corresponding 
to  the  theory  presented  in  Ref.  1. 

1.  Input  Format 

In  the  program,  the  input  data  must  specify  the 
frequency,  wire  radius,  wire  conductivity,  the  parameters  of 
the  exterior  medium,  coordinates  of  the  points  to  describe 
the  shape  and  size  of  the  wire  configuration,  a  list  of  the 
wire  segments,  and  the  indicators  for  the  various  outputs. 
Table  1  is  the  input  data  necessary  to  analyze  a  half-wave 
dipole. 

2.  Output  Format 

In  the  original  form,  the  only  outputs  which  could 
be  requested  by  the  input  data  stream  are  the  following: 


a.  Antenna  Problems 

(1)  Current  Distribution  on  the  Structure. 

(2)  Input  Impedance. 

(3)  Radiation  Efficiency. 

(4)  Near-Zone  Field. 

(5)  Far-Zone  Field. 

b.  Backscattering  Problems 

(1)  Absorption  Cross  Section. 

(2)  Scattering  Cross  Section. 

(3)  Extinction  Cross  Section. 

(4)  Complex  Elements  of  the  Polarization 
Scattering  Matrix 

c.  Bistatic  Scattering  Problems 

Echo  Area. 

Table  2  is  an  example  of  the  output  data  available  for 
data  of  table  1. 

3.   LIMITATION 

Although  the  program  can  analyze  a  structure  with  up  to 
50  segments,  55  points  and  60  dipoles  modes;  it  can  not 
analyze  a  structure  in  the  presence  of  a  finite  ground 
plane. 
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Ill*   MODIFIED  COMPUTER  PROGRAM 

A.  Input  Format 

As  illustrated  in  table  1  the  format  for  the  input  data 
cards  is  not  self  explanatory.  This  format  can  be  determined 
by  referring  to  the  FORMAT  statements  of  the  program  of  Ref. 
2.   Since  the  modified  program  is  directed  toward  the 
student,  the  input  data  format  was  changed  to  allow  free 
format.  Reference  2  was  written  in  a  form  which  permitted 
modifications  to  allow  flexibility  in  specifying  input  data 
for  the  analysis  program.   Appendix  B,  titled  "User's 
Manual",  discusses  the  input  data  cards  necessary  for  proper 
execution  of  an  analysis  problem.   Appendix  B  is  self- 
contained  and  may  be  used  independently  of  the  remainder  of 
this  document. 

B.  Output  Format 

In  the  original  computer  program,  the  absence  of  labels 
encumbered  the  output  data  and  lessened  the  usefulness  of 
the  program.  To  improve  the  usefulness  of  the  modified 
version,  detailed  labels  were  added  to  the  output  data.  As 
with  the  input  data,  Ref.  2  was  written  in  a  form  which 
enabled  modification  to  allow  more  specific  output  data  for 
the  analyzed  problem.  With  the  addition  of  the  polar 
plotting  package,  the  far-zone  electric  field  intensity 
polar  radiation  and  reradiation  patterns  can  be  plotted.  A 
sample  problem  can  be  found  on  page  120  in  Appendix  B, 
User's  Manual. 

C.  Finite  Ground 

To  enable  the  student  or  the  engineer  to  have  an 
improved  analysis  program,  the  finite  ground  effects  were 
added  to  ASAP.  The  theory  corresponding  to  the  ground 
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effects,  which  utilize  Fresnel  reflection  coefficients,  is 
discussed  in  Appendix  A,  titled  "System  Manual".   Also 
discussed  in  Appendix  A  is  the  modified  computer  program  and 
the  corresponding  theory.  The  electro-magnetic  theory  was 
developed  in  Refs.  1,  2,  and  3;  and  it  is  restated  with  its 
corresponding  computer  code  to  assist  in  the  understanding 
of  the  methods  applied.   Appendix  A  is  self-contained  and 
may  be  used  independently  of  the  remainder  of  this  document. 


IV*   CONCLUSION 

The  addition  of  ground  effect  technigues  to  the  original 
program  did  not  alter  the  accuracy  or  the  computational 
capabilities  of  the  program.  The  ground  effect  technigues 
utilized  the  results  of  the  original  program  and  modified 
these  results  to  account  for  the  effects  of  the  presence  of 
the  finite  ground. 

To  verify  the  numerical  results  of  ASAP,  the  input 
impedances  of  both  a  horizontal  and  a  vertical  dipole  were 
compared  to  the  solutions  of  the  exact  form  of  the 
Sommerfield' s  eguation.   As  can  be  seen  in  table  3  the 
finite  ground  treatment  of  ASAP  agrees  favorably  with 
Sommerfield • s  solutions.  The  ASAP  finite  ground  results  are 
also  in  excellent  agreement  with  the  previous  computer 
solutions  of  Refs.  4  and  5. 


V«   RECOMMENDATIONS 

Although  the  program  is  a  general  analysis  tool  for 
students,  several  future  modifications  will  enhance  the 
program  as  a  design  tool  for  engineers.  These  items  include: 
varying  the  wire  radius  on  the  structure;  incorporation  of 
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non-radiat ing  elements  such  as  transmission  lines;  varying 
the  wire  insulation  radius,  conductivity,  and  dielectric 
constant;  ana  a  geometry  generation  package  such  as  dipole 
array  or  helix.  One  ma~jor  change  that  would  both  improve  the 
speed  and  reduce  the  core  requirement  is  that  of  symmetry. 
No  attempt  was  made  to  utili2e  the  symmetry  in  the 
admittance  matrix  when  the  ground  plane  is  present.  If 
symmetry  were  applied,  the  structure  size  limit  with  the 
ground  plane  present  would  be  approximately  that  of  the 
structure  without  the  ground  plane. 


VERTICAL  DIPOLE 
FREQUENCY    3MHZ 
LENGTH    .5  WAVELENGTH 
RADIUS    .005  METERS 
DIELECTRIC  CONSTANT  (RELATIVE) 
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CONDUCT IVITY 
.1 

.001 

.00001 


HEIGHT/WAVELENGTH 

.25 

.30 
.35 
.45 
.25 
.30 
.35 
.45 
.25 
.30 
.35 
.45 


ASAP 

108.43+J  63.22 
92.43*J  37.73 
82.59*J  34.76 
74.0UJ  40.11 
98.74*J  63.52 
86.15*J  39.02 
80.15*J  38. 0B 
75.42»J  41.79 
97.94+J  64.86 
86.12*J  40.85 
80.82+J  38.89 
75.93*J    41.55 


EXACT* 

113.7  *J  56.22 

93.97*J  32.62 

83.56*J  29.74 

74.87*J  35.34 

99.86»J  46.43 

87.40*J  34.04 

8l.ll*J  33.16 

76.32*J  37.03 

96.68*J  49.11 

87.37*J  35.94 

81.79*J  33.99 

76.83+J  36.78 


HORIZONTAL  DIPOLE 
FREQUENCY    3  MHZ 
LENGTH    .5  WAVELENGTH 
RADIUS    .001  METERS 
OIELECTRIC  CONSTANT  (RELATIVE) 


10 


CONDUCTIVITY 
.1 


.001 


.00001 


HEIGHI/WAVELENGTH 

.5 
.3 
.1 
.5 
.3 
.1 
.5 
.3 
.1 


ASAP 

72.95*J  28.53 
103. 70O  58.27 
30.16*J  69.54 
74.23*J  35.84 
94.16*J  51.32 
54.52*J  62.89 
76.13*J  35.99 
9  1.8  5* J  5  3.21 
54.30»J    56.81 


EXACT* 

78.25»J  23.91 

117.5    »J  52.05 

33.200  73.65 

80.66*J  31.84 

105.3    *J  46.01 

61.31*J  64.33 

82.790  31.73 

103.0    *J  48.56 

61.03+J  57.27 
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SYSTEM  MANUAL 

INTRODUCTION:       The  Antennas-scatterers  Analysis  Program 
(ASAP)  for  thin  wire  structures  in  a  homogenous  conducting 
medium  performs  a  freguency  domain  analysis  of  antennas  and 
scatterers.   The  program  is  applicable  in  the  presence  of  a 
ground  either  perfect  or  finite.  This  appendix  will  describe 
the  computer  program  which  accomplishes  this.   Although  the 
program  was  written  for  the  IBM  360  computer  system  it  can 
be  executed  on  another  system  with  minor  modifications. 

A  piecewise-sinusoidal  expansion  is  used  for  the  current 
distribution.   The  matrix  eguation  ZI  =  V  is  generated  by 
enforcing  reaction  tests  with  a  set  of  sinusoidal  dipoles 
located  in  the  interior  region  of  the  wire.   Since  the  test 
dipoles  have  the  same  current  distribution  as  the  expansion 
modes,  this  may  be  regarded  as  an  application  of  Galerkin's 
method.  Sumsey's  reaction  concept  was  most  helpful  in  this 
development,  and  therefore  the  formualtion  is  known  as  the 
"sinusoidal  reaction  technigue". 

The  main  routine  and  each  subroutine  is  discussed 
separately  in  this  appendix.  The  writeups  for  the 
subroutines  are  arranged  alphabetically  by  subroutine  name 
after  the  main  program.  Each  of  the  discussions  includes  the 
purpose  of  the  subroutine,  brief  description,  and  a  listing. 
After  the  subroutine  writeups  is  a  table  of  the  more  common 
symbols  used  in  this  program. 

The  input  data  and  program  limits  are  discussed  in 
detail  in  the  next  appendix  titled  "USERS  MANUAL". 
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GROUND  EFFECTS:      In  the  modified  antenna  analysis 
computer  program  finite  and  infinite  ground  effects  were 
added  by  using  the  reflection  coefficient  technique.   The 
method  in  which  tnis  technique  was  used  required  the 
generation  of  an  image  structure.   In  this  section  the 
reflection  tecnniyue  will  be  discussed  in  detail. 

In  order  to  apply  ground  effects  to  the  electric  field, 
the  field  for  the  image  structure  was  first  calculated  as  if 
a  ground  were  not  present.  Then,  the  field  was  decomposed 
into  parallel  and  perpendicular  components.  (A  parallel 
component  is  the  component  which  is  parallel  to  the  plane  of 
incidence.  A  perpendicular  component  is  one  which  is 
perpendicular  to  this  plane.  The  plane  of  incidence  is  the 
plane  containing  the  normal  to  the  reflecting  surface  and 
the  incident  ray.) 

Z  1 
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Consider  an  image  monopole  with  the  electric  field  in 
the  1  direction.   The  ray,  6,  is  a  vector  which  is 
perpendicular  to  1  and  passes  thru  the  point  of  interest.  To 
apply  reflection  technique,  the  plane  of  incident  must  be 
found.  It  is  advantageous  to  define  a  new  coordinate  system 
(a#/3#z)  where   and   are  parallel  to  the  xy  plane  with  a   in 
the  plane  of  incident  and  /3  perpendicular. 


If  the  direction  cosines  (cosx,  cosy,  and  cosy)  are 
known,  it  can  be  shown  that  the  components  of  the  field  in 
the  a/3  (xy)  plane  have  the  following  relationship: 


II 


sin  <p        -cos  i> 
where  <)  =  arctan  (cosy/cosx)  . 


cos  $   sin  <|> 


Now  the  reflection  coefficients  for  the  interface  can  be 
applied  as: 

„N(R)      »  J1 

E       =  R  E 

II 

1(H)  J_ 

E       =  R  E 


where  R    and  R   will  be  defined  later  in  this  section, 
II      1 

Applying  the  matrix  equation  above  yeilds: 
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*(R) 

E 

X 

= 

(R) 

E 

y 

cos    <J>        sin    (J) 


sin    $        -cos    <t> 


I  I  (R) 


±(R) 


(the  square  matrix  is  unique,  in  that,  the  inverse  is  equal 
to  the  oriqinal  matrix) .  Since  the  imaqe  direction  is 
opposite  to  the  oriqinal  monopole,  that  is. 


(1  X  z)  =  -  (1   X  z)  . 

oriqinal  imaqe 

the  z  component  of  the  field,  which  is  in  the  plane  of 
incident,  is  qiven  by: 


(R) 


-  R    E 
II   z 


From  electro-maqnetic  theory  the  reflection  coefficients 
for  the  fields  in  medium  (1)  at  the  interface  with  another 
medium  (2)  are  defined  as: 

for  parallel 


cos 


9  -  \/  6»  -  sin*  e 


H  cos   9    +     \/  &•    -   sin2    9 

and    for    perpendicualr 


R 


-    &•    cos   9    +     \/  &•    -    sin2    9 


&•    cos   9    +     \/  S«    -    sin2   9 


where   9    is    the   anqle   of    incident   as   measured   from    the    normal 
to    the   interface   and 


&•    =    (S       +     a  /j6j)/(6      +    a  /ico) 

2  2  11 


where  the  subscripts  correspond  to  the  mediums  above. 
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To  determine  the  relationship  between  E   ,  R   and  R  ,  R 

I  I    J_       V    H 

a  perfect  ground  (£   =0,  a  =  oo )  was  investigated. 

r 

limit  R  -   -1 
H 

limit  R   =  +1 
V 

But,  for  a  perfect  ground  the  contributions  to  the  field 
from  the  image  monopole  would  be  equal  to  the  field  of  the 
original  monopole  but  opposite  in  sign  due  to  the  chosen 
reference  direction, 

I  I  (R)         II       II 

E       =  R    E    =  -  E 

II 

J_(B)  J_       X 

E       =  R    E    =  -  E 

X 


therefore 


R   =  R 
II     H 


R    =  -  R 
J_        V 


In  summary,  the  contribution  to  the  electric  field  of  a 
monopole  over  a  ground  plane  at  a  given  point  is  given  by: 

(R)      (R)  (R)  (R) 

E    =  E    cosx  +  E    cosy  +  E    cosz 
x  y  z 

where 

(R) 
E     =  R   E  cosx  +  (R    -  R  )  E  cosx  cos2  $ 

x  J.  MX 

+     (R         -    R    )     E    cosy    sin    $    cos    <t> 
I  I  X 
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(R) 
E  =   R         E   cosy    -    (R        -   R    )     E    cosy   cos2    $> 

y  M  1 1        _l_ 

+    (R        -   R   )     E    cosx   sin   <t»    cos    <J> 
II         _L 

(R) 

E  =    -    R         E   cosz 

z  I  I 

where   E   is   the  field   without    the   ground   plane  present   and 


-   &•    cos  e   +     \/  &•    -   sin2   e 

a      =  — — 

I  I  &•    cos  e   ♦     \/  &•    -    sin2   e 


COS 


6   -     y &'   -  sin2    e 


X  cos   9    +     v/  &•    -   sin2    0 


&•     =    &       -    j  (CT/&  cu) 
r  o 
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MAIN 


PURPOSE:       to  control  the  input,  output,  and  the  flow  of 
calculations. 

METHOD:        The  main  program  controls  the  flow  of  the 
required  calculations  by  calling  only  a  few  subroutines. 
These  subroutines  in  turn  call  other  subroutines  which 
actually  do  the  required  calculations.  The  order  of  the 
calling  sequence  is  diagramed  after  the  listing  for  the  main 
program. 

The  DIMENSION  statements  at  the  beginning  of  the  main 
routine  provides  the  required  storage  for  a  wire  structure 
with  up  to  50  segments,  60  nodes  and  60  dipoles  without  the 
presence  of  a  ground  plane.  If  a  ground  plane  is  present 
one-half  of  the  reserved  storage  is  required  for  the  image, 
therefore  a  wire  structure  with  up  to  25  segments  and  30 
nodes  can  be  analyzed. 

NM  denotes  the  actual  number  of  monopoles  (segments) , 
INM  is  the  corresponding  dimension,  and  the  dimension  for 
CG#.  VG,  and  ZLD  is  twice  INM.  The  second  subscript  for  MD 
always  has  a  dimension  of  4  to  correspond  to  the  number  of 
segments  meeting  at  a  given  node. 

N  denotes  the  number  of  simultaneous  linear  equations 
and  ICJ  is  the  corresponding  dimension.  The  dimension  for  C 
is  (ICJ  *  ICJ  +  ICJ)/2. 

In  the  statements  above  statement  4,  the  initial 
conditions  and  defaults  are  established.  After  calling 
subroutine  READ  to  determine  the  input  parameters,  the  IF 
statements  output  the  parameters  to  be  used  for  the 
calculations.   In  the  DO  LOOP   ending  at  statement  1,    the 
the  input  data  of  the  structure  geometry  is  stored  in  order 
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to  recall  if  the  structure  is  to  be  moved  for  ground  plane 
calculations. 

After  the  image  structure  is  generated  and  structure 
location  is  moved,  subroutine  SORT  is  called  to  determine 
the  dipole  modes.   Prior  to  calling  SGANT,  the  load  and 
generator  information  is  established. 

Subroutine  SGANT  is  then  called  to  calculate  the 
elements  of  the  impedance  matrix.  If  FEEDS  or   GENERATORS 
are  specified  by  the  input  data  stream,  subroutine  GANT1  is 
called  to  solve  for  the  current  distribution  due  to  these 
forcing  functions. 

In  the  DO  LOOP  ending  with  statement  29,  subroutine 
GNFLD  is  called  to  calculate  the  near-zone  field  for  the 
current  distribution  of  the  subroutine  GANT1. 

The  subroutine  GFFLD  is  called  for  the  far-zone  field  of 
the  current  distribution  of  the  subroutine  GANT1  in  the  DO 
LOOP  ending  at  statement  35.  The  subroutine  GFFLD  is  called 
again  in  DO  LOOPs  ending  at  statements  42  and  51,  if 
bistatic  and  backscattering  calculations  are  reguested  by 
the  input  data  stream. 

CALLS  TO:  GANT1 
GFFLD 
GNFLD 
POLPRT 
READ 
SGANT 
SORT 
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D  I  ME  NS  I  UN    XI60I,     YI60I,     Z(60),     XGI60),     YG(60I,     ZGI60I  1 

DIHbNSIKN     11160).     12(601,     13(601.     JAI631,     JBUO),     KFLAGI30I  2 

DIM1NS10N    CPHI(SOO),     CTHFTI500)  3 

DIMENSION    i)AIY  1(3601,     DAfY2(i60),    DATY31360I,     DATY4(360I  4 

niHlNSUiN    01501,     UlSOli     IBI50),     ISC  I  50 1  .    M0|50,4I,    NOI50I  5 

DIMENSION    120(60),     KGENI60)  6 

DIMtNSION    XNPI50I,     YNPI50I,     ZNPI50)  7 

COMPLEX    CI  18301  8 

COMPLEX    CDATl  I  500  1  ,CDAT2l  500  ),CDAT3(500)  .CDAT4I500)  9 

COMPLEX    C J(60I,£P<60I  . EPP ( 6 0 1 . E T t 60 1 , E T T  160  I  10 

COMPLEX    CGOI  50)  ,SG0(60)  ,CG(  lo6)  ,VG(i6oi  .ZLOtlOO)  11 

COMPLEX     VOIM60J  ,  ZLL0(60)  12 

COMPLEX    fcPPS.EPTS,ETPS,ETTS,EX,EY(EZ  13 

COMPLEX    EP2,EP3,EP4,ERR,ETA,GAM,Yll,ZIl,ZS  I* 

OAIA    PI ,TP/3. 1*159,6.28318/  15 

OAIA    E0.U0/8.854E-12, 1.25666-6/  16 

NGEN    »    -I  17 

IGRO    -    -1  18 

LOAD    -    -1  19 

BM    *    -1  20 

1CAR0    "0  21 

AM    *    -1  22 

IFLAG    »     0  23 

VOLT) II     •    I 1..0. I  2* 

HGI    •    0.  25 

NM    «    0  26 

NP    '    0  27 

MSG    -    0                                                                                                         ,'  28 

SIG2    -    -I.  29 

102    *    -1.  30 

SIG3    ■=    -I  31 

ER3    -     I  32 

T03    «    0.                                                                                             .  33 

CMM    '    50.  3* 

ER2    -     1 .  35 

FMC    *    300.  36 

INM    =    50  37 

ICJ    -60                                                               •  38 

WRITE     (6,741                                                  .     .  39 

*  40 

DO    2     I' 1 ,30  41 

KFI.AGU  I    -    -1  42 

43 

DO    3    J"l  ,  INM  45 

ISCU)     -    0  46 

VGI Jl     •     ( .0, .01  47 

ZLDUI     •    (  .0,-3)  48 


J  J  ■  Jl 

INM 

VGUJI 

=  (  .0 

ZLDI J  J  1 

-  (  . 

NFFP  - 

0 

NBI  P  « 

0 

NPAP  » 

0 

AFFP  ■ 

1000. 

AFFT  - 

1000. 

ABIP  - 

1000. 

ABI  1  >= 

1000. 

ABAP  ■ 

1000. 

ABAt  - 

1000. 

STFP  » 

1. 

,  .01 
0..0) 


49 
50 
51 

w 

54 

55 

56 

57 

5B 

59 

60 

61 

62 

KNM     =0  63 

CALL     "FAD     (  IA,  IB,  I8ISC.ICAR0, IGA1N, IGRO,  1NEAR, INT  ,  ISCAT .  IWR,  IFLAG.  64 

IKFLAG,KGEN,10AD, LZO,MSG,NBAP,NBI P,  NF  FP,  NGEN  ,NM,  NP,  ABAP, ABAT, AFFP, A  65 

2FFT,AlllP.A1)IT,A1,BM,CMM,bR2,ER3,ER4,FMC,HGT,PHAF,PHAI,PHlF,PHll,PH  66 

3SF.PHS1 ,  tHAF, Thai, TH IF, IHII ,THSF,THSl,SlG2,SIG3,SIG4,T02.T03, VOLT,  67 

4X.XNP.Y ,YNP, I ,Zl  lO.ZNP.STEPl  68 

WRITE     (6,561  69 

IF    1MSG.LT. I )    GO    TO    5  70 

IF     IMSG.tO.II     WRITE     (6.70)     KFLAGOO)  71 

IF    IIFLAG.EU.4J    GO    TO    I  72 

5  IF  (IFLAG.  Fl).  5)  STOP  73 
IF  (AM.LT.OI  WaiTE  (6,1271  74 
IF  (AM. L  t .01  GO  TO  6  75 
IF  ( (NM.GT.OI .AN0.INP.GT.01 I  GO  TO  7  76 
WRI  IE     16,1161                                                                                                                                                                                 77 

6  IF  I IFLAG. EU.l  )  GO  TO  1  78 
MSG  -2  79 
GO    TO    4                                                                                                                                                                                   BO 

7  WRI TE  (  6, I  14  I  81 
WRI FE  I  6, 1  I  3)  82 
WRITE  (6.1121  81 
IF  IKFLAGIll.EU.il  WRITE  (6,83)  FMC  B4 
IF  (KFl AGI 2) . EQ. I i  WRITE  (6.B4I  AM  85 
IF  IKFL AGt 3I.E0. 1  I  WRITE  (6,851  CMM  86 
IF  (KFLAOI20) .NE  . I  )  WRITE  (4,871  87 
IF  (KFl A&( 41 .E3. I j  win  16,861  88 
IF  (KFl  AG  4)  . Eg.  I  )  WRITE  (6,88)  BM  89 
IF  IKFl AClS) .t 0. I  I  WRITE  16,891  SIG2  90 
IF  (KFl AGI6I.F 0. I i  WRITE  16,901  ER2  91 
IF  IKFLAi.(7l.ia.||  WRITE  (6,91  TD2  92 
IF  (KFlAGlal.Nt.il  WR1TI  I6,92|  93 
IF  ifcFlAGI9l.EQ.lj  WPITf  (6,931  SIG3  94 
IF  (KFLAGIIOI.tO.il  WRIlt  16,941  ER3  95 
IF     IKFLAGI  II  1.E3.1  I     WRITE     (6,95l     T03                                                                                                            96 
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I  F 
IF 
IF 
IF 
IF 
IF 


WRI  IE 
I  F     (K 


I  I  I 

UKl 

IE 


WR|  IE 
IF     (K 

IF  IK 
WRITE 
WRITE 
WRI  TO 
WHITF 


I! 

I  F 

IF 

1  I 

I  F 

IF 

IF 

IF 

IF 

If 

I  F 

IF 

FHZ 

OMEGA 

IF    (S 


IF 

IF 

IF 

IF 

IF 

IF 

ETA 

GAM 

IF    IK 

NPG 

NMG 


FL  AG( 

1  G  V  0  . 

IGPU. 

ll»U. 

lid'). 

FL  At.  t 

I  6,  I 

FL  ACI 

I  I    I  .X 

I  6,  I 

FL  Al,( 

FLAG! 

I  6,  I 

FL  AG( 

FL  AG  I 

16,1 

16.1 

(6,9 

16,1 

FLAG! 

FL  Aoi 

FL  AG  I 

F  L  AG  I 

FLAG! 

FL  AG  I 

FF  P.,. 

Fl   T.L 

BAP.L 

BAT. I 

Bl  P.L 

Bl  T.L 

IBISC 

FMC» 

-     TP 

IG2.L 

02.  LT 

IG3.L 

03. LT 

GRD 

GRO.G 

FLAG! 

CSQR 

OMEG 

FL  A&T 

NP 

NM 


NE.l  I  WRITE  16,12?) 
I  .AND.  (KFLAGI  251  .EQ.l  I  I 
I .AND. IKFLAGI  25) .EO.l I  1 
I .AND. (KFl AGI 251 .EU.ll  I 
I .AND. I KFLAGI  251  . Ew.  I  I  I 
EO.l  I     WRITE     (6.1211     INT 

C3.ll     WRITE     16,96 
I  I  I  ,YI  IBI  1  I  I  ,  Zl  IB 


WRITE 
■*R  I  T  E 
WRITE 
WRITE 


(6,1231 
16, 1251 
(6, 12*) 
16,1261 


ER*,SIG* 
HGT 


HR  ITE 

WRITE 


I .GT  .0  1    WRITE     (6,12 
I.GT.O)    WRITE     16,97 
I 
I 


I      ll,IA(ll,X(IA(lll,YIIA(   I  I  I  ,  Z(  IAI   I 

1 1 1 )  ,  I  ■  1 ,  NH) 

91     ILZDI  I  I .ZLLDI  I  I , I'l.LOAD) 
81     ILZDI  I  I.ZLLOI I  )  ,  l'l.LOAOl 

01     (KGENI  Il.VOLTI  11,1-l.NGENI 
I     (KGENI I I.VOLTII liI-ltNGENI 


.N 

.t 
I  .E 
I  .1 
I.E 

I  .E 

500 
503 
500 
500 
500 
500 

r.o 

E6 
HZ 
0.1 

) 

0.  I 

.  I 

II 

I.G 

uo/ 
csa 

I.N 


I  WPITE 
i  WRI  IE 
I  WRITE 
I  WRI  It 
I  WRITE 
I  WRITE 
WRITE 
WO  I  TE 
WRITE 
WRITE 
WRITE 
WRITE 


.  I 
.1 
I  .AND.  I  ISCAT.LT 


16.1101 

(6,99) 

(6,100)     PHA|, PHAF , THAI , THAF, STEP 

16,101)    Prill, PHIF.THI |,  THIF,  STEP 

(6,1021     PHSI ,PHSF .THSI.THSF.STEP 

16.103)     IXNPl   I  I.YNPI  I  I  .ZNPtl  )  ,  1-1,  I  NEAR) 

(6,1651     AFFP 

(6,10',)     AFFT 

16,109)     ABAP 

16,108)     ABAT 

I  6,  1071     ABIP 

(6,106)     ABIT 

01)    WRITE    (6,73) 


>T.L 


EP2-ER2»E0«CMPL 
EP2    ■    CMPLXI ER2» 

EP3«ER3»E0»CMPL 
EP3    ■    CMPLXI  ER3» 
EP*    •    CMPLXI ER*» 
ERR    "    EP*/EP3 
T  .01    WRI  TE     16,121  I     INT 
EP3) 

RT(-U0»EP3I 
"     1  I    GO    TO    9 


XI  1..-TD2I 
E0.-SIG2/0MEGA) 
X(  1..-T03I 
E0.-SIG3/0MEGA) 
E0,-SIG*/OMEGA1 


00   8    1-1. NPG 


97 

98 

99 

100 

101 

102 

103 

I0<> 

105 

106 

07 

08 

09 

10 

II 

2 

3 

<> 

5 

6 

17 

18 

119 

120 

121 

122 

123 

12* 

125 

126 

!2' 

128 
129 
1  JO 

Hi 
133 
13* 
135 
136 
137 
138 
139 
1*0 
Ml 
1*2 
1*3 
1** 


c 

c 

8 

X&l II  ■  XIII 
YGI 1)  ■  Y( II 
ZGI 11  ■  Zl II 

9 
10 

DO  10  1  '1 ,NPG 
XIII  •  XGl 1 1 
Yll 1  -  YGI 1 ) 

nil  •  zgi  1) 

NP    »    NPG 

NM    .     NMG 

IF     (  I&B  O.LE .01     GO     TO     15 

SFT    UP     IMAGE     FOR    GROUND    PLANE 

ZM1N    •     Zl 1 ) 

I  ML    •   0 

K     ■    0 

on   11    i = i >np 

IF     (Zl  I  I.LT.ZMIN)     ZMIN-ZUI 

Zl  II     ■     Zl  I  MH&T  

||     IZI I  I .GT.l .E-601    GO    TO    11 
IWL    ■    IWL»1 

11  CONTINUE 

IF     (ZMI N.LE.HGT I     GO    TO    12 

WRI IE     (6,1171 

IF     ( IFLAG.EO.l I     GO    TO     1 

IF     I IFLA&.E0.2)     STOP 

HSG   «    2 

GO    TO   * 

12  00    13    J=l,NM 

K     «     J»NM 

If'iiai  jkgtIiwli    iaiki-ia(ji»np-iwl 

13  1B(K)  -  IBU)»NP-IWl 


c 
c 

IWLP  -  IWL»I 

1* 

DO  1*  1 -IWLP 
J  .  |  »NP-  |WL 
XIJ)  ■  X( 11 
YUI  -  Vl  1  I 
ZIJI  •  -Zl II 

NP 

KNM  •  NMtl 
NM  -  2»NM 

1*5 
1*6 

1*7 
1*8 
1*9 
150 
111 
152 
153 
15* 
155 
156 
157 
158 
159 
160 
161 
162 
163 
16* 
165 
166 
167 
168 
169 
170 
171 
172 
173 
17* 
175 
76 
77 
78 
79 
80 
81 
82 
83 
8* 
85 
86 
87 
88 
89 
90 
91 
192 


26 


15 


17 


18 
l<3 

20 
21 


HP  -  2» 
CALL  SO 
IF  (MAX 
WRITE  I 
IF  1  IFL 
IF  I IFL 
MSG  *  2 
&CJ  10  * 
IF  (MIN 
WR1IE  ( 
IF  (|FL 
IF  ( IFL 
MSG  «  2 
GO  10  * 
WRITE  ( 
IF  (MAX 
1  12  '  I 
IF    1L0A 


HP  -  I  W I 

RT     (  I  A,  Id  .  I  1  >  I  2,  I3,JA,JB,MQ,N0,NN,NP,N,MAX,HIN,ICJ,  INM  I 

.LE.*)     GO     10    16 

t,m 

AG.Cd.l I     GO    TO     1 
AG.EU.2l     STOP 


.GE.l I    GO    TO    IT 
6 ,  7  2  I 

AG.FO.ll     GO    TO     1 
AG.Ew.2l     STOP 


6,561 

.'.!.<.. OK. H1N.LT.1.0R.N.GT.ICJ)     GO    TO    5* 

D.GT.OI    30    TO    19 


DO    18     l=l,NM 

aoi  i )   ■   lo.  ,o.  i 

IF    INGEN.GT.O)     GO     TO    21 


UO    20 
VGI  1  I 


=  1  iNM 
I0..0. 


KN    .    NM 

IF     (IGRO.GT. 0)     KN    ■    NM/2 

J    =     I 

ANTENNA     CALCULATIONS 

IF     (LOAO.LE.O)     GO    TO    24 

IF    IKFlAGI2tl.GT.OI    GO    TO    22 

00   23    JM.KN 

22  00   23    I  -I, LOAD 
■    LZDI I i 

II  IAIJI.EQ.K.)  .AND.  IKFLA&I l*I.GT.O) )     2L0IJI-/LL0I 
(KFL AG! 2*1 .GT.O)     2  L0(  K  1  -  2LL0  (  I  I 

I  I KF LAG  I I  *  I. GT.O I. AND.  (  I GRD.GT  .0  1  )     ZLD(J»KNI-ZID 
(  (KFIAG I2*).GT.0).AND.<  IGRO.GT  .0)  )     I  LD  I  K  »KN)  ■  210 
( IGRO.LE.OI    30    TO    23 
(  ( K.  IE  .  I WL I .AND. I KFL AG< I  4  I. GT.O I )    ZLL)(J»KNI-(0., 

II  IA I  J). IE. I  ML  I. ANO.t KFL AG! 2* I. GT.O) I     ZLD1K»KN)« 

23  CONTINUE 


I) 


0.1 
(0..0.) 


2*    IF    INGEN.LT.O)    GO    TO    27 


193 
19* 
195 
196 
197 
198 
199 
200 
201 
202 
203 
204 
205 
206 
207 
208 
209 
210 
211 
212 
213 
214 
215 
216 
217 
218 
219 
220 
221 
222 
223 
22* 
225 
226 

I21 
228 
229 
230 
231 
232 
233 
23* 
235 
236 
237 
238 
239 
2*0 


VGI J1«V0LT(  I  I 
'I     VGU»KNI-VG(.JI 


,0.1 


KN    -    NM 

IF     I  IGRD.GT.O  I     KN     "    NM/2 

J     ■     1 

IF    (KFL AG(23).GT.OI    GO    TO    25 

DO    26     J>I , KN 

25  00    26     1=1  ,IIGEN 
K     -     KCl  Nl   I  I 

IT     II   IA  I  Jl  .EU.K I  .ANO.IKHAGtm.&T.OM 
IF     (KFLAGI  231.31  .01     VGIK) -VOLTt  I  I 
IF     I  (Kf  LA&( I  3I.GT.0I .AND. I  IGRO.GT. 01 
IF    I IGRO.LE.OI     JO     TO    26 

IF     ( (K. IF  .  IHLl .ANO. (KFLAGI 131 .GT .0)1     VG I  J »KNI - 1  0 . , 0  I 
IF     (I  IA  Ul  .1  E.IWL)  .ANO.IKFLAGI  23I.GT.0I  I     VG  (  K  f  K.N  I  -  (  0  . 

26  CUNTINUE 

27  CALL     Sl.ANT     <  IA,  I  B,  INM  ,  INT  ,  I  SC,  I  I  ,  12,  I  3,  JA,  JB , MD , N , NO.NM , NP , AM ,BM, C 
1  ,CGn,CMM,[j.EP2.EP3.ETA,FHi,GAM,SGD,X,Y,2,iLD,2S,ERR,lGRD) 

IF  IN.(.  T.O)  GO  TO  28 
IF  (IFLAG.EU.2I  STOP 
MSG    ■    ? 

IF     IIFlAC.Ea.il     GO    TO    1 
GO    TO    '. 

28  I  F    (NG( N.LE.OI    GO    TO    36 
WR1 TE     ( 6,  f5) 

MR  I  IE  16,  761 
MR  I  IE  (6,77) 
WRITE     I  6,82) 

CALL    .,.'.'<  I  I     I  I  A,  IB,  INM,  IWR,  I  1,  12  ■  13,  I  12,  J  A  ,  JB,  MO,  N  ,  NO.NM,  AM.C  ,  C  J  ,CG 
1  ,CMM.D,EFF  ,GAM.r,&,Cr,U,SGO,VG,Yl  1  ,  21  1  ,  2L D , 2S , I GRO I 
WR1 IE     16,571     EfE ,GG,2ll 
NEAR    Fl El  0 

IF     (INfAR.LE.O)     GO    TO    30 
WR  I  TE     16,751 
WRITE     I  6, 781 
WHI  TF     16,77) 

(10    29     l>l,  1NEAR 
XP    •     XNP(  I 
YP    .    YNPI  |  I 
IP    ■    (NCI  I  ) 

CA1  L     GNFI  0    (  IA.IB,  INM,  II  . 
I X , Y . I «  XP  r VP, IP, E  X , E ¥ , E  2 , 1  &R 
WRITE     (6,58)     XP, YP, IP 
WRITE     I  6,591     EX, EY.E2 

29  CON! INUE 


12, 13, MO 

GRO.ERR  I 


,N,ND.NM,AM,C&0,SGO,ETA,GAM,CJ,D, 


2*1 
2*2 
2*3 
2** 
2*5 
2*6 
2*7 
2*8 
2*9 
250 
251 
252 
253 
25* 
255 
256 
257 
258 
259 
260 
261 
262 
263 
26* 
265 
266 
267 
26B 
269 
270 
271 
272 
273 
2  7* 
275 
276 
277 
278 
279 
280 
2B1 
282 
283 
28* 
285 
286 
207 
288 


27 


FAB    FIELD 

30 

IF     I l&A IN.LE .01 

DO    31     1=1, 360 
PATY1 (  1  1     •    0 
OAIY?( 1 1     .    0 
UATY3I 1  i     -    0 

GO    TO    36 

31 

OATY*l  1  1     *    0 

WRI1E     It, 751 
WRI  TE     1 6, 791 

ni'  1  It     16,771 
WRITE     (6,82) 
1  NC    "    0 

NPL     =    -  1 

IF     (KKAGII6l.EQ.il     WRITE     (6,69) 

I  F     INFFP.EQ.l  I    GO     TO    32 

NPHA     *      (  PMAF-PHAI  l/STEPH 


NTHA    »     ( THAf -THAI 

1/STEPM 

GO    TO    3* 

32 

IF     (AFF I.&T.500.) 
NPL    ■     1 
NPHA    -     360 
NTHA    ■     1 
PHAI     «     0. 
THAI     *     AFFT 
STEP    •     1. 
GO    TO     3* 

GO    TO    33 

33 

NPl     ■     2 
NPHA    -     1 
NIHA    -     360 
PHAI     -     AFFP 
THAI     •     0. 
STEP    •     1. 

3-. 

PH    '     PHAI -STEP 
00    35    KM, NPHA 
PH    ■    PH»S I E  P 
TH    =     THAI-SUP 
00    35     1 -l.NTHA 
PHSPH    ■     0. 
PhSIH    -    0. 

IH 


-    TH»SIEP 
IF     UICRO.GT.O) 
CAl 


CALL    GFFl  0     I  I  A,  | 
I •C,CGD|CG,CJ,CMM.u 

2&PP,GTT,PH,sto,SCSP,SCST,SPPM,SPIH,STPM,STTM,TH,X,Y,l,ZLD,7.S,ET».C. 
3AM, ERR, tGROI 


289 
290 
291 
292 
293 
29* 
295 
296 
297 
298 
299 
300 
301 
302 
303 
30* 
305 
306 
307 
308 
309 
310 
311 
312 
313 
31* 
315 
316 
317 
318 
319 
320 
321 
322 
323 
32* 
325 
326 
327 
328 
329 
330 


ANO.l ( TM.GT .90) ,ANO.( TH.LT .270) I )     GO    TO    35  331 

B.INC.  INM.JWR.I  1,  12,1  3,  I  12,M.01N,NDINM,»M.ACSP|ACST  332 

i.O.ECSP.ECST.EP.tf  .EPP.ETT.EPPS.EPTS.ETPS.ETTS.GC,  333 


EIMAG    •    CABS(ETTS) 


333 

33* 
335 
336 


IF 
IF 
IF 
IF 


35 


36 


3/ 


NPL  .EO.l  ) 
(NPl  .  EU.2I 
(NPl .EQ.2 I 
(KFl  AG(  16)  . 
60) 


EPMA&    =     CAHSIEPPSI 

irdlMAG.GI.I.E-32)  PHSTH=5  7.295779»ATAN2<AIMAG(ETTSI.REAL(ETTS) I 
I  KEPMAG.G'.  l.E-  32)  PHSPH=5  7.29  57  79»ATAN2(AIMAG(EPPSI,REAL(EPPSI) 
IE     (NPl.EJ.ll     OATYHKI  =EPMA& 

DATY2(K(=ETMAG 

[IATY1I  ll-EPHAG 

0ATY2I  I  )  *ETMAG 
,NE .11    GO    TO    35 
WHITE     (6,60)     TH.PH.GTT, GPP, E T T S , E IHAG, PH STH , EPPS , EPMAG , PHSPH 
CONTINUE 

WH 1 Tt     ( 6,561 

It      (NPl  .1  E  .0  )     Gil     TO     36 

CALL     PllPRT     (l.UATYl) 

CALL     POLPRT     (2.0ATY2) 

BAIK    SCATTERING 

IF     ( ISCAI.LE.O)     GO    TO 

WRITE     (6,75) 

16,80) 

I  6,77) 

(6,82) 


5* 


WRI 
WP  I 
WRI 
L 


TE 

TE 

TE 

0 


NPL 

INC    = 


-I 
I 


f     (NBAP.EO.  I  )    GO    TO    37 
NPtll     '     i  PHIF-PHI  |  l/STEP*  1 
NIHI      •=      (  TH|f  -TH1  |  1/StEPM 
IE     I IWR.LE.OI     WRITE     (6,621 
GO    TO    39 
I  F     IAeAT.GT.500. )    GO    TO    38 


*i 


N°L     = 

1 

NPHI 

•     360 

NIHI 

-     1 

PHI  1 

•    0. 

THI  1 

•     ABAT 

SIf  P 

■     1. 

GO    TC 

1     39 

NPL    ■ 

2 

NPHI 

•     I 

NIHI 

-     300 

PHI  1 

•     ABAP 

THII 

•     0. 

SIEP 

-     1. 

PH    - 

PH  1  1 -STFP 

00    *2    C«l ,NPH| 

PH    • 

PH»SIEP 

TH    - 

TH1I-SIEP 

337 
338 
339 
3<>0 
3*1 
3*2 
3*3 
3** 
3*5 
3*6 
3*7 
3*8 
3*9 
350 
351 
352 
353 
35* 
355 
356 
357 
358 
359 
360 
361 
362 
363 
36* 
365 
366 
367 
368 
369 
3  70 

3" 
372 

373 

3  7* 

375 

3  76 

377 

378 

379 

380 

381 

382 

383 

38* 


23 


OC  4  2  I •! ,NTH| 

In  -  THtSIEP 

IF  (  I  li.RD.f.I  .01  .ANO.I  I  TH.GT.90)  .AND.<TH.LT.270)  I  )  GO  TO  42 

I  •  L»l 

CALL  CFFLD  (  I  A,  16,  INC.  INM,  IWR,  I  1  .  12,1  3,  I  I  2, MO ,N, NO, NH, AM. AC SP ,ACST 
l,C.C(.t),tG.CJ,CMM,0,E  CSP,ECST,EP,Er,EPP,ETT,EPPS,EPTS,ETPS.EfTS.GG, 
2ii'(,,Gll  ,PH,SGD,SCSP,SCST,SPPM,SPTM,  STPM,STTM,TM,X,Y,7.,ZLO,ZS,ETA,6 
3AM.IMH, 1GR0I 

IF  (  IkO.O!  .01  GO  TO  40 

IF     INPl.LT.O)     NSITE     (6.631     PH , TH , SPPM, SP TM, STPM, STTM.AC SP , ACST , ECS 
I  P.ILSI,  SCSP.SCST 
3    C  Fill  (LI     -     PH 

CTML  Mil     .     Ill 

CUAM IL I     ■    CPPS 

C0AT2I l I     ■    EPTS 

CDAI M I  1     •    t  IPS 

CUAI4I  I  )     ■    FII5 

IF     1NPL  .NE.I  I     SO     TO    4  1 


-.1 


42 


43 


44 


45 


LAKM  EPPSI 
C AbSI  E  PIS) 
CAUSIE  TPSI 
CABSIE ITS) 

C  ABSIEPPSI 
CABSI.EPTSI 
CABS!  E  TPSI 
CAbSIETTS) 


0ATY1 Ik  I 
UAIV2IKI 
DATY3IKI 

Uaiv4(k  ) 
o(i  TO  42 
ciATru  I  | 
OAl Y2( I  I 

0ATY3I  I  ) 

ua  I  v.  1 1  I 

CUNT INUE 


WK| IE     I t, B>) 

IF     INPL  .LE.O)     r,0    TO    43 

CAIL     PfirRT     I  7  ,  DA  TV  1  ) 

CAIL    POLPRT     I8.DATY2) 

CALL    PULPRT     (9,DATY3) 

CALL    POLPRT     llOiDAmi 

I F     IKFLAGI  I  71 .NE.I)    GO    TO    45 

WRITE     (6,64) 


DO 
MR  I 


I'l.L 
(6,65) 


CPHK  I)  .CTHETUI  tCOATll  II.  COAT  21  1),CDAT3(  1  LCD  A  I  41  II 


BISTATIC    SCATTERING 
IF     (  IB  I SC. IF .01     GO    TO 
WRITE     (6,75) 
WRITE     (6,81) 
WRITE     (6,771 
WRITE    16,82) 


54 


385 
386 
387 
388 
389 
390 
391 
392 
393 
39* 
395 
396 
397 
398 
399 
400 
401 
402 
403 
404 
405 
<i06 
407 
408 
409 
410 
411 
412 
413 

414 

415 
416 
417 
418 
419 
420 
421 
422 
423 
424 
425 
426 
427 
428 
429 
430 
431 
432 


WRI  U 
WPI  TE 
L  =  0 
INC  - 
NPL 


47 


lb, 61  I 
16, 8?) 


1 


CPHI (LI .CTHETIL) 


IF     (N8I  P.EO.  1  )    GO    TO    46 
NPHS    «    (PHSF-PHSI  1/STEPU 
NIHS    -    I  THSF-THS  I  l/STEPM 
IE     (  IWR. LE.O)     WRI TE     16,67) 
GO    TO    48 
IF    (A8I T.GT.500. )    GO    TO    47 


NPL     =     1 

NPMS    •     360 

NIHS    ■     1 

PHSI     -     0. 

THSI     "     A3  I T 

STEP    ■     1. 

GO    TO    s8 

1    NPl     -     2 

NPMS    ■      1 

NIHS    -     360 

PHSI     -     AIM 

P 

THSI     -     0. 

STEP    ■     1. 

1    PH    »     PHSI- 

STEP 

on   51    KM  . 

NPHS 

PH     =      PH  t^  I 

EP 

TH    c     THSJ- 

SUP 

IF     l( ILRD. 

GT  .0)  .AND. 

,(  (  TH.G 

00    51     1=1, 
TH    =-     TH»S1 

NIHS 

EP 

L     «     L»l 

CALL     GFFLO 

(  IA , 1 B, INC , INM, 

1  iCCGO.CG, 

2GPP.GI  I  ,PH 

CJ.CMM.O.ECSP.EC 

.SC-D.SCSP  , 

SCST.S 

3AM, ERR, IGROI 

IF     (IWM.GI 

.0)    GO    TO 

49 

IF     (NPL. LI 

.0)    W*ITE 

1 6,63) 

CPt'lUI     - 

PH 

C  IHE  ML  )     • 

IH 

COAII (LI     ■ 

EPPS 

C DA  12 (i  |    . 

EPIS 

C0AI3IL)     • 

E  IPS 

CUAT4(l  )     • 

ETIS 

IF     (NPl .NE 

.11    GO    TO 

50 

0AIY1 { •  t     • 

CABSIEPP! 

.90) . AND. 1TH.LT. 2701 ) 1     GO    TO    51 


IWR,  I   I,  I  2,  I  J,  I  I  2  ,  Md,N  ,NO~,NM,  AM.  ACSP  .ACST 


PH,TH,SPPM,SPTM,STPM,STTM 


433 
434 
435 
436 
437 
438 
439 
440 
441 
442 
443 
444 
445 
446 
447 
448 
449 
45J 
451 
452 
453 
454 
455 
456 
457 
458 
459 
460 
461 
462 
463 
464 
465 
466 
467 
468 
469 
470 
471 
472 
473 
4  74 
475 
4  76 
477 
478 
479 
480 


29 


50 


51 


52 


53 
54 


56 
57 


c  abs ( e  p  r s  J 

CABS! E IPS) 

cabsiet  t s i 
2)  in  io  51 

CABSl EPPSI 

cabsi eptsi 
cabsi e  ptsI 
cabsietts) 


OAIY2IK I  . 

DAI Y3(K  )  * 

OA1Y4IK)  « 
IF     (NF'L.NE 

DA  I y I ( | »  * 

DAI Y2(  I  I  « 

DA  I Y3I  I  I  « 

DATY4I  |  |  • 
COM  INUE 


WR  I  IF     I  6,82  I 
IF     INPl .LE .01     GO    TO    52 
CALL     PliLPRT     (3,DATYl) 
CAI  L     PC'LPRl     I4,UATY2I 
CALL     PJLPRT     I5,DATY3I 
CALL     Pf!LPRT     I6,0ATY4| 
I  f     IKFLA&I  18  I.NE.l  I 
WRITE     (6,661 


"L 


TO.   5* 


00    53 
WRI  IE 


I'LL 
16,651 


CPHMII  .CTHET1  II, COAT  Ml  I,  COAT  2  I  II,  COAT  3  I  ll.CDATAI I ) 


II FL  A&.FO.l | 
I  IFLAG.EQ.2l 


00    55    I »l, 
KFLA&III      • 

KFLA&I8  I     ■ 
KFLA&I iO) 
KFLAGI 261 
IF     I  I  FLAG, 
IF     ( IFl AG. 
GO    10    4 


10 


■     1 
=     1 

Eg. 31 

EQ.6I 


GO    10 
STOP 


WRI  TE 
MR  I  TJE 


16,68  1 
16,115) 


F  OF.  MAT 
FORMAT 


<  1H0I 
I IOX, < 


.      THF     RAOIATION    EFFICIENCY     IS     '  ,  F  1  5. 7/ / 1  OX,  '  THE     TIME-AV 
J7    i0!   .fP?^    1//V     IS     ''':15.7//10X,'THE     ANTENNA     IMPEOANCE     IS     !,F15. 


FOR 


58i  *T?fl5Toi}2?,l  T?1  .  N|*?-f  IEUO    ELECTRIC    FIELO     INTENSITY    AT     THE    OBSERV 
lATION    POINT    •fEll.5,>,'iElUSi>,,tEll.S,<     (X.Y.Z    RESPECTIVELY)    IS: 

"ix^:.!F2?x5:7:?-:^!F5iL'7/;/j;iFl5-r/20x,,EY-,Fi5-7''  ,''f"-"» 

60iF6pmaT     I3X,FS.I,2X,F5.1,3X,EI0.4,2X,EJ.0.4,2<3X,3(E10.4,2XI,F6.1,1X 
61lF?SHpH|iJ*Fl5!l??     rHESTAl%F5"//"'N&    ™E     '  NC  '  0ENT  '  '  ™l  '  '  PL*NE    WAVE 


4B1 
-.82 
483 
484 
485 
",86 
487 
488 
489 
490 
491 
492 
493 
494 
495 
496 
497 
498 
499 
500 
501 
502 
503 
504 
505 
506 
507 
508 
509 
510 
511 
512 
513 
514 
515 
516 
517 
518 
519 
520 
521 
522 
523 
524 
525 
526 
52T 
528 


62  F 
IT 
24 
3 
4* 

63  F 

64  F 
I  I 

\\ 

65  F 

66  F 
1 

2 
3H 

67  F 
ID 
2 

68  F 
10 

69  F 
IN 
2X 

70  F 
IP 

71  F 
1U 
2R 
3 

72  F 
IS 
2R 
3 

73  F 
1 

74  F 
I 

2 
3 
4 
5 
6 

75  F 
Jb    F 

77  F 
l» 

78  F 

79  F 


OR  MAI 
I NC I  I  U 
X,  31  5X 
TUl TA' 
,3I5X, 
UPMAT 
OR  MAT 
It  I  UlN 

WAVE* 
llf  TA'/ 
UF  MA  I 
[IRMAT 
MATRIX 
•  PHI -P 
I  Till 

URMA  I 
FMT-  SC 

•  THE  1 
ORMAT 
/DP    CH 
ORMAT 

,.'.!>. 

21  dX, 

I'lRMA  T 

i n» »•• 
iikma  t 
r  • 

u 

ORMAT 
D 

ORma  T 
CAl  L  */ 
ORMAT 
T50, •• 
T50,  "> 
150, •• 
T50,  •«• 
T50, •• 
T50,'» 
ORMAT 
(JHMA  T 
URNAT 

I) 
ORMA  T 
III- MA  I 


I  •  INC  10 
N'  ,  11  14, 
, 'CROSS 
.IX, 'PHI 
lPHl • , 7X 
UX,2(F5 
IT54,'B» 
SCAtTE  * 
,12  3.  'PH 
'  Pi-i| 
I  IX  ,2(F5 
(  154. 'Bl 
•/•  ftBSE 
HI • , T49, 
TA' ,4X  ,4 

I'    ,)r,ins 

AT  I  I  R  t  [)  I 
A- PHI  '  ,1 
I 1H1.5X, 
ANGES'// 
I54X,  •  El 
•  THI  T  A  •  , 
•REAI  •  ,  B 
I  IOX,  '♦• 

I40K, ■ • 
■/40X, • • 
•' /40X,' 

••  I 
I40X, '♦ 
'  /4UX , '  • 
•■/40X,' 

••  I 

I30X, • A 

/50X, 'PC 

(• l',T50 

OH 

ANTE 

M 

NAV 

• , T86.*« 
I • I ' , T50 
O50,«»« 
1150, ••' 


ENT •  . T  2 

• SCATTE 
SEC  TIlJN 
-PHI ■ , 3 

,  •  THF  tA 


1,  l»l  , 

K.S 


CAM 
NG  MA  T 
-PHI • , 
HF  I  A'  , 
I  ,  Ml  , 
T Al  1C  ' 
VAT  ION 
PHI  -IH 
3X, 'RE 
ATiON' 
/  '  PH 
8,' IHE 
CONT IN 


/,'  ECHD    A 

RlNG' /' 
' ,6X1 /' 
X, 'PHl-TH 
'.4X11 
101  E10.4, 
ER ING'/' 
Rl  X'  /  •        P 
T49, • PHI  - 
3X.4I  3X,  • 
2X, 41  Fl  I  . 
/T3/,'ELE 
,T50.'  I  I 
F  TA1  , t 76, 
AL '  ,8X,  •  I 
,  T27,  TCH 
I        THF TA* 
TA-THE  TA' 
LIE    fXECUT 


RfA    SI&MA',lt6,,ABSORPTI0N,,190,,FX 
PlANE',T25,MINCIDENl-SCAI7ERED)',l 
WAVE     •  ,52X,  3  I  IOX, 'FOR' ,  UXI /•        PhI 
f.  TA'  ,4X,  '  THE  I  A-PHl  ■  ,2X,  •  THETA-THETA 


2X11 

INCIDENT'  ,T3  7,  'ELECTRIC     FIELD     POLAR 
LANE'  ,T49,'(  INC10ENT-SCATTEREDI'/3X 
THETA',T7S,'IHE1A-PHI',T102,'THETA- 
RF AL '  ,RX  ,  '  IMAi'  ,8X11 
5, 2X. Ell . 5,3X1  I 

C7RIC     FIELD     PJl AR I 7AT IUN     SCATTERING 
NCIUtNI-SCATTFREOI'/'  POINT ' ,14X 

•THETA-PHi',TlOl,'THETA-THETA'/» 
MAG' ,8X1 ) 

O    ARLA    SIGMA'/'  PO I  NT •  , T2 5 . • (  I NC I 

T14,'PHl-PHI',T24,'PHl-IHETA',T37. 


ION  WITH  THE  FOLLOWING  ADDITIONS  AN 


P 


FCIRIC  FIELO  1NIFNSI  IYV5X,  'DEGREES' ,11X, 'POWER  GAI 
4;x,'PHI'/lX,'THtTA',3X,'PHI',fX,'rHETA,,8X,'PHI',l 
<  ,  •  I  MAC,  '.BX.'mAGN'.SX.'PHASE'II 

•  ..munR     IN    DATA    CARD    NJMBER     ',12.'        EXECUTION    STOP 

A    WIRF     SEGMENT     MAYNOT     BF     SHARED    BY    MORE     THAN    FO 
UIPOLE     MOOES CHECH    DESCRIPTION    DATA    CA 

•  EXECUTION     STOPPED 

AN  ISOLATED  WIRE  MUST  HAVE  AT  LEAST  TWO  SEGMENT 
AND  THREE  POINTS CHECK  DESCRIPTION  DATA  CA 

•  EXECUTION  STOPPED 

BACKSCAT IF Rl NG  CALL  MUST  BE  INCLUDED  FOR  A  BISIATIC 
3UCSI  IGNORIO'  /////I 
,  371 • ••  I/T50,  '•'  , I86,'*'/ 
lO  STATE  UNI VERSllY  ••/ 

NNA  ANALYSIS  PROGRAM         ••/ 
OOIF ICO  FOR  USE  AT  ••/ 

AL  POSTGRADUATE  SCHOOL       ♦'/ 
17  JULY  1974  ••/ 

• /ISO, 371' •• II 
,291 ••'I/T50,'*' .178, ••' > 
.1 IX, 'ANTENNA'  ,T  78,'*' I 
,9X,,CAlCULATlONS',T78,,»,/T50,,»,,T78,'»'/T50,29l' 


150, '••  ,9X,'NFAR  F  1EI(V  .778.  ••'  I 
150, ••• ,9X,'FAR  FIELD' , f 78, '••» 


529 
530 
531 
532 
533 
534 
535 
536 
537 
538 
539 
540 
541 
542 
543 
544 
545 
546 
547 
548 
549 
550 
551 
552 
553 
554 
555 
556 
557 
558 
559 
560 
561 
562 
563 
564 
565 
566 
567 
568 
569 
570 

512 
573 
5  74 
575 
576 


30 


80  r 

81  f 

82  f 

83  f 

8*  F 
85  F 
66  F 
8/  F 

88  F 

89  F 

90  F 

91  F 

92  F 

93  F 
9*  F 

I 


95 

96  F 
1. 
22 

97  F 
I 

98  F 

99  F 

100  F 
1 
21 

101  F 
11 
21 

102  F 
1* 
?I 

103  F 
11 

104  F 
105 
1  Uo 
10/ 
108 
109 
1  10 
HI 

112  F 

113  F 
114 
115 


T50, 
150, 
//// 
130. 
t  30, 
I  30, 
130, 
I  30, 

no, 
r  30, 

T30, 
130, 
T  30, 
I  30, 
T30, 


ORMAT 
IIHHA  1 
llflAl 

urma  r 

PR  MA  I 

URMA  T 
OrMA  I 
OPMA  T 
ORMAT 
till"*  I 
UP  MAT 
I.  P  M  A  T 
(!>.NA  r 
I  IP  MA  I 
(1PMA  T 

HI. 51 
OPXAI  (T30, 
ORMAT  (T50, 
18XI/T21, «N 
112. 5Xf Ell. 
OhmaT  (T50, 
'Rf  AL'  ,  7X,  ' 

ormat  t  r so, 

UP  HAT  (T30, 

uPHAt  t  r  30 . 

•A  NO  THETA 
50, •  IN  SIF  P 
OR  MA  I  1130, 
,'  AND  THFT 
50, • IN  STIP 
ORMAT  (T30, 
,F5. 1 , '  AND 
SU.'lN  STEP 
OPMAT  (T30, 
.5,5X1  I  I 
UPMAT  IT30, 
ORHAF  (130, 
URMAT  (T30, 
ORMAT  (  T30, 
UPMAT  (T30, 
CRMAT  1130 


*  • 


,7X, 


'BACKSCATTER ING'  ,T  78 
'BISTATIC     SCATTERIN" 


?l 


,T78, ■••I 


FREQUENCY     (MHZ )• ,T81 ,E1 1.5  I 

Hl-tE    RADIUS    (METERS)  '.TBI, ell. 51 

MIRE    CONDUCTIVITY    I MEGAHHJS/MC TERI '. TBI. Ell. 5) 

WIRE      INSULATED     INU/YES)  •  ,  (85,  '  YES' I 

WHE     INSULATED    (NG/YESI  ■  ,TB5.'N0'I 

INSULATION    RADIUS     (METERS)  •  .18  1, Ell. 

INSULATION    CONDUCTIVITY     (  MHOS/ ME  TE  R  1 

INSULATION    DIELECTRIC    CONSTANT     (RELA 

INSULATION    LOSS    T ANIENT •  , T 8  1 , E  1 1 . 5  I 

EXTERIOR     MEDIUM'  , TBI  ,'FRFE     SPACF'I 

EXTERIOR    MfcDIUM    CONDUCTIVITY     (MHOS/M 

EXTERIOR    MEDIUM    UIELECTPIC     CONSTANT 

OSS  TANGENT'  ,T81,E1 
2( 'NOD 
13X.-2 


51 

•.TBI, Ell. 51 

1  IVE)  ', TBI, F  11.5) 


•FXTER I  OR 
•Wl RE  SIR 
U.' ■ 3X.2( 
5 , 4  X  .  E  1 1 . 

•  Antenna 

IMAG 1NARY 
■•'  ,  6X,  • 
•ST  <UC  TU« 
•FAR  FIEL 
VARY ING  F 
S  OF  •  ,F5 
•BACKSCAT 
A  VARYING 
S  CF  •  ,F5 
•Bl  STATIC 
THETA  VA 
S  OF  ',F5 
•NtAR  FIE 


MED  I 
UCTUR 

•  NO. 
5.4X, 
FEEDS 
•/IT4 
OUTPU 
E  CUR 
OS  FJ 
ROM  " 

T  E  &  I :( 

FROM 

•  It1 
SCAT 

RYING 

LOS  F 


UM  L 
E'// 

•  ,  9X 
Ell. 

•  /T4 
1.12 
T  RE 
RENT 
R  PH 
,F5. 
OtGR 
G  FO 

,F 
DLC-R 
TERI 
FRO 
OEGR 
OR  F 


.•X'.HX, 
5,1X1) ) 
0,  'NODE'  , 
,6X  ,21  4X, 
QUESTED'  , 
S') 

I  VARYING 
I,'  TO 
FES.'  ) 
R  PHI  VAP. 
5.1  ,'  TO 
EES.' ) 
NG  FOR  PH 
M  •  ,F5.1, 
EES.' I 
OLLOWING 


,4X,  i 
'Y'  ,1 

16X, 
Ell. 
T  78, 

FRO 
F5.1 

Y  ING 
•,F5. 

I  VA 
•  TO 


•VOLTS 
51)1 
'•'  I 

M'.IX 
I 

FROM 
1/ 

RYING 
'.F5. 


ETERI '.TBI, Ell. 5) 
(RELATIVE)', TBI. 

1.5) 

E' ,19X, 'LOCATION' 

•  ,7X1/1 121,  12, 5X, 

*/T41,'N0.',12X, 


F5.1,'  TO  SF5.li 
•■F5.lt1  TO  «,F5. 


FROM  'tF5.lt'  TO 

1/ 


POINTS  (X,Y,Z)'/50IT40,3IE1 


F 

F 
ID 


ORMAT 
ORMA  T 
ORMAT 
ORMAT 
ORMAT 
ORMAT 


(T30, 

(//) 

(T50. 

(T50, 

(150, 

I  IOX, 


•PLOT     FOR  FAR    FIELD     THE  T  A  '  ■  ,  F  5.  1  1 

•PLOT     FOR  FAR     FIELO    PHI='.F5.1> 

'PLOT    FOR  BISTATIC    SC  AT  T  ER  I  NG- F  OR    THE 

•PLOl    FOR  BISTATIC     SCATTERING    FOR    PHI 

•PLOT    FlU  BACKSCATTERING     HETA-',F5.1 

■PLOT    FOR  BACKSCATTERING    PHI«',F5.1I 
•NO    OUTPUT    OR    PLOTS    REQUESTEO'I 


TA-',F5.1I 
-'.F5.ll 


A    CHANGE 


••• ,T78,'»'/T50.291 '••! ) 
'•*,     8X, 'INPUT    DATA     '.T78." 
29( '♦•  )/T50.'»'  ,T78,  •♦■) 
•SINCE     THIS    DATA   BLOCK   OOES 
MO     IS    ASSUMEO'I 


'I 

NOT  HAVE  A  TERMINATION  CAR 


577 
578 
579 
580 
581 
582 
583 
58* 
585 
586 
5B7 
588 
589 
590 
591 
592 
593 
594 
595 
596 
597 
5  98 
599 
600 
601 
fr02 
603 
604 
605 
606 
607 
608 
609 
610 
611 
612 
613 
614 
615 
616 
617 
618 
619 
620 
621 
622 
623 
624 


116  FOPMAI 

I     STRUCT 
2»    FXECU 

117  FCPMAT 

1     UNI)    PL 

118  FORMAT 

I     ,'REAL 
1  19    FOPMAT 
IX, 'Pi  AL 

120  FOPMAT 
IX, -PEAL 

121  FORMAT 
lS'/r30, 
2./T30.I 

122  FORMAt 

123  FOPMAT 
12*    FOPMAT 

I     T30,'G 

125  FORMAI 

126  FORMAT 

127  FORMAT 

!••• 1 1 
END 


I//10X.40I • •• 1 /IOX, ' 
URE •/ IOX ,• MUST  BE  SI 
T ION    STOPPED    ••♦ • I 

(//UXj'NO  PART  OF  THE  WIRE  STRJ 
ANE . ' / IOX , • ••••EXE  CUT  I JN  STOPPED 
(  T50,  '  STRUCTUKF  LOADS'  /TW,  •  NODc 
■,7X, • IMAGINARY' /( T4 I, 12 ,6X, 21 4X 
I T  50. '  STRUCTURE    LQAOS'/I39,'SEGM 

,,/X,'IMA0lNARY,/IT41.12,6X,2(4X 

( t 50, 'AN  IF  UNA  FEEDS' /139,'SEUMEN 
•,7X,'IMAG1NARY'/IT4|,|2,6X,2(4X 
(//T30,'THE  NUMBER  OF  INTERVALS 
•IN  THE  IMPEDANCE  MATRIX  WITH  SI 
3,'.  IF  CLOSED  FORM  INTEGRATION 
(T30,'GR0UN0  PLANE  ( NU/ Y E S  )  ' , T85 
TBOi'GROUNO  PLANE  I NO/Vt S) ' . TB5 
1T30, "GROUND  DIELECTRIC  CONSTANT 
ROUND  CONDUCTIVITY  (MHOS/METER)' 
(T30,'GR0UN0  PLANE' ,T83, 'PERFECT 
(T30. 'ANTENNA  HEIGHT  (METERS)', T 
(//IOX. 401 '••I/IOX,' THE    HIRE    RAD 


THE     DESCRIPTION    AND    THE     GEOMETRY    OF     THE 
ATED     IN    THE     FIRST    DATA     BLOC  K .'/  IOX  ,« »•• 


CTURE    CAN    LIE    BELOW    THE    GRO 
•*»»>) 

•  ,  16X,' OHMS' /T41,  'NO.  •  ,12X 
,E1  I  .5)  ) I 

ENT,,UX,,0HMS'/T41,'NO',12 
, Ell  .5)1  I 

t' ,14X,'V0LTS'/T41,'N0.',IZ 
.611  .511  I 

FOR    CALCULATING    THE     ELEMENT 

MPSONS-RULE     INTEGRATION    IS' 

IS    REQUIRED    SET    INT-O'///) 

,  • NU  '  I 


VE  I  ',  TBl.f  11.5/ 


•YES'  I 

(RELAM.. 
•  TBI  , Ell. 51 
1  ) 

81. Ell. 5) 
IUS    MUST    BE    STATE0'/10X,40( 


625 
626 
627 
628 
629 
630 
631 
632 
6  33 
634 
635 
636 
637 
638 
6  39 
640 
641 
642 
643 
644 
645 
646 
647 
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BLNK 


PURPOSE:       to  compress  data  to  the  left  by  removal  of 
the  blank  spaces  on  the  input  data  cards. 

METHOD:         MI)  character  is  compared  to  the  blank;  and 
if  it  is  true,  the  A  (1  +  1)  character  is  shifted  to  the  A  (I) 
position. 


CALLED  BY, 


READ 


CALLS  TO 


NONE 


simnnui  ine  ulnk  i a) 

OlMtNSlON  M60) 
UATA  BLANK/'  •/ 
K    •    0 

00     I     I "I, BO 
J    ■     l-K 
At Jl     •     At  I  I 
1     IF     U<  II.  60.  BLANK)     K-KH 

IF     IK.EO.OI     RETURN 
AI8I-KI     •     BLANK 
RETURN 
ENO 
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CBES 


PURPOSE:        to  calculate  the  quantity  E01  where 

B01  =  J  (z)  /  J  (z)  . 
o        1 


METHOD:         If  the  absolute  value  of  the  argument  for  the 
Bessel  function  is  less  than  12,  B01  is  calculated  via  the 
power  series  expansion  for  the  Bessel  function  in  the  DO 
LOOP  ending  at  statement  3.  If  greater  than  12,  the 
asymptotic  expression  is  utilized  at  statement  4.  If  the 
magnitude  of  the  complex  part  of  the  argument  for  the  Bessel 
function  is  greater  than  20,  B01  is  set  to  (0.,-1) .   If  the 
complex  part  of  the  argument  is  negative,  the  sign  of  E0 1  is 
changed  prior  to  returning  to  the  calling  program. 

CALLED  BY:      SGANT 
CALLS  TO:        NONE 


SUMROUT  INE    CBES     U.B01  I  1 

COMPLEX    ARG.CC.CS.EX  2 

COMPLEX     BOI ,lt TfcRMJ,TERMN,M22*, JN(2)  3 

DATA    Pi  /3.  l<>159/  * 

IF     I  CAR  SI  I  I  .GE  .  1  2  .01     GO    TO    4  S 

FACTOR    »    0.0  6 

TEkMN    =     (0. ,0. I  7 

W/24    '     -0.  25«2«Z  B 

1ERHJ    «     I  1  .0,0.01  9 

10 

DO   3   NP  =1,2  II 

N    =    NP-1  12 

JN1NPI     -     TEBHJ  13 

M    •    0  1* 

M    =    M»  1  15 

TERMJ    *     TERMJ«MZ2<./FLOATlM»(N»MI  )  16 

JNINPI     *    JNINPItTERNJ  17 

IF     INP.NE. I  I     CO     TO    2  18 

FACTOR     =    FACT03 » 1 .0/FL0ATIM1  19 

TEHMN    ■     HPMN»TERMJ»FACT0R  20 

ERROR    =     CABSITFRM.J)  21 

IT     (  ERROR. GT.  1.0E-10I     GO    TO     1  22 

TEBHJ    =    0.5»Z  23 

2* 

801    •    JMII/JNI21  25 

RETURN  26 

Y    ■    A1MAG12)  27 

IF    IABSIYI .cr.20.1    GO    TO    5  28 

ARG    ■     (  .0,  1  .  I  'I  29 

EX    «    CEXPIAR&I  30 

CC    -    EX»1  ./EX  31 

CS    ■    (  .0,-l.l»<EX-l./EX)  32 

BOI    -    ICS»CCI/(CS-CCI                       .  33 

RETURN  3* 

BOI    •    (.0.-1. I  35 

IF     (r. IT. 0.1     BOI     •     1.0,1.)  36 

RETURN  37 

END  38 
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DSHELL 

PURPOSE:        to  calculate  the  mutual  impedance  term 
contributed  by  the  dielectric  insulation  on  the  surface  of  a 
thin  wire. 

METHOD:         The  contribution  to  the  impedance  matrix  is 
calculated  utilizing  the  equation  below 


(8   -  5)  ln(b/a) 

z    =  -  .1-2 . _      ?•  (1)  F'(D  <il  , 

mn       2n  jw&S  J  m      n 

2         m,  n 


J  m 


where  z    is  defined  in  subroutine  SGANT,  &      is  the 
mn  2 

dielectric  constant  of  the  insulation,  b  is  the  outer  radius 

of  the  insulation,  a  is  the  inner  radius,  S  is  dielectric 

constant  of  the  external  medium,  and  F  is  the  sinusoidal 

expansion  function. 

CALLED  BY:      SGANT 
CALLS  TO:        NONE 


SUBROUTINE    DSHELL    ( AM, BH, DK.CCOS . SGOSi EP2 . EP.ET A ,C*N |P 1 1 iP 1 2 )  1 

COMPLEX    CGDS,SGOS,EP2,EP,Erl,GAM,PU,P12  ,GO,CST  5 

DAIA    PI /3.1*l5V/  3 

GO    -    GAM*OK  t, 

CSI    -    ( EP2-EP)»ETA»ALOG(BM/AMI/(*.»PI»EP2»SGOS*SGDSI  5 

Pll    ■    -CSr»(GO«SGOS«CGOS)  6 

P12    •    CST»(G0»CG0S»SGOS)  7 

RETURN  B 

ENO  9 
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EQUAL 


PURPOSE:       to  determine  position  (location)  of  the  equal 
symbol  on  input  data  card. 

METHOD:        The  character  search  begins  in  the  column 

passed  to  the  subroutine.   On  returning  to  the  calling 

program,  the  argument  passed  is  the  column  following  the 
equal  symbol. 


CALLED  BY: 


READ 


CALLS  TO: 


NONE 


SUBROUTINE  EQUAL  IN) 
COHHQN  /If    A180) 
DATA  EQUAL/'-1/ 
K  •  N 

00  I  I-K.80 
N  -  IM 

IF  (6(1). EQ. EQUAL)  GO  TO  2 

1  CONTINUE 

N  ■  I 

2  RETURN 
ENO 


10 

ii 

13 
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EXPJ 

PURPOSE:       to  calculate  the  exponential  integral  with 
complex  limits. 

METHOD:        The  exponential  integral  is  defined  as: 

V2 


./■ 


_e dv  =  E  (V1)  -  E  (V2)  +  j2nn  , 


where  the  integration  path  is  the  straight  line  from  V 1  to 
V2  on  the  complex  v  plane  and 

.00 

E  (z)  =    f  e   dt  . 


The  integration  path  is  a  horizontal  line  in  the  w  plane 
or  an  inclined  straight  line  from  V1  to  V2  the  v  plane.   The 
integer  n  is  zero  unless  this  path  intersects  the  negative 
real  v  axis  at  a  point  between  V1  and  V2.   When  there  is 
such  an  intersection, 

a)  n  =  1  if  Im(V1)  >  Im  (V2) 

b)  n  =  -  1  if  Im  (V1)  <  Im(V2)  . 

The  term  j2mi  is  caluclated  below  statement  12. 

CALLED  BY:      GGMM 
CALLS  TO:       NONE 
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SUIiKOUTINE     EXPJ     (VI.V2.W12)  ) 

COMPLEX    FC.bl5.St I.UC.VC.V1 .V2.W12.Z  2 

OIMFNSKIN    VI21I,     W<2li,     1)116),     E(lo)  3 

0*  14    V/0. 2228*66  7F00.0. 1 188932 1 F01 ,0.2992 7363E01  ,0.57751*36E01,0.9  * 

18  3?*<>  r*E0t,J.l59B2df*EO2,0.933078l2E-0l,O.*9269l7*E00,0.12l5595*EO  5 

21. 0.226 99*95 EOl.O.  366  762  2  7E0I, 0. 5*253 36 6E01, 0.756591 62 E01.0. 10120 2  6 

32  8EO2,0.l3l3O232EO2,0.166  5**O8E02,O.2O7  76*7  9E02.O.2562389*E02,O.31  7 

**0/519t02.0. 385306 83E 02, 0.*802608  6E 02/  8 

0414    W/0.*5B96*6OE0O,0.*l  70008  3  E  00,  0.  1  1  3  3  73  38FOO  ,0  .  1  0399  19  7E  -01  ,0.  9 

126101  72  OE  -0  3,  0.8  9  85",  791  E-06,0.  2  1  8  23*8  7E00  ,0  .  3*22  I  01  7E00,  0.26302  758  10 

2E  00,  0.  1  26  *2  5  82E  00,0.*  020686  5E  -01,  0.B5638778E-02.  0.1212*  36  IE  -0  2,0.  I  11 

3  1  167**0E-03,O.6*59926  7E-O5,0.2  226  316  9E-O6,O.*227*30*E-O8,O.392l89  7  12 

*JE-l0,u.l*56  5l52E-12,O.l*B3O2  7OE-l5,O.16  065  9*9E-19/  13 

04 T A    0/ 0.22* 958*2  6 02, 0.7** 11 568 E 02, -  0 . *  1 *3 1 5 76E03 , -0 . 78 75*3 39E02 ,0  I* 

1.  1125* /**E 02,0. 16021 76  I E 03, -0. 2 3 86 2 1 95E 0 3. -0.5009*6 8 7E 03, -0.6 8*8 78  15 

25*1  02. 0.1225*  77  8E  02, -0.10161  97  6E02,-0.*7219  59  1E01  ,  0.  79  7  2  9681  E01  ,-0  16 

3. 210695  7*E0  2,0.2  2  0*6*90EOI ,0.897282**E01/  17 

04 T4    E/0.2lli3l07E02,-0.3  7959  78  7E03,-0.97*892  20E02,0.l2  9006  7  2E03,0  IB 

1.1  79*9 226EO^,-0.12910931E03, -0.55705 5 7*E03,0. 1 352 *801E 02 ,0 . 1 *696  72  19 

21EO3.0. 179*9 5 2 8E 02 ,-0.3  29  8 101 *E 00 ,0.3 10288 J6E02, 0.8 165 765 7E0 1,0.22  20 

3  2  3696  1E02,0.3912*89  2E02,0.81636  799E01/  21 

2    =    VI  22 

23 

00    12    JIMM.2  2* 

X  »  RE4L1  Z)  25 
Y  =  AIMAG(Z)  26 
E15  «  1.0..0I  27 
48  -  CAUSlZI  28 
IF  (4B.E0.0.)  GO  10  11  29 
IF  IX.GE.O. .AN0.AB.GT.10.)  GO  TO  10  30 
YA  -  ABSIYI  31 
IF  (X.LE.0..4N0.YA.GT.10. )  GO  10  10  32 
IF  <YA-X.GE.17.5.0R.YA.GE.6.5.0R.X*YA.GE.5.5.0A.X.GE.3.I  GO  TO  2  33 
IF  (X.LE.-9.1  GO  TO  6  3* 
IF  (YA-X.GE.2.5)  GO  TO  7  35 
IF  <X»YA.GE.1.5I  GO  TO  3  36 
N  =  6.»3.»AB  37 
E 15    -    1 ./IN-I.I-Z/N»»2  38 

1    N    ■    N-l  39 

E15    •    I ./IN-l.)-Z»E15/N  *0 

IF    IN.GE.3I    GO    TO    1  *1 

E  15   •    Z»El5-CMPLX(.5772l6tALOG(ABI,ATAN2<Y,XM  *2 

GO    TO    11  *3 

.2    Jl    -    1  ** 

J2   -    6  45 

GO    TO    *  *6 

3    Jl    -    7  *7 

J2   •    21  48 

*    S    *    I .0..0I  *9 

YS    »    Y«Y  50 

51 

DO    5    l=Jl,J2  52 

XI  -  VI  I  I  »x  53 
CF    *    Wl  I  l/(X  l*XI *YS)  5* 

5  S    =    S»CHPLXIXI»CF,-YA*CF>  55 

56 

GO    10    9  57 

6  T3  =  X»X-Y»Y  58 
T*  »  2. »X»YA  59 
T5  =  X«  T3-Y4»F*  60 
T6  =  X*  T*tYA«T3  61 
UC    -    CMPI XIOI 11 ) »0I I2)*X»0I 13)»T3»T5-E( 12I»YA-E( I3)*T*,E(11 )»E(I2I  62 

1*X»EI 13I*T3«T6»DI 12I*YA»DJ 13I«T*)  63 

VL    =    CiPlXI  D(  1*)  »0(  15  l»X*DI  16I«T3»T5-EI  151*YA-E(  16I*T*,E(  14I»E(15)  6* 

1*X»E( 16 )*T3»T6»0{ 15>*YA+0(  16)»T*I  65 

GO    TU    8  66 

7  T  3  =  X»X-V»Y  67 
I*  =  2.*X»YA  68 
15  '  X«  T3-YA«  T*  69 
T6  *  X«T*«Y«»T3  70 
H  ■  *•  I5-YAM6  71 
T8  =  X»  T6fYA»T5  72 
19  *  X»  I  7-Y4»T8  73 
T 10  =  X«T8»Y4»T I  7* 
UC    =    CMPLXI01 1) »D(2l»X»0(3)*T3*0(*l*T5»0(5)»f7*T9-(E(2)*YA»EI 3I*T*  75 

l»E(*l»I6»El5l»T8),E<ll«E(2)»X»E(3l»T3*E(*>»T5»E(5>*T7»T10MDl2l»YA  76 

2»0( 3l»T*»D(*)»T6»OI5l*I8) )  77 

VC    =    CMPLXID(6I »0I 7)«X»D( 8I»T3»0(9I»T5»0( 10  I  * T 7»T9-l E ( 7  I »YA» E ( 81 »T  78 

I4»C(9IM6»E< l0)«r8),E(6)»E(7l»X*E(8l»T3»E(9l«T5»E(10l»T7»T10»ID(7)  79 

2»YA»D(8)»r*»D(9)*T6»0H0)»T8))  80 

8  EC  =■  UC/VC  81 
S     =    EC/CHPLXIX.YA)  82 

9  EX  «  EXP1-X)  83 
T  »  FX»CMPLX(COS(YA) ,-SINIYA) I  8* 
EI5  ■  S»l  85 
IF  IY.IT.0.1  EI5  •  C0NJGIE15I  86 
GO    TO    11  87 

10  E15  -  .409319/1  7 ».1930**I». *2  183  1 /(2  »  1 .  02666  )  ». l*7126/IZ»2. 56788)»  88 
l.206  335E-|/(Z»*.9  00  35l». I 0  7401E -  2/1 1* 8. 1 821 51 ». 1 58654E-4/I Z» 12. 734  89 
22l».317031E-7/(Z»19. 39571  90 

E15    -     E15»CEXP(-Z)  91 

11  IF  <JIM.EO.ll  W12  •  E15  92 

12  Z  •  V2  93 

9* 

Z    -    V2/VI  95 

TH    -    A1AN2IAIMAGIZI,REAL(Z) ) -AT AN2 ( A I  HAG < V2 I , RE AL < V2 > I ♦ AT AN2 I Al  MAG  96 

I  (VI  I.PE  4L (VI  I)  „, 

AB    •    ABM TH)  'J 

IF    (4B.IT.1.  I     TH    •     .0  22 

|03 
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GANT1 


PURPOSE:       to  consider  the  wire  structure  as  a 
transmitting  antenna  and  calculate  the  input  impedance  and 
current  distribution. 

METHOD:        If  a  wire  antenna  is  driven  by  a  voltage 

generator  v  located  at  one  of  the  current  sampling  points 
i 

1   and  if  displacement  currents  are  neglected,  Ampere's  law 
i 


yields 


V   =  v   F  (1  ) 

m     i   m   i 


where  F  is  the  sinusoidal  expansion  function.  Thus,  the 

excitation  voltages  V   will  vanish  everywhere  except  where 

m 

v   is  not  zero, 
i 


The  DO  LOOP  ending  with  statement  50  uses  the  delta-gap 
model  defined  above  to  determine  the  excitation  voltage 
CJ(I)  for  all  the  dipole  modes.   These  are  stored 
temporarily  in  CG(I).  Then  subroutine  SQROT  is  called  to 
obtain  a  solution  of  the  simultaneous  linear  eguations. 
SQROT  stores  the  solution  (the  loop  currents)  in  CJ  (I)  . 

In  the  DO  LOOP  ending  at  statement  80,  the  complex  power 
input  and  input  impedance  (s)  are  calculated.  The 
time-average  power  input  (PIN)  is  the  real  part  of  the 
complex  power  input. 

Subroutine  RITE  is  called  to  make  the  transformation 
from  the  loop  currents  to  the  branch  currents.  If  IWR  is  a 
positive  integer,  RITE  will  write  out  the  list  of  branch 
currents. 
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Finally,  GANT1  calculates  the  radiation  efficiency  by 
calling  subrouinte  GDISS  to  obtain  the  time-average  power 
dissipated  in  the  lumped  loads  and  the  imperfectly 
conducting  wire. 

CALLED  BY:      MAIN 

CALLS  TO:       GDISS 

RITE 
SQROT 


SUBROUTINE     GANT  1     (  I  A , I B , I NM  ,  IWR , I  1 , I  2 , I  3 , I  I  2 . JA ■ JB . MOi N , NO.NM, AH,C  1 

1  ,CJ,CG,CHM.D,EFF,GAM,GG,CGO,SG0,VG.YlI,ni,/.I.D,is.lGRDI  2 

COMPl  EX    YY  3 

COMPLEX    C<ll,CJIll.CGO(ll,SGD(ll,VG(l),ZlO(l>  .Y II ,  II 1 , 1  S  ,GAM,CG(  1  )  * 

DIMENSION  0(11.     IaId.     IBll),    JAlll,    JBdl  5 

DIMENSION     lllll,     12(11.     1311),     MDKNH.4I,    N0(1)  6 


DO    3    I'l.N  q 

CJ(  I  )     •    I .0..01  10 

K    ■    JAI  II 


c 


[\ 


c 
c 

c  13 

00  2  KK=l,2  It 
KA  "  UIKl  15 
KB  •  IBIKI  6 
JJ  -  K  [I 
FI  -  |.  8 
IF  1KB. EO. 12(111  GO  TO  1  19 
IF  (KB.  tQ. II  I  II  I  FI  — 1.  20 
CJl I)  •  CJI I)»FI*VG(JJ|  21 
GO    TO    2  7? 

1  IF  (KA.  EQ.I3I  II  I  FI  — |.  2| 
JJ  •  K»NM  ?<. 
CJI  II    ■    CJI  ll*FI«VG(JJI  25 

2  K  -  JBI  I)  26 
C  27 
C  28 

3  CONTINUE  29 
C  30 
C  31 

C  33 

DO  *  I'l.N  34 

*  COlll  >  CJI  II  35 

C  36 

CALL  SQROT  (dCJiO.I12.NI  38 

1  12  «  2  39 
Yll  -  I .0, .01  40 
NNN  >  N  41 
IF  (IGRO.GT.OI  NNN  »  N/2  42 


•>i 


C  a 

DO  6  I ■ l.NNN  45 

NN  ■  IAIJBI  I  I  I  46 

YY  -  CJI  I l»CONJG(CG( III  47 

IF  ICABSI YYI .L(. 1.E-20I  GO  TO  5  48 

I  I  1     ■     I  ./YY  49 

-i    I  l  I      I  '■  ,  n  )     NN,  M  1  50 

5  Yll     '     Yl l«YY  51 

6  CONTINUE  52 
C  53 
C  '                 .                            5* 

IF     I  IMK.GT .01     WRITE     (6.71  55 

CALL    RITE     I  I  A.  IB,  INM,  IWR,  I  1,  12,  I  3 , MO .NO , NH , C J ,CG ,  IGROI  56 

GO    -    REAL(Yll)  57 

7  I  I    -     l/Yll  58 

PIN    •    GO  59 

CALL    GDISS     (AM,CG.CMM,0,OISS,GAM,NM,SGO,  7.L0.2SI  60 

PHAO    -    PIN-OISS  61 

EFF    «     1  00.  .PRAfj/P  IN  62 

RETURN  63 

C  64 

C  65 

7  FOPMAT     I50X, 'ANTENNA     BRANCH    CURRENTS'I  66 

8  FORMAT     IIOX.'THE     INPUT     IMPEDANCE     AT    NOOE     '.IS,'     IS'.FIS.T,*     ♦    J',  67 
IFI5.7//I  68 

ENO  69 
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GDISS 

PURPOSE:       to  calculate  the  time-average  power 
dissipated  in  the  imperfectly  conducting  wire  and  in  the 
lumped  loads. 

METHOD:  The  time-average  power  dissipated  by  the  wire 
is  calculated  in  the  DO  LOOP  ending  at  statement  1  utilizing 
the  equation  below: 


«    V     * 

P   =    s    /   I  I  dl 
d    2na   / 

0 


where  R   is  the  surface  resistance  of  the  wire  and  a  is  the 
s 

radius  of  the  wire. 


The  power  dissipated  by  the  lumped  loads  is  calculated 
by  the  DO  LOOP  ending  at  statement  3.  If  the  wire  is 
perfectly  conducting,  CMM  <0,  the  first  calculation  is 
by-passed. 

CALLED  BY:      GANT1 
CALLS  TO:       NONE 


40 


SUPPIjUTINE    GUISS     <AM,CG,CMM,D,DISS,GAM,NM,SGO.ZLDtISI 

COMPLEX    CGI1  I.SGDIM  ,ZL0(1I  tCJA,CJB,GAH,Z.S 

DIMENSION    01 1 1 

DMA    PI  II.  I<>159/ 

DISS'     .0 

If    (CMM.LE.O.I    GO    TO    2 

»LPH    *     BtAKGAMI 

BETA    •     A  IHAGIGAH I 

RH    •    REAL  I ZSI/I 4.»PI*AM) 

DO     I    K-l.NM 

DK    -    Dl M 

DEN    «    CA13SISGDIK I |»»2 

EAD    •     EXPI  AL  PM»UK  I 

CAO    ■     ( EAD»1 ./EADI/2. 

CBO    «     COSIBETA'DKI 

SAU    "     OK 

If  I ALPH.NE.O. I  SAO-I  EAO-l ./EADI/ 12.»ALPH> 

SBD  ■  DK 

IF  (8ETA.NE.0.)  S8D=S I N< BE T A*DK 1 /BETA 

FA  «  RH«( SAD»CA0-SBD»CB0I/DEN 

FB  *  2.»RH«ICAD»SBD-SA0«CBD)/DEN 

C JA  «  C  GIKI 

L  "  K»NM 

CJB  '  CG(L) 

1  DISS  -  01 SS»FA»<CABS1CJA)»*2»CABS<CJB)»»2I»FB*(REAL<CJA>»REALICJBI 
1*AIMAGILJA)»AIMAG(CJB> ) 

2  00  3  J'l ,NM 
K  -  J»NM 

3  OISS  •  DISS»REALULD(  Jl  I*  I  CABS  I  CGI  Jl  )««2  UREAL  WIDl  K  1 1*  I  CABSI  CGIK  I 
1)*»2) 

RETURN 
END 


I 
2 
3 

* 

5 

6 

7 

B 

9 

10 

II 

12 

13 

It 

15 

1? 
18 
19 
20 
21 
22 
23 
2* 
25 
26 
27 
28 
29 
30 
31 
32 
33 
3* 
35 
36 
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GFF 


PURPOSE:       to  calculate  the  far-zone  field  of  a 
sinusoidal  electric  monopole. 

METHOD:        If  an  electric  line  source  has  length  d  and 

endpoints  at  (x  ,y  ,z   )    and  (x  ,y  ,z  ),  then  the  coordinates 

111         2   2   2 

of  any  point  on  the  source  are 

x  =  x   +1  cosx 
1 

y  =  y   +  1  cosy 
1 

z  =  y   ♦  1  cosz 
1 

where  cosx,  cosy,  cosz  are  the  direction  cosines  of  the  1 

axis,  and  1  is  the  distance  along  the  source  measured  from 

the  endpoint  (x  ,y  ,z  ).   Let  the  current  distribution  on 

11   1 

the  monopole  be 


I   sinh  /(d  -  1)  +  I   sinh  yi 

1(1)  =  I ' . - ? 

sinh  /d 

where  I   and  I   are  the  endpoint  currents.  The  far-zone 
1       2 

field  of  this  source  is 

E   =  (cosx  cos6  cost  ~  cosy  cosG  sin<>  r  cosz  sinQ)  E 

E   =  (-cosx  sin$  +  cosy  cos$)  E 
Q  1 

where 

-/r 
r?e  Kgd  yf<D 

E  =   —  [  (e    -  g  sinh  /d  -  cosh  yd)     I  e 

1   Unr  (1-g2)  sinh  yd.  i 

yd  Kf<2) 

+  (e    +  g  sinh  yd  -  cosh  yd)    I   e      ] 
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fCij  =  x  sin©  cos*  +  y   sin©  sin*  +  z   cos9 
1  1  1 


f<2)  =  x  sin©  cos*  +  y   sin9  sin*  +  z   cos© 
2  2  2 

g  =  cosx  sinQ  cos*  +  cosy  sin©  sin*  +  cosz  cos9 


and  (r#9,*)  are  the  spherical  coordinates  of  the  observation 
point. 

In  this  subroutine  the  range  dependence  has  been 
suppressed.  The  far  field  vanishes  in  the  endfire  direction 
where  GK  =  0.  If  a  ground  plane  is  present  (IGRD  >  0)  the  E 

equation  above  is  decomposed  into  the  x,    yr  and  z  components 

and  the  reflection  coefficients  are  applied  before  E   and  E 

g      <|> 

field  components  are  returned  to  the  calling  program. 
CALLED  BY:      GFFLD 


CALLS  TO:       NONE 
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SUBROUTINE    &FF     (  XA.V4,  /A.  X B , Y B , J B , U , C GO  .  SGD , C TH , S T H, CPH , SPH , G AM. E I  1 

1A.ET1  ,f  12  ,EP1  ,EP2,  IgRD.ERRI  2 

CCMPLEX     ERR.RV, RH.RR.EX.EY, EZ.EE  J 

COMPLEX     M I ,ET J,EP1,EP2,GAM,ETA  4 

COMPLEX    r.O,CGD,SGO,£GD  5 

COMPLEX    EGF A,EGFB.EGGO,ESA,ESB  6 

COMPLEX    CST  7 

f P    ■     12.56657  8 

XAB    =    XB-XA  9 

YAB    •     YB-YA  10 

2AB    -     JB-/A  11 

CA    »     XAB/O  12 

CB    -    YAB/O  13 

CG    ■    ZAB/O  1* 

G    «     ICA«CPH»CB»SPH)*STH»CG*CTH  15 

GK    »     l.-G»f.  16 

Ell    <     I.O..OI  17 

E  T2    -    (  .0,  .01  18 

EPl     -    ( .0, .01  19 

EP2    =    1 -0..01  20 

IF     IGK..LT  .  .001  I     GO    TO    3  21 

FA    *     IXA«CPH»YA»SPMI«STHWA»CTH  22 

F8    -    IXB»CPHf YB»5PH)»SIH»2B»CTH  23 

EGFA    •    CEXP(GAM»FA|  2* 

EGFB    -    CEXPtGAM'FBl  25 

ECGU    «    CEXPI GAM«G»DI  26 

CSI    =    EIA/(GK»SGD»FPI  27 

ESA    =    CSI«F<".FAMEGGD-G*SGD-CGDI  28 

ESB   ■=    CST»EGFB»( l./EGGD»G*SGO-CGDI  29 

IF    (IGRO.LE.OI    GO    TO    2  30 

RV    ■    (- 1. ,01  31 

RH   '    i -  1 . ,91  32 

IE    (IGRO.EO.l)    GO    TO    1  33 

RR    ■    CSCRT (ERR-STH»STH)  3* 

RV    -    -( ERR»CTH-RRI/I ERR»CTH»RR)  35 

RH    ■     (CTH-RR)/tCTH»RRI  36 

1     EX    •    CA»£SA  37 

EY    -    CB«ESA                                                                                        .  38 

il    ■    CC'ESA  39 

EE    •    < tX»SPH-EY*CPHI»IRH-RVI  40 

EX    ■    EX»RV»EE»SPH  41 

EY   ■    EY»RV-EE»CPH  42 

11    -   -EZ»RV  43 

ESA   «    EX»CA»EY»CB»E2»CG  44 

EX   ■    C»»ESB  45 

EY    -    CB»ESB  46 

E2    •   CG»ESB  47 

EE    -     IEX»SPH-EY»CPH)»IRH-RVI  48 


EX    .    EX»RV»EE»SPH 

EY   ■    EY»RV-EE»CPH 

E2    •    -E7»RV 

ESB   •    EX«CA»EY»CB»E2»CG 

Z    I    '    (JA,SPH,CB»SPHI»CTH-CG»STH 

P    -    -CA»SPH»CB»CPH 

ET1    -    T»ESA 

ET2    -    T»ESB 

EPl    •    P»ESA 

E P2    -    P»ES8 
3    CONTINUE 

RETURN 

END 


49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
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GFFLD 


PURPOSE:       to  calculate  the  far-field  for  the  thin  wire 
structure. 

METHOD:        The  far-field  for  the  structure  is  calculated 
from  the  loop  currents.  The  loop  currents  are  either  the 
currents  produced  by  the  transmitting  antenna  calculations 
of  subroutine  GANT1  or  the  currents  produced  by  an  incident 
plane  wave. 

If  the  incident  field  is  generated  by  a  distance  source 

with  spherical  coordinates  (r  ,9  ,*  ).    the  excitation 

0   0   0 

voltages  induced  by  a  incident  plane  wave  are 
V   =  /  F  E   dl 


1      =  f F  E 
m    /mi 


m 
where 


E   =  E   exp (Y   r  f   )  dl 
i     0  0 


where  E   is  a  vector  constant,  r   is  a  vector  from  the 
0  0 

coordinate  origin  to  the  distance  source,  and  r  is  the 

radial  vector  from  the  origin  to  the  observation  point. 


The  field  E   is  generated  by  test  dipole  m  when 

M 

radiating  in  the  homogeneous  medium.  Using  the  vector 

potential,  the  field  at  the  distance  point  (r  ,0  ,  <J>  )  is 

0   0   0 


/' 


]wu  e    ° 

E   =  -  /F   exp(/  f  7    )  dl 

m      4nr        /  a  ° 
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where  the  radial  component  is  to  be  suppressed.  From  the 
above  equations, 

4nr    /r 

V   = 2-  e   o  e  e 

m      jwu         o   m 

If  an  antenna  gain  calculation  is  desired,  INC  is  set  to 
zero.   PH  and  TH  denote  the  spherical  coordinate  direction 
of  the  distance  observation  point.  The  phi-polarized  (EPPS) 
and  the  theta-polarized  (ETTS)  components  of  the  electric 
field  intensity  are  returned  to  the  calling  program. 

If  INC  =  1,  a  backscattering  calculation  is  desired.  In 
this  case  PH  and  TH  denotes  the  incident  angles  for  the 
incident  plane  wave.  These  are  also  the  spherical 
coordinates  of  the  distance  source.  The  outputs  returned  to 
the  calling  program  include  absorption,  extinction,  and 
scattering  cross  section  for  each  polarization;  scattered 
electric  field;  and  echo  areas. 

If  INC  =  2,  a  bistatic  calculation  is  desired.  In  this 
case  PH  and  TH  denote  the  spherical  coordinate  of  a  distance 
observer.  Since  this  calculation  uses  the  induced  loop 
currents  (EP  and  ET) ,  a  backscattering  call  must  preceed 
this  calculation.  The  outputs  returned  to  the  calling 
program  consist  of  the  scattered  electric  field  components 
and  echo  areas. 

EPP(I)  and  ETT  (I)  denote  the  phi-polarized  and 
theta-polarized  far-zone  fields  of  dipole  mode  I  with  unit 
terminal  current.  In  a  backscattering  situation,  the 
excitation  voltages  EP(I)  and  ET(I)  are  obtained  by 
multiplying  EPP  and  ETT  by  the  constant  CJI.  Then  calls  are 
made  to  SQHOT  which  stores  the  solution  (the  induced  loop 
currents)  in  EP(I)  and  ET(I).  RITE  is  called  for  the  branch 
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currents  CG ( J) ,  and  GDISS  is  called  for  the  time-average 
power  dissipated  in  the  imperfectly  conducting  wire  and  the 
lumped  loads.  This  power  is  denoted  PDISS  and  TDISS  for 
phi-polarized  and  theta-polarized  incident  waves, 
respectively. 

In  scattering  problems,  the  incident  plane  wave  has  unit 
electric  field  intensity  at  the  origin.  GGG  denotes  the 
time-average  power  density  of  the  incident  wave  at  the 
origin.  ACSP  and  ACST  denote  the  absorption  cross  sections 
for  the  phi  and  theta  polarizations. 

PIN  and  TIN  denote  the  time-^average  power  input  to  the 
wire  structure,  delivered  by  the  eguivalent  voltage 
generators  VP  and  VT  at  the  terminals.  PIN  and  TIN  apply  for 
the  phi  and  theta  polarizations,  respectivity.  The 
time-average  power  input  is  regarded  as  the  sum  of  the 
time-average  power  dissipated  and  the  time-average  power 
radiated  or  scattered  by  the  wire.  ECSP  and  ECST  denote  the 
extinction  cross  sections  and  SCSP  and  SCST  denote  the 
scattering  cross  sections. 

The  distance  field  is  calculated  in  the  DO  LOOP  ending 
with  statement  7  for  scattering  situations,  and  in  the  DO 
LOOP  ending  with  statement  9  for  the  antenna  situation. 

The  radar  cross  sections  (echo  areas)  SPPM,  SPTM,  STPM, 
and  STTM,  are  defined  as 

2   2ar 
a   =  limit   4nr   e     S   /  S 
r  si 

where  S   and  S   denote  the  time-average  power  densities  in 
s      i 

the  scattered  and  incident  fields  evaluated  at  the  origin. 

For  an  antenna,  the  following  definition  is  employed  for 
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the  power  gains: 

2   2ar 
G  («,♦)  =  limit   4nr   e    S(r,0,<t>)  /  P 

p        r  i 

where  P  ,  GG,  denote  the  time-average  power  input  and 
i 

S(r,6,<fr)  is  the  time-average  power  density  in  the  radiated 

field.   GPP  and  GTT  denote  the  power  gains  associated  with 

the  phi-polarized  and  the  theta-polarized  components  of  the 

field,  respectively. 

The  use  of  the  variables  JFLAG  and  KFLAG  are  described 
in  subroutine  SGANT. 

CALLED  BY:      MAIN 


CALLS  TO:       GDISS 
GFF 
RITE 
SQROT 
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<C     GFFLD     IIA, IB. INC.INH.IWR.Il,  12.13,  U2tHDiN,NO.NH,*M.ACS 
.CGO.CG.CJ.CMM.O.ECSP.ECST.EP.ET.fiPP.ETT.EPPS.EPts.ETPS.ET 

>,gt i,PH,s&o,scsp,scsr,sppH,spiM,srph,siiH,tH, x.r.i ,  zld.zs 


13(1) ,    NO(l),     MDl INM.4) 


subrout ine 
ip.acst.c.cc 

2T  S,GG  ,i.PP,  Oilirnm 
3,FTA,GAM,ERR, IGRD) 

COMPLEX  ERR 

COMPLEX  CJI.ET1.ET2,EP1,EP2, EPPS, ETTS, EPTS.ETPS.ZS.VP.VT 

COMPLEX    CII),CJ<ll,fcP(ll,ETill.EPP(l).ETTUI,ZL0Il) 

COMPLEX     ET»,GAM,CGO<  ll.SGO     1  I  ,CGI  1 

OIMENSION     [All),      IBID.     Till),      12(1). 

DIMENSION    0(1).     XIII,     YllI,     ill) 

DATA    PI ,IP/3. 14159,6.28318/ 

CJ1 

&GG 

THR 

CTH 

SIH 

PHR 

CPU 

SPH 


-4.»PmETA»GAM) 
REALM  ./ETA) 
.01 74533*TH 
COSI  THRI 
SINI THR I 
.01  7<.533»PH 
COSt  PHR) 
SIN(PHR) 


DO     I     1=1, N 

EIT( | |     •     ( .0, .0) 

EPPI II     -    ( .0..0I 


00    3    K=l ,NH 
KA    =     IAIK) 
KB    =     IB(K) 
NGRO    *     IGRD 
IF    IK.LT.NM/21     IGRD«-1 

CALL    GFF     {X  (k.A)  ,Y  (KAI  .  Z1KAI  ,X(  KB  I  ,Y(KBI  ,  ZIKB)  ,D( 
1TH. STH1CPH,SPH,6AM,ETA,ETl,ET2,EPl,EP2t  iGRO.ERRl 
IGRD    •    NGRD 
NDK    «    NOIKI 


DIKI ,CGD(KI,SGD(KI,C 


DO 


=  1,ND 
C.  II) 


NDK 


3    I  l< 
1     *    MD( K, 
Fl     -     1. 
IF    1KB. EQ.I2I III    GO    TO 
IF     (KB. EO. Ill  II  I     Fl— I. 
EPPI  I  I     •    EPPI  I  1  *F  1*EP1 
ETTII I     ■    ETT( I)»F[*ETl 
GO   TO    3 

IF    (KA.E0.I3I  II  I     Fl  — I. 
EPPI I)    ■    EPPI ll»FI*EP2 
ETTII I     ■    ETTI l)»FI*ET2 
CONTINUE 


EPPS    ■     I  .0..OI 

E  TTS    ■     (  .0..0I 

IF     ( INC.EQ.OI    GO    TO    8 

IF     I  INC .EQ.2I     GO    TO    6 


DO    4    1=1. N 

ETUI     «     ETTI  IMCJI 

EP(  I  I     -     EPPI I l»CJI 

CALL     SCROT     (C,EP,0,I12,NI 
112"=    2 

CALL    SOROT     (C,ET, 0,112, Nl 

IF     (IMR.GT.OI     MR  i  TE     (6,10)     PH.TH 

IF     (IMR.GT.OI    MRITE     (6,111 

CALL    RITE     I IA, IB, INM, IHR.I 1.12, I3.MD.N0, NM.EP.CG, IGRD) 

CALL    GOISS    (AM,c6,CMM,O.POI S  ,  & AH , NH , SGO , Z L6 , Z S I 

IF    I IMR.GT.OI    MRITE    4  6,121 

CALL     RITE     I  I  A,  IB,  INM,  1MR,(  I  .  12,  I3,M0.ND,NM,  ET,CG,  IGRD) 

CALL     GDISS     I AM,CG,CHM,D,tDI S ,GA M , NM , SGO , ZL 6 , Z S I 

ACSP    -    POIS/G&& 

ACST    ■     TOIS/GGG 

PIN    •    .0 

TIN    «     .0 

DO    5    l«WN 

VP    =    cj  l»EPP(  I  I 

VI    -    CJI'ETTI I ) 

PIN    ■     P INfREAL IVP»CONJG(EP( I ) I ) 

TIN    -     T  1N»REAL(VT»C0NJG(ET(  I  I  I  I 


ECSP 

■ 

PI N/GGG 

ECST 

■ 

TIN/GGG 

SCSP 

■ 

FCSP-ACSP 

SCST 

■ 

ECST-ACST 

EPTS 

■ 

l.0,.0) 

FTPS 

- 

1  .0,-01 

00    7     |« I.N 

EPPS    *    EPPS»EP(  l)»EPP( I  1 
EPIS    -    EPTS»FP(  I  I  »FTT(  I 
ETTS    •    ETTS»ET(  I  |*ETM  I 
ETPS    •    ETPS»ET(  I  l»EPP( I  I 


SPPM  •  2.  »TP»  (CABS 
SPTM  •  }.*TP*  (tOS 
STPM  •  }.♦ tP»(CABS 
STTM  .  2.»TP«(CABS 
RETURN 


EPPS|»»2 
E  PT  5  •♦? 
ETPS  •♦? 
ETTSI»«2 


8  00    9    I"  1  ,N 

ETTS    -    ETTS»CJ(  I  I'ETTI  I 

9  EPPS    -    EPPS»CJ(  I )«EPPI I  I 


APP 
ATT 
GPP 
GTT 
RETURN 


CABSIEPPSI 
CABS(ETTS) 

4.«PI»APP»APP»GGG/GG 
4.»PI»ATT»ATT»GGG/GG 


10  FORMAT  (IOX, 'BRANCH  CURRENTS  ASSOCIATED  MI T H  PLANE-MAVE  SCATTERING 
I  FOR  THE  INCIDENT  ANGLES,  PHI ■ •  ,F5.  1  . '  *N0  THETA- • , F5. 1 // I 

11  FORMAT  U*X, 'CURRENTS  INDUCED  BY  ?HE  JH1  POLARIZED W*VE'  I 

12  FORMAT  I44X, 'CURRENTS  INDUCED  BY  THE  THETA  POLARI ZED  M AV E ' ) 


3 

* 

5 

6 

7 

8 

9 
10 
11 
12 
13 
I* 

IS 
17 
18 
19 
20 
21 
22 
23 
2* 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 

49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
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91 
94 
95 
96 

97 
98 
99 
100 
101 
102 
103 
104 
105 
106 
107 
108 
109 
110 
I'll 
112 
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GGMM 


PURPOSE:       to  calucate  the  mutual  impedance  between  two 
filamentary  monopoles  with  sinusoidal  current  distribution. 

METHOD:         As  stated  in  subroutine  SGANT,  the  mutual 

impedance  of  coupled  dipoles  may  be  expressed  as  sum  of  four 

monopole-monopole  impedances.  This  subroutine  calculates  the 

mutual  impedance  with  closed-form  expressions  in  terms  of 
exponential  integrals. 

For  skew  monopoles  it  can  be  shown  that  the 
monopole-monopole  mutual  impedance  is  given  by: 


t  -z  z 

i+j  n  mm 

z        =    (-1)  B[e  (F        -e  G        +e        G) 

ij  j1  12  22 

-t  -z         z 

n  mm 

-e     (F   -e     G    +  e    G   )  ] 
j2         11         21 

where  m=  2/i,  n  =  2/j  and 


B  = 


16  n  sinh  d   sinh  d 
1        2 


The  functions  F   are  defined  by: 
ik 

gz   cos  |J 
i 
F    =2  sinh  d   e  E  (R   +  gz   cos  9  -  gt) 

ik  i  i     i 


k 
where  g  =  (-1)  ,  d   and  d   are  the  lengths  of  the  monopoles 

being  considered.  The  functions  G   are  defined  as  follows: 

ik 
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G    =  E(R   +  gz   +  q't  -  jg»)  +  E  (R   +  gz   +  q»t  +  jg») 
lk       2      2  2      2 

-  E  (r   +  qz   +  q't  -  jq")  -  E  (H   +  qz   q  q't  +  jq") 
11  11 

i  k 

where  q  =  (-1)  ,  q'  =.  (-1)  ,  and  q"  =  qb  +  q'c  with 

b  =  c  cos  f  and  c  =  d/sin  f.  The  angle  V  is  the  angle  formed 

by  the  apparent  intersection  of  the  two  monopoles.   This 

will  be  discussed  later  in  detail. 


ik 
In  the  above  equation  for  G   ,  t  denotes  the  position  of 

an  observation  point  somewhere  on  monopole  2.  R   and  R   are 

12 

the  distances  from  the  endpoints  of  monopole  1  to  this 

observation  point.  Finally,  the  E  functions  are  defined  as 

follows: 


E  (a  +  jg")  =  e 


ai+jg» 


a  +jg" 
2    -yv 

_£ dw 


where  a  and  g"  are  real  quantities  with  dimensions  of 

length,  a  is  a  function  of  t,  a   =  a (t  ) ,  a   =  a (t  )  and 

112      2 


Y   -    jw  \/u6.   The  integral  above  is  evaluated  by  subroutine 
EXPJ. 


U- 


I 
S2 


S1! 


-D---- *- 


T1 


PARALLEL 


t 


T2 


S2  , 


S1  1  p 


/ti 


SKEW. 


T2 


To  explain  the  input  data  for  GGMM,  refer  to  the  above 
figure.   If  the  monopoles  are  parallel,  then  the  new 
coordinate  system  is  defined  such  that  the  new  z  axis  is 
parallel  to  the  monopoles.   The  coordinate  origin  may  be 
selected  arbitrarily.   S1  and  S2  denote  the  z  coordinates  of 
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the  endpoints  of  the  test  monopole,  T1  and  T2  are  the 
coordinates  of  the  endpoints  of  the  expansion  monopole,  and 
D  is  the  perpendicular  distance  (displacement)  between  the 
monopoles.   The  mutual  impedance  of  parallel  monopoles  is 
calculates  in  the  last  part  of  GGMM  below  statement  5. 

For  skew  monopoles,  let  the  test  monopole  s  lie  in  the 
xy  plane  and  the  expansion  monopole  t  in  the  plane  z  =  D.  (D 
is  the  perpendicular  distance  between  the  parallel  planes.) 
If  the  monopoles  are  viewed  along  a  line  of  sight  parallel 
with  the  z  axis,  the  extended  axes  of  the  two  monopoles  will 
appear  to  intersect  at  a  point  on  the  xy  plane.  Let  s 
measure  the  distance  along  the  axis  of  the  test  monopole 
with  the  origin  at  the  apparent  intersection.   S1  and  S2 
denote  the  s  coordinates  of  the  endpoints  of  the  test 
monopole.   Similarily,  let  t  measure  the  distance  along  the 
axis  of  the  expansion  monopole  with  the  origin  at  the 
apparent  intersection.   T1  and  T2  denote  the  t  coordinates 

of  the  endpoints  of  the  expansion  monopole.  Let  s  snd  t"  be 

unit  vectors  parallel  with  the  positive  s  and  t  axes, 

respectively.   Then  CPSI  =  s»^  =  cos  f.    The  monopole  lengths 

are  d   and  d  . 
s      t 

The  output  data  from  GGMM  are  the  impedances  P11,  P12, 
P21,  and  P22.  In  defining  these  impedances,  the  reference 
direction  is  from  S1  to  S2  for  the  current  on  monopole  s, 
and  from  T1  to  T2  for  the  current  on  monopole  t.  In  the 
impedance  P   ,  the  first  subscript  is  1  or  2  if  the  test 

dipole  has  terminals  at  S1  or  S2  on  monopole  s.   The  second 

subscript  is  1  or  2  if  the  expansion  dipole  has  terminals  at 

T1  or  T2  on  monopole  t.   The  monopole  lengths  d   and  d   are 

s       t . 

assumed  positive  in  defining  the  input  data  CGDS,  SGD1  and 
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SGD2. 


For  parallel  monopoles,  CPSI  =  1  or  -1.  S1,  S2,  T1,  and 
T2  are  cartesian  coordinates  for  parallel  monopoles  and 
spherical  coordinates  for  skew  monopoles.  For  skew 
monopoles,  the  radial  coordinates  S1,  S2,  11,  and  T2  tend  to 
infinity  as  the  angle  fJ  tends  to  zero  or  n.   Therefore,  if 
the  monopoles  are  within  4.5°  of  being  parallel,  they  are 
approximated  by  parallel  dipoles. 


CALLED  BY: 


GGS 


SGANT 


CALLS  TO: 


EXPJ 


.SUBROUTINE  GGHM  1SI.S2,T1.T2,0,CG0S,SGDI .SG02.CPSI  ,m,G«,PUl„2 

COHPLEX  El  2  J  |  FI2.5  tL*EKL,ECi',Esl'ES2.En,Er2,EXPi,EXPB 

0S0    ■    0»D 


SGDS    -SG01 
IF    JS2.1T. Sll 

SGOT    -    SG02 

if   iT2.ir.ni 


S30S    -    -SGDl 


ES2  -  CEXP(G4M«S2) 
EM  -  CEXPIGAMM1) 
Er2/oCEXP(GAM.T2. 

opsi  .  cpsi 
km  »   ri 

ID2    .    T2 

CPSS    "    OPSI'OPSI 

cn^«  oo/DsoRTi l.oo-cpssi 

BO    -    CD-DPSI 

B    ■    BO 

CEXP(G»«.tHPLX(.o!c     I 


EC 


00    1    K-1.2 


■  °9  '  '•■' 

1    EIK.U    - 


0,-01 


[si  ■   roi»roi 

TS2    •     r02»T02 
OPg    •    00*00 
SI    •    SI 

00   *    1-1,2 

Fl     -     l-l(»»I 

SOI    -    SI 

SIS-    S0l*S0l 

SU    .    2.«SD!«I0t«Of>S| 

kl2    *„2.»SOI»t02»0PSI 

*i  -  oso»T(opo.5is»fsj-sni 


i 

2 
3 
* 
5 
6 
T 
8 
9 
0 
II 
2 
3 
<. 
5 
6 
7 
8 
9 
20 
21 
21 
3 
_* 
25 
26 
?T 
28 
29 
30 

IS 

3* 
35 
36 

n 

39 
♦  0 

XI 

*3 

it 

46 

*r 

*B 
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R2    -    DSQRMDPQ»S1StTS2-ST2I  *9 

EK    -    EB  50 

51 

00    3    K*l,2  52 

FK    *    I-  I (•»*  53 

SK    *    FK»SOI  5* 

EL    •    EC  55 

56 

DO    2    I' I,  2  5  7 

FL    *    I-  ll»«L  58 

EKL    =    EK»EL  59 

XX    *    FK«PD»FL»CO  60 

TL1     ■=     FIM01  61 

IL2    =    FL»T02  62 

RRl     =    «  1  f  SKML  I  63 

RH2    =    F.2.SKML2  6* 

CAll     EXPJ     (GAM'CMPLXI  HR|  ,-XXI  ,GAM»CNF1X(  RR2.-XX1  ,  EXPAI  65 

CALL     EXPJ     (GAM'CIPLXI  RR1  ,XXi  ,GAM»CMPIXIRR2,XX1  .EXPBI  66 

EIK.L)     =    E IK, LI »F I  » < E X  PA* EKL'E X PB/EKL I  67 

2  EL     *     I. /EC  68 

69 

3  E«    '    1./E8  70 

71 

20    «    SDI'OPS I  72 

2C    •    20  73 

EGZI    •    CEXP(GAM»2C1  7* 

RRl    -    R1»20-TD1  75 

RR2    *    R2»ZD-TD2  76 

CALL    EXPJ     (CAMtRRl ,GAM»RR2, EXPB I  77 

RRl    •    Rl-ZOHOt  78 

PR2    *    P.2-2U»I[)2  79 

CALL     EXPJ     (GAM'ORl ,GAM»RR2,EXPA)  80 

F(i.l)     »    2.«SGDS»EXPA/EG2l  81 

F(I,2I     «    2.»SGDS»EXPB»EGZI  82 

*    SI    -    S2  83 

B« 

CST    -    E TA/I 16.»PI»SG0S»SGDT I  85 

,2)/ES2l»FT2M-Fll,2l-E<2,ll»ES2»  86 


: TA/I 16.»PI»SG0S»SGDT I 
P  I  1    •    CSTXIFI  1,1  )»E(  2,2t»ES2-E<  1,2 
1E<  1,1 1/ES2I/EI2) 
P12    •    CSt»< (-FI 1,11 


1)-EI2.2I*ES2»EI1,2)/ES2)*ET1»IF(1,2I»EI2,I)*ES2-  88 

1EI1, 1I/ES2I/ETI I  89 

P21    ■    CST»U-FI2,11-EI2,2I»ES1»E(1,2)/ES1I»ET2MF(2,2I»EI2,I  »»ESl-  90 

lEll.il/ESll/ET2)  91 

P22  ■  CST»((F(2.1  I»E(2,21»ES1-E< 1 ,2  I /ES I  1 *E  1 1 ♦ ( -F < 2 , 2 ) - E I  2 , 1 I »ES 1 ♦  92 

IE  I  I,  ll/ESI  l/FU  I  93 

RETU&N  9* 

5  IF  ICPSI.LT.O.I  GO  TO  6  95 

TA  •  Tl  96 


18  '  r2  11 

GO  TO  I  98 

6  TA  =  -Tl  99 
TB  •  -12  100 
SGDT  •  -SGOT  101 

7  SI    •    SI  •                                     J02 
C  103 

00    9    1*1,2  10J 

TJ    •    TA  105 

C  106 

DO    8    J=l,2  (07 

2 1  J    =     TJ-SI  108 

R    -    SURMOSUWI  J  *  Z  I  Jl  109 

M    *    R»2 I  J  110 

IF     WIJ.LT.O.l  M    •    0SQ/IR-2IJI                                                                                                                            111 

V    =    R-2  I  J  112 

IF     (2IJ.GF.0.I  V    •    0SQ/1R.2IJ)                                                                                                                   Hi 

IF     (J.EO. II     VI  "    V                                                                                                                                                               11* 

IF     (  J.EQ.l  I    Ml  =    W                                                                                                                                                   1}5 

EG2II ,JI     ■    CEXP(GAM»2I  Jl  }16 

8  T J    -     TB  117 


118 
I  19 
120 


CALL    EXPJ     IGAS»V1  ,GAM»V,GPI  I  II 
CALL    EXPJ     <GAM»H1,GAM»W,GHI  I  II 

9    SI    ■    S2  }|1 

CSI    -    -ETA/I 8.»PI»SGD5»SG0TI  123 

PI  I    -    CSt'(&MI2)»EG2(2,2l»GPI2l/EG2(2,2l-CG0S»<GH( 1 1 »EG2 I  1 ,2  I »GP I  1  12* 

'p12&*    Csf»l-GM(2|»£GZ  I2,1I-GPI2I/EC2(2.  1  I  »CGOS»(  GM(  1 1  »EG2  U.  1 1  »GP(  126 

1 1 1/EG2I  1,111)  127 

■       P21    ■    CSf*IGNIll*EG2( 1.2l»GPUI/EG2t 1 , 2  I -CGOS*. GM( 21 »EG2 (2.2>»GP<2  128 

P22    •    cil»i-GMU  l»EG2U,ll-GPIl  I /EG2C 1 , I  I »CGOS»l GUI  2 l*EG2 < 2, 1  I »GP(  130 

12I/EG2(2,1I II  |31 

RETURN  132 

END  U3 
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6GS 


PURPOSE:       to  calculate  the  mutual  impedances  between 
two  filamentary  monopoles  with  sinusoidal  current 
distributions. 

METHOD:        The  monopole-monopole  mutual  impedance  as 
defined  by  SGANT  is  calculated  using  the  equations  defined 
in  subroutine  GNF.  The  endpoints  of  the  axial  test  monopole 
s  are  (XA,YA,ZA)  and  (XB,YB,ZB),  and  the  endpoints  of  the 
expansion  monopole  t  are  (X1,Y1,Z1)  and  (X2,Y2,Z2) .   DS  and 
DT  denote  the  lengths  of  monopoles  s  and  t,  respectively, 
CAS,  CBS  and  CGS  are  the  direction  cosines  of  monopole  s, 
and  CA,  CB  and  CG  are  the  direction  cosines  of  monopole  t. 

The  effects  of  ground  for  vertical  co-linear  monopoles 
are  applied  in  a  slightly  different  manner  than  mentioned 
previously.  As  With  self  impedance  calculations,  the  test 
monopole  and  the  expansion  monopole  are  laterally  displaced 
by  the  wire  radius.  This  lateral  displacement  is  used  to 
determine  the  angle  of  incident.  This  technique  is  applied 
at  statement  8. 

If  INT  =  0,  GGS  calls  GGMM  for  the  closed  form  impedance 
calculations.  Otherwise  GGS  calculates  the  mutual  impedance 
via  Simpson' s-rule  integration  with  the  following  number  of 
sample  points:    IP  =  INT  +  1.  If  the  monopoles  are  parallel 
with  small  displacement,  GGS  calls  GGMM  to  avoid  the 
difficulties  of  numerical  integration. 

Since  the  point  (X,Y,Z)  of  subroutine  GNF  lies  on  the 
expansion  monopole  t,  T  is  the  integration  variable  and  is 
measured  from  (X1,Y1,Z1).   C1  is  the  current  at  T  for  the 
mode  with  terminals  at  (X1,Y1,Z1),  and  C2  is  the  current  at 
T  for  the  mode  with  terminals  at  (X2,Y2,Z2) .  C  denotes  the 
Simpson • s-rule  weighting  coefficient. 
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Below  statement  7,  GGS  performs  some  analytic  geometry 
in  preparation  for  calling  GGMM.  The  remainder  of  this 
section  is  concerned  with  this  preparation. 

Let  s  denote  a  unit  vector  in  the  direction  from 

(XA,YA,ZA)  toward  (XB,YB,ZB) .  Also  let  T  denote  a  unit 

vector  from  (X1,Y1,Z1)  toward  (X2,Y2,Z2).  Then  s.E  = 

cos  9  =  CC  where  0  is  the  angle  formed  by  the  axes  of  the 

two  monopoles.  Let  monopole  s  lie  in  one  plane  P   and 

s 

monopole  t  lie  in  another  parallel  plane  P  .  CAD,  CBD  and 

CGD  are  the  direction  cosines  of  the  unit  vector  3"  =  f  x  s  / 

sin  0  which  is  perpendicular  to  both  planes.  To  obtain  the 

distance  DK  between  the  two  planes,  a  vector  TZ        is 

ii 

constructed  from  (XA,YA,ZA)  to  (X1,Y1,Z1)  and  take  DK  = 


k  .a. 
1 1 


A  line  is  constructed  from  (X1,Y1,Z1)  to  the  test 
monopole,  such  that  the  line  is  perpendicular  to  the  test 
monopole.  SZ  denotes  the  s  coordinate  of  the  intersection  of 
this  line  with  the  test  monopole,  and  the  cartesian 
coordinates  of  this  intersection  are  XZ,  YZ,  and  ZZ.  The 

direction  cosines  of  s  x  3  are  CAP,  CBP,  and  CGP. 


From  the  point  (X1,Y1,Z1)  in  plane  P  ,  a  line  is 
constructed  perpendicular  to  the  point  (XP1 , YP1,ZP 1 )  in  the 
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plane  P  .   This  line  is  parallel  with  3  and  has  length  DK. 
s 

Let  I  represent  a  vector  from  (XZ,YZ,ZZ)  to  (XP1 ,YP 1 , ZP1)  . 


P1  denotes  TJ  (s  x  "3)  .  S1  and  T1  are  defined  in  subroutine 
GGMM. 


CALLED  BY: 


SGANT 


CALLS  TO; 


GGMM 


SUBRD 

is, oi, 

COM  PL 
COMPL 
COMPL 
COMPL 
COMPL 
COMPL 
COMPL 
DMA 
CA  ■= 
CB  ' 
CG  « 
CAS  - 
CBS  - 
CCS  ■ 
CC 
IF 
IF 
SI 
IF 
INS 
IF  I  I 
IP  - 
DECT 
T  •  . 
OS*  • 
Pll  » 
P12  - 
P21  - 
P22    - 

AMS    ■ 
SGN    - 


UUNE 
SGOI  , 
EX    EX 


OGS    ( XA,YA,ZA,XB,Y8,ZB,XI,Y1,Z1,X2,Y2,22,AM,DS,C&0S,SGD 
NT  ,E  rA,GAM,Pli,Pi2,P21,P2  2,ERR,lGRD) 
,EY1,EX2,EY2,EZI,EZ2 

,P12,P21,P22,EJA,EJB,EJ1,EJ2,ETA,GAM,C1,C2,CST 
.CGOS,SGDS,SGDT,ERl,ER2,ETl,ET2 


I  I 


FP/12 

(X2-X 

(  Y2-Y 

H2-1 

(  XB- 

(  YB- 

1  ZB- 

CA'CA 

CG.LE 

BSICC 

IXI-X 

NI.LE 

2*(  I 

ns.lt 

INS»1 

=  or/ 
o 

cc»o 
i.o, 

I  .0, 
(  .0, 
(  .0, 

AM-i 
-1 


,EXX,EYY 
.PX.PY.PZ 

I,RP.2,RR3,RR4,RHl  ,RV  I  ,  RH2 , RV2 , RN3 ,RV3 ,RH4, RV4 
.56637/ 
I  l/OT 
ll/OT 
I  l/OI 
XAJ/DS 

YAI/DS  |. 

ZAI/OS 
S»CB»CB 
..003) . 
I.GT..9 
AI»C»S» 
0  1    GO 

Nr/2| 
.21   Ins 

ins 

FLT 
.01 
.01 
.01 
.01 
M 


S*CG*CGS 

AND. (C&S. IE. .003). AND. I IGRO.GT.OII    GO    TO    1 

971    GO    TO   6 

(V1-VAI*CBS*IZ1-Z»)*CGS 

TO    7 


00   5    IN-1, IP 

ZZ1    -    SZ 

ZZ2    ■    SZ-OS 

XXZ    -    XI »r»CA-XA-SZ»CAS 

YYZ    ■    Y1»T«CB-YA-SZ«CBS 

III    -    Zl • T»CG-ZA-SZ»CGS 

RS    •     XXZ»»2»YYZ»«2»ZZZ»»2 

Rl    •    SQRT(RS»ZZ1»»2) 

EJA   -    CEXPI -GAM»R| ) 

EJ1     ■     EJA/Rl 

R2    •    SURI(RS*ZZ2**2I 

EJB    -   CEXP<-GAw»R2> 

EJ2    ■    EJB/R2 


1 

2 

3 

* 

5 

6 

7 

8 

9 

0 

1 

2 

13 

4 

5 

6 

7 

8 

9 

0 

21 

22 

23 

2* 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 
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EM    ■    E  JA»SOOSWZ  l»EJl»CGDS-ZZ2»EJ2  *9 

ER2    -    -fcje»SG0S»Z22»EJ2»CG0S-2Zl»EJl  50 

FAC    «     .0  51 

IF     IRS.GF.AMS)     FAC     ■     I C A*XX Z »C B»YYZ »CG» I  II I /OS  52 

EM    =    CCMt  J2-CJ1*CGDSI»FAC*ER1  53 

ET2    -    CCME  J1-EJ2»CG0S  UFAC*ER2  5* 

IF    < IGRD.LT.OI    GO    TO    *  55 

RV1     •     1  -1.  ,01  56 

Rill     -     1  -1.  ,0)  5f 

RV2          I  -  1 . . 0 1  58 

RH2    =    (  -1.  ,01  59 

IF    (IGRU.EO.ll    GO    TO    2  60 

XG1    •    Xl»T»CA-XA  61 

YG1    «     >  l»MCB-YA  62 

2G1    ■     I  1 ♦ I*CG-2 A  63 

XG2    =    X  1«T»CA-X9  6* 

Y&2    «    Y1»T»CB-YB  65 

102    «    ZIMKG-ZB  66 

RG1    «    SQRMXG1*XG1»YG1*YG1I  67 

R&2    -    SURT IXG2»XG2»YG2»YG2I  68 

IT1    «    A  TANI RGl/ZGl I  69 

TT2    -    A IANIRG2/ZG2I  70 

CtHl    '    COS(  Mil  71 

SSTH1     *     SIN! til I*SINITT1)  72 

C1H2    -    C0SIIT2I  73 

SSIH2    =    S1N(  T  r  2  1  •  S 1  N(  TT2I  7* 

RR|    •    CSQRTI  ERR-SSIH1  I  75 

RH1     ■     ( CtHl -RR1 I / (CTHlfRRl I  76 

RV1     •     -!FRK«cTH|-RRll/(ERR»CTHl»R(lll  77 
RR2    ■    CSuMIERK-SStHjl 

RH2    -    ICIH2-RR2|/(CIH2»RR2I  79 

RV2    ■    -»ERR«CTH2-RR2I/1ERR»CTH2»RR2 I  80 

RG    -    S0RT1  IXB-XA)*(X8-XA)M  YB-Y A ) •(YB-YA I  I  81 

CPH   -    0  82 

SPH   -    0  83 

IF    (RG.LT.1.E-32I    GO    TO   3  8* 

CPH   •    ( XB-XAI/RG  85 

SPH    «     IY-B-YAI/RG  86 

EXX    »     ETI'CAS  87 

EYY    ■    E  T1»CBS  88 

EE    ■    (EXX«SPH-EYY»CPH)»(RH1-RY1)  89 

EX1     •     EXX»RV1 tEE'SPH  90 

EY1     ■     EYY*RV1-EE*CPM  91 

Ell    «    -EI1*RV1*CGS  92 

ET1    ■    EXI»CAS»EY1*CBS»EZ1»CGS  93 

EXX    -    ET2»CAS  9* 

EYY    -    ET2»C8S  95 

EE    ■    (EXX»SPH-EYY»CPHl»lRH2-RV2l  96 


EX2    •    EXX«PV2»EE»SPH  07 

EY2    «    EYY"RV2-f E«CPM  98 

EZ2    ■     -FT2»CGS»RV2  99 

EI2    ■    EX2»CAS»EY2»CBS»EZ2»CGS  100 

*    C    -    3.»SGN  101 

IF    ( IN.EU.I .OR.IN.EQ.IPI    C-l.  102 

E&D    *    CEXP(0AM»(0T-TI  I  103 

CI    =    C»IFGU- I./EGD1/2.  .                                  10* 

E&D    -    (  EXPI(.AM»T  I  105 

C2    •    CMEGD-1./EGDI/2.  106 

Pll    •    PI  I  »  b  T 1»C|  107 

P12    -    PIM  I  I  »C  2  108 

P2I    ■=    P2l»ET2»Cl  109 

P22"P2?»ET2»C2  110 

T    ■    T»0EIT  III 

SZ    *     SZ»OSZ  112 

5  SON    •    -S&N  l|  J 

C  115 

CST    •     -EIA»JELT/ 13.»FP»SG0S»SG0TI  116 

Pll    '    CSI«PI1  i|7 

P12    =    CSI»P12  1 18 

P21    -    CSI«P?t  I  19 

P22    ■    CM»P22  120 

RETURN  121 

6  S2 1    ■    (xl-XAI»CASMYl-YA)»CBSMZl-ZAi»CGS  122 
DR1    >    SQRT<(Xl-XA-SZl»CAS)«»2»tYl-YA-SZl»CBSl»*2MZl-Z»-SZl*CGS>«»  123 

12  I  124 

SZ2    ■    SZl»DT»CC  125 

0R2    '    SQRT(IX2-XA-SZ2»CASI»»2MY2-YA-SZ2»CBS»»»2MZ2-ZA-SZ2»CGS1»»  126 

•  2  I  127 

000    «     (DRI»OR2l/2.    ■  128 

IF     lOOU.&r ,20.*AM.AN0. INT.GT.OI  GO    TO    1                                         .                                                        129 

IF     lulm.l  T  .AMI     DOO    •     AM  130 

CALL    GGMM     I  .0,DS,SZ1,SZ2,DDD,CGDS,SGDS,SGDT  ,  I  .  .  E  T  A  ,G  AM  ,  P  I  1  ,  P  I  2  ,  P  2  1  131 

1.P22I  132 
If    I  I  i. m i.i  i  .,i|    RETURN 
IF    (  IGPO.GT.ll    GO    TO    8 


33 

Pll    -    -PI j     """     ""  35 

P12    ■    -PI2  36 

P21    ■    -P21  37 

P22    ■    -P22  38 

RETURN  )9 

SS    •    SORT  U.-CC»CCI  ll 

CAD    •     ICGS«CB-CBS»CG1/SS  *2 

CBO   -     (CAS»CG-C5S»CAI/SS  U3 

C&O   ■    ICBS»CA-CA5»C8I/SS  1*4 


5B 


OK    =    tXl-XAI*CAD»<Yl-YA)*CBOHZl-ZA>*C30  1*5 

OK    •     ABSIOKI  1*6 

IF     (OK.LT.AH)     OK    •    AM  1*7 

XZ    »    XA»SZ»CAS  Us 

YZ    ■    YA»S2«CBS  1*9 

11    •     ZA»SZ»CGS  150 

XPI     -    Xl-OK»CAD  151 

YP1    -    Y1-DK«CBD  152 

ZPl    •    Zl-OK»CGO  153 

CAP    .    CBS«CGO-CGS»CBD  15* 

CBP    •    CGS»CAO-CAS»CGO  155 

CGP    -    CAS»C0O-CBS*CAO  156 

PI    «    CAP»(XP1-XZI»CBP»CYP1-YZ)»CGP»IZP1-ZZI  15/ 

II    ■    PI/5S  158 

SI    •    Tl»CC-SZ  159 

CALL    GGMN    ( S 1 , S 1  *D S , I  I , T 1 »D T ,OK , CGOS , SGOS, SGDT , CC , E T A ,GAM , PI  1 , PI  2 ,  160 

1P21.P22)  161 

RETURN  162 

169 

)    AHS    '    AM«AM  16* 

RG    =     (X1-XAI»IXI-XA)»(Y1-YAI»(YI-YA»  165 

IF     IRG.LT.AMS)     PG    •    AMS  166 

OG   *    S(IRTI(21-ZAI*tZI-ZAI»RG>  lo7 

CPH    -    ABSIZI-ZAI/OG  168 

SSPH    ■    RG/Ol'.»OG  169 

PR1    ■    CSCRTCERR-SSPHI  170 

RV1    •    -(EHR»CPH-RRII/(ERR»CPH»RR1)  171 

Pll    '    P11»RV1  172 

PG    =    (Xi-XBI*(X1-X8)HY1-YR1»IYI-YB>  173 

IF     IRG.LT.AMSI     RG    =    AMS  17* 

OG   *    S0RTIIZI-ZBI*IZ1-2BI»RGI  175 

CPH   =    ABSIZ1-ZB1/DG  176 

SSPH    «    RG/OG«OG  177 

PR1    =    CS0RTIERR-SSPH1  178 

RV1    -    -iERk»CPH-RRl)/(ERR»CPH*RRl)  179 

P12    ■    PL2*RV1  180 

RG    -    (X2-XA>»(X2-XA»»(Y2-VAI*<Y2-VA1  181 

IF     (HG.LT.AMS)     RG    -    AMS  182 

OG   -    SORT((Z2-ZAI»IZ2-ZA)»RG»  163 

CPH   «    ABS<Z2-ZA)/0G  IB* 

SSPH    •    RG/OG»OG  185 

RR1    •    CSCjnri  ERR-SSPH)  186 

RVl    ■    -1ERR*CPH-RR1I/IERR*CPH»RR1 I  187 

P21    •    P21*RV1  188 

RG    -    (X2-XBI»(X2-XBl»(Y2-YBt*<Y2-YB>  189 

IF    (RG.LT.AMS)    RG    •    AMS  190 

OG   •    SQRTH  Z2-ZBI»tZ2-ZB)»RG>  191 

CPH   •    ABSIZ2-ZBI/0G  192 


SSPH   •    RG/DG*DG  193 

RR1    •    CSQRT(ERR-SSPH)  19* 

RVl    •    -CERR*CPM-RR1 )/<ERR»CPH»RRll  195 

P22   •    P22*RV1  196 

RE7URM  197 

END  198 
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GNF 


PURPOSE:       to  calculate  the  near-zone  electric  field  of 
a  sinusoidal  electric  monopole. 


z 
i 
i 


E   4 
z 


^  E 


z   ' 

2    *•" 


I 

d 


z  w 
1 


2  > 


:>•' 


•• .-- 


METHOD: 


An  electric  line  source  is  located  on  the  z 


axis  with  endpoints  at  z   and  z   as  shown  in  the  above 

12 

figure.  Let  the  electric  monopole  have  the  following  current 
distribution: 

I   sinh  /(d  -  1)  +  I   sinh  yi 

1(1)  =  -u_! . ? 

sinh  yd 


where  I   and  I   are  the  endpoint  currents-  y  is  the  complex 

12 

propagation  constant  of  the  medium,  d  =  z   -  z   is  the 

source  length.  The  cylindrical  components  of  the  field  are 
E(6)  =0  and 

E(C)  =  ^= [  (I   e    i  -  I   e    2)  sinh  yd 

4ne  sinh  yd     »  2 

-ya 

+  (I   cosh  yd  -  I  )  e   l  cos  9 
12  1 

-yR 

♦  (I   cosh  yd  -  I  )  e    2  cos  e  ] 
2  1  2 J 
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rj  -/R 

E(z)  =  [  (I   -  I  cosh  yd)   e 

an  sinh  yd     »     2  R 

2 

-/R 
+  (I  -  I   cosh  yd)  e   1  ] 
2     1  R    J 

1 

where   is  the  intrinsic  impedance  of  the  medium  and  where 
(6,<J>,z)  denote  the  cylindrical  coordinates  in  a  coordinate 

system  centered  at  the  endpoint  of  z  . 


These  expressions  exclude  the  field  contributions  from 
the  point  changes  at  the  endpoints  of  the  line  source,  since 
these  charges  disappear  when  two  monopoles  are  connected  to 
form  a  dipole. 

Let  the  coordinate  s  measure  distance  along  the  test 
monopole  with  the  origin  at  (XA,YA,ZA) .  From  any  point 
X,YrZ,  a  line  is  constructed  perpendicular  to  the  monopole. 
SZ  denotes  the  s  coordinate  of  the  intersection  of  this  line 
with  the  monopole.  The  length  of  the  line  is  the  radial 

2 
coordinate  6,  and  RS  denote  6  .  R1  and  R2  are  the  distances 

from  (XA,YA,ZA)  and  (XB,YB,ZB)  to  the  point  (X,Y,Z)  . 

In  the  statements  above  statement  1,  the  above  equations 
are  solved;  and  after  statement  1,  the  cartesian  components 

(E  #E  ,E  )  of  the  field  are  determined.  If  a  ground  plane  is 
x   y   z 

present  (IGRD>0)  the  reflection  coefficients  are  applied  to 

the  cartesian  components  before  returning  to  the  calling 

program. 

CALLED  BY:      GNFLD 


CALLS  TO:       NONE 
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SUBROUMNE    GNF     (  X A , Y A , Z A. XB , YB , ZB ,X , Y , I , AM, OS , CGOS . SGOS i ET»t C»H, EX  1 

1  l,f Yl,t Zl ,txz,E Y2.EZ2, IGRO.ERRl                        •  2 

CCMPLEX     ERR.KV1 ,RH1 ,RV2,RH2,RR1 ,RR2,EE  3 

CtlMPifcX    E JA.E JB,EJl,EJ2,ERl,ER2,ESl,ES2,SGDS,GAM,CST,CGOS,ETA  * 

COMPLEX    EX1.EYI , E Z 1 ,E X2 , EY2 ,E Z 2  5 

DATA    PI/3. 1*159/  6 

CAS    =     I  XB-XAI/DS  7 

CBS    *     ( YB-YAI/OS  8 

COS    =     ( ZB-ZAI/OS                                                                                                                                              .  9 

SZ     =     I X-XAI«CAS»IY-YA)»CBS» IZ-ZAI»CGS  10 

ZZ1    ■=    SZ  11 

III    *    SI-OS  12 

XXZ     ■     X-XA-SZ»CAS  13 

YYZ     =     Y-YA-SZ»CBS  1* 

III    =    Z-ZA-SZ»CGS  15 

RS    =    XXZ»»2»»YZ»«2»ZZZ»*2  16 

Rl     *    SURT(RS»ZZl*»2l  17 

E  JA    *     CEXP|-GAM»R1  I  18 

EJl    -    tJA/Pl  19 

R2    =    SO«MRS»ZZ2»»2l  20 

EJB    •=    CEXP(-GAM»R2I  21 

EJZ    =    EJB/R2  22 

ESI    ■    EJ2-tJl«CGDS  23 

ES2    ■    EJ1-EJ2»CG0S  2* 

E«l    '    i  .0. .0)  25 

ER2    =    (  .0,  .01  26 

AMS    =     AM»AM  27 

IF    IRS. LT. AMS)    GO    10    1  28 

CIH1    «    ZZ1/R1  29 

CIH2    =■     lll/«l  30 

ER1     >    I E JA«SG0S»EJA»CGDS»CTH1-E JB«CTH2I/RS  31 

EK2    -    l-EJB»SGOS»E JB »CGDS »C TH2-E JA»C TH1 » /RS  32 

1    CSI    -    f  TA/C.»P|»SGOSI  33 

EX1    =    CST»IES1*CAS«ER1*XXZI  3* 

EY1    ■    CSI»IES1»CBS»ER1»YYZ I  35 

EZ1    ■    CSr»(ESl»CGS»ERl»ZZZI  36 

EX2    -    CST»(ES2»CAS»ER2»XXZ|  37 

EY2    -    CST»(ES2»C8S»ER2»YYZ  38 

EZ2    -    CST»(ES2«CGS»ER2»ZZZi  39 

IF     tlGRO.LE.Ol    RETURN  *0 

RV1    -     1-1. ,01  »l 

RH1    -    1-1. ,01  *2 

RV2    •    1-1. ,01  *3 

RH2    ■    I -1..0I  *♦ 

IF    I IGRO.EO. 1  I    GO    TO    2  45 

Kl     ■     SORT! IXA-X)*(XA-X)»(YA-Y)»1YA-YI I  46 

R2    -    SORT! (XB-X1 *IXB-X) ♦! Y8-YI»IYB-Y) I  *7 

TH1    •     ATANIR1/I  ZA-ZII  *8 


TH2    -     A  TAN(B?/(  ZB-Z)  1 

RRl     ■    C SURT ( ERR-S INI TH1  1 • SI N( TM1  I  I  *9 

55f  ;  £?KI!fMTiw!*Mi?il81iSJIi  I? 

52 

SPM    -    0  57 

i-'dJRG:LI*'-e-3*'    GO    TO    3  58 

CPH    •     ( XB-XAI/RG  59 

SPH    -     IYB-YA1/RG  60 

lx,Mxxi:Rvpr;iE'ispSH,*,RM,-Rvi'  i\ 

III :  lira;!!**"'  U 

ix2^,fxx^2:Rsvp^I^;psH'MRH2-RV^,  tl 

cTI  "  |?2:fV2-EE»CPH  67 

EZ2  •  EZ2*I-RV2I  68 

RETURN  f| 

ENO  ?5 

Tl 
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GNFLD 

PURPOSE:       to  calculate  the  near-zone  electric  field 
intensity  at  a  given  point. 

METHOD:        This  subroutine  calls  GNF  for  the  near-zone 
field  of  each  wire  segment,  and  sums  over  all  segments  to 
obtain  the  near-zone  field  of  the  wire  antenna.  FI  is  used 
in  a  manner  similiar  to  FI  of  subroutine  SGANT.  CJ(I)  is  the 
loop  currents  calculated  by  subroutine  GANT1. 

The  use  of  the  variables  JFLAG  and  KFLAG  are  described 
in  subroutine  SGANT. 


CALLED  BY: 


MAIN 


CALLS  TO: 


GFF 


SUBROUTINE    GNFLO    (U.IB,lNM.U.I2,.3.HO.N.NO,NM,*H.CGO.S6D.ET».G»H 

K!io^ii!i"DiB!iN0iiii».  «2">.  •*"..  °'i»- xm' vu'-  m 

ex  •  ( .6. .oi 

EY    *    ( .0, .0) 

El    "     ( . 0, .01 


DO 

»» 
KB 
NOhO 


KM.NH 
UIKI 
IBIKI 
IGKO 


If    CK.tT.NH/Z1     IGRD.-l 


l(K|,S&OIKI,ETA,G»M,EXl,EYl.Ell.EXZ,EYZ,£22.IGR0.Ef<RI 


IlKBI,XP,YP,II>.»H,D(KI,CG0 


IGRO 
NDK    • 


NGRD 
NOIKI 


DO    2     Il-liNOK 
I     •    HOIK,  I  II 


FI 
IF 
IF 

EX 
EY 

11 

GO 
IF 
EX 

EY 

tl 


(KB*.  EQ.I2I  III    GO    TO 
1KB.  EQ.  11(11]  .f|"-l« 
«    EX»FI*EXl*CJIII 

•  EY»F1*EY1*CJI  I  I 

•  E2»FI*EZ1*CJ(I) 

(2»?EQ.I3<III    FI--I. 
-    EX»FI»EX2»CJII I 
.    EY»FI»EY2»CM I | 

•  EZ«FI*EZ2*CJl  li 


2   CONTINUE 


RETURN 
END 


t 
2 
3 

* 
5 
6 
7 
8 
9 
10 

,2 
13 
1* 
15 

ll 

IB 
19 
20 
21 
22 

li 

25 
26 
27 
2B 
29 
30 
31 

?! 

3* 
35 
36 
37 

38 
39 

*0 
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LEFT 


PURPOSE:        to  determine  position  (location)  of  the  left 
paren  symbol  on  the  input  data  card. 

METHOD:         The  character  search  begins  in  the  column 
passed  to  the  subroutine.  On  returning  to  the  calling 
program  the  argument  passed  is  the  column  following  the  left 
paren  symbol. 

CALLED  BY:      READ 
CALLS  TO:       NONE 


SUBROUI INE    lEf T     (Nl  1 

COMMON    /»/    4(801  2 

DAT*    PLEFT/M  •/  -3 

K     -     U  * 

5 

PO    1     I-K.80  6 

n  •   i»i  7 

IF    Oil  I.EO.PLEFT)    GO  TO    2                                                                                                                         8 

1  CONTINUE  9 

10 

N    -    1  .11 

2  PETURN  12 
END  13 
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LINECK 


PUflPOSE: 


to  insert  grid  characters  on  the  polar  plot. 


METHOD:        The  peroid  character  (ISM  (2))  is  inserted  in 
the  proper  position  in  the  statements  above  statement  U.  In 
the  statements  after  statement  4,  the  grid  numbers  labels 
are  inserted  on  the  horizontal  axis. 


CALLED  BY: 


POLPLOT 


CALLS  TO: 


NONE 


SUBROUTINE    LINECK     IX, Y) 

THIS    SUBROUTINE     INSURES   ALL    GRID   CHARACTERS    LIE    ON    THE    POLAR    GRID 

CCHHLiN     ISYH.LINE 

INTFOER     Y 

DIHtNSIUN     ISYHW4I,     L1NEI130I 

IF     IY.tQ.0l     GO    TO    3 

K    =    0 

IF     IX. L  T.  10.01     GO    TO    5 

SET   UP   AREAS    OF    •PERIOD"    POLAR    GRID    POINT    CHARACTERS 

I     •     INTIXI 

I     •     IABSI  I  1 

I    «    ABSIX1 

IF    llt-ll.GT.0.5)     l«I»l 

IF    W.LT.lO.O.DR.J.GT.lll.OI    GO    TO   2 

IINIIII      •      ISYMI2I 

LlNKfeO)     -     ISYHI3I 

LINEI62I     -     ISYMlli 

K    •    K»l 

IF     IK.E0.2I     GO    TO    2 

I     •     12?- I 

GO   TO    I 

LINE  161  I     ■     ISYHI2I 

IF  IV.NE.OI  GO  10  5 


3  DO  4  K=  11 

111 

LINEIKI  • 

ISYNI2 1 

c 

*  CONTINUE 

I 

C 

F  ILL  IN  GRIO  NJMBER 

LINE! Ill  ■ 

■  ISYMI7) 

LI  NCI  20  1  • 
L INEI 21  1  • 

•  1 SVHI 1 0 1 

■  ISYHI5I 

LINEI22I 

•  ISYHI  1  1  1 

LINEI30I  • 

•  ISYMI9I 

LINEI31 1  • 

•  ISYMI5I 

L INEI 32 1  • 

■  ISYHI III 

L INt 1 40  1  • 

•  ISYHI8I 

L  INK  M  I  • 

•   ISYHI5I 

LINEI42I 

•  ISYMI1I) 

LINEI50)  • 

•  1SYNI7I 

L INEI 51  1 

•  1SYH15) 

LINEI52I 

>  ISYHIUI 

L  1  NE  1  6  1  I 

•  ISYHI  1  1 
■  ISYHI 7) 

L INEI TO  1 

L INEI 71  1 

■  ISYMI5I 

LI  NCI  72  1 

•  ISYHI  |  1) 

L  INF  1  80  1 

■  1SYHI8I 

l inei ai  i 
LiNEiaii 

•  ISYHI5I 

•  ISYHI  I  I  ) 

L  1  NE  (  90  1 

■  JSYHI91 

LINE  191 1 

•  ISYHI5I 

LINEI921 

•  ISYHI  11  ) 

LINEIIOOI 

•  ISYHIIOI 

LINE  110 11 

•  ISYHI5I 

•  ISYH  111 

LINEI102I 

LINE! 11 II 

■  ISYHITI 

5  CONTINUE 

RETURN 

END 

I 

2 
3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

1* 

\l 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
3* 
35 
36 
37 
38 
39 
40 
41 
♦  2 
43 
44 
45 
46 
47 
48 


49 
50 
51 
52 

n 

55 

56 
57 
58 
59 
60 
61 
62 
6) 

» 
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NUMB 


PURPOSE: 


to  place  degree  numbers  on  the  polar  plot. 


METHOD:  The  current  line  which  is  being  printed  is 
passed  to  the  subroutine  in  the  calling  argument.  If  this 
line  contains  degree  numbers,  these  numbers  are  placed  in 
the  correct  position  by  the  IF  statements. 


CALLED  BY: 


PTPLOT 


CALLS  TO: 


NONE 


SUBRUUI INE    NUMB     ( Yl 

THIS    SUBROUTINE    PUTS    DEGREE    NUMBERS    ON    POLAR   GRIO 

COMMON     ISYM.UNE 

INTEGER     V 

OIMFNSION     ISYMI14),     L  t  NE  I  I  30  ( 


IF  IV. Nf 
I  I NE  I  3  > 1 
LINE  (  3<.) 
MNEI  351 

H'lllB/l 

L INEI 6B I 
L INE I  891 
IF  I  Y.NE 
l INEI 1?|  • 

iiMini    ■ 

UNI  I  1<.  |  . 
I INEI 108)  ■ 
I.  INE  I  1091  = 
LINE! 11 0)  = 
IF  (Y.NE. 01 
L INEI 71  = 
I INt I Bl  ■ 
L INE (91  = 
L  INE  I  II  31 
I  INt (  II  4  1 
L INt  (  1  I  5) 
IF  IY.NE. 
L INF ( li  I 
MNt(  I  II 
L INEI  14  I 
LINCI 1081 
I  INE  I  109 


3  71     GO     TO 

■  I  S  Y  M  I  7  I 
=  ISYMI8I 
=     ISYMI6I 

-  ISYMI6I 

*  I SYMI  1 2  I 
«     ISYMI6) 
211    GO    TO    2 

•  I SYMI 71 

■  I  SYMI  t  I  | 

-  I SYMI 6  I 
I SYMI 6  1 
ISYMI 9) 
ISYMI6I 

.  GO  TO  3 
I SYMI 71 
ISYMI 1 3) 
ISYMI6) 

■  I SYMI 61 
»  ISYMI6) 
=  ISYMI6I 
211  GO  TO 

I SYMI8I 
ISYMI 71 
ISYMI6I 

■  ISYMI9I 
.  -  ISYMI91 

LINEI 1101     •     ISYMI6) 
IF     IY.NE. -37)    GO    TO 


I INEI 33 
LlNEI  3<> 
I INFI 35 

I    INI   (  o7 
L INEI 88 
L INEI 89. 
CONT  INUE 
RETURN 
END 


ISYMI 
I  SYMI 
I  SYMI 
ISYMt   . 
ISYMI6I 

isymUI 


3 
4 
5 
6 
7 
8 
9 
10 
II 
12 
13 
I* 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
3* 
35 
36 

11 
39 
40 
*1 
♦  2 
43 
44 
45 
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NUMBER 

PURPOSE:       to  convert  alpha-numeric  numbers  to  floating 
or  fixed  point  numbers. 

METHOD:        After  initially  determining  the  sign  of  the 
number,  the  DO  LOOP   ending  at  statement  fa  scans  each 
character  beginning  at  HI.  The   DO  LOOP  ending  at  statement 
3  terminates  the  outer  DO  LOOP  if  the  character  being 
compared  is  not  an  alpha-numeric  number.  The  DO  LOOP  ending 
at  statement  5  converts  the  alpha-numeric  number  to  an 
actual  number.   Below  statement  7,  the  multiplier  correction 
is  applied  to  the  floating  point  number  before  returning  to 
the  calling  program. 

CALLED  BY:      READ 
CALLS  TO:       NONE 


SUBROUTINE  NUMBER  (N1,N2,X,IXI  1 

COMH(,N  /A/  4(801  } 

DIMfNS [UN  B( 101  3 

DMA  B/'O','  l'1,2*,a3<t>A>*t>5al,A*ttT>f'a»*t*9»/  * 

DATA  AMNUS. PLUS, POINTS1-' ,'♦•,'.'/  5 

OATA  AK.  ,J4,«U/'K  '  ,'M'  ."U1/  6 

N  ■  Nl  J 

NSIGN  •  0  a 

I  I  ■  -1  9 

IX  «  0  10 

I  SET  •  0  11 

IF  IAIN). EQ. PLUS)  N*N»1  17 

IF  (AINI.NE.AMNUS)  GO  TO  1                            .  13 

NSICN  ■  1  14 

N  ■  N»l  15 

1  00  6  l*N,B0  17 
IF  (All  I. NE. POINT)  GO  TO  2  18 
I  SET  -  1  19 
GO  TO  6  70 

2  IF  (ISET.E0.lt  11-11*1  21 

DO  3  K=1,10  23 

IF  (All  I.EQ.BIKII  GO  TO  <>  2* 

3  CONT INUE  25 

GO  TO  7  |* 

4DO!K'I,IO  ?9 

KK  "  K- I  30 

IF  IAI I I.EO.B(K) )  NUrtB-KK  31 

5  CONTINUE  32 

IX  -  NUMB»10»IX  3* 
N2  -  l»l  35 

6  CONTINUE  36 

7  IF  (NSIGN.EQ.l)  |X  •  -IX  38 
Y  -  IX  39 
IF  II 1.LT.0I  II  •  0  40 

X  -    V/(  10** I  I  I  41 

IF    (AIN2I .EO. POINT)    N2-N2»l  42 

IF     jAIN2l.E0.AK)     X    •    X*1000.  43 

IF     (A(M2|  .Fli.AH        X    •     X'0.001  44 

IF    <AtN2).E0.AU      X    ■    X»0. 000001  45 

|FI|A(N2I.EQ.AX).OR.IA(N2I.EQ.AH).OR.IAIN2).EO.AUI)    N2-N2H  46 

Nl    •    Hi  47 

RETURN  48 

ENO  49 
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POLPRT 

PURPOSE:       to  control  the  plotting  of  the  polar  plot. 

METHOD:        This  subroutine  is  the  main  subroutine  in  the 
polar  plot  package  and  is  responsible  for  calling  the 
various  subroutines  of  the  package. 

The  scale  factor,  S,  must  be  changed  according  to  the 
printer  characteristics.  The  scale  factor  in  this  subroutine 
is  set  for  ten,  10,  characters  per  inch  for  the  abscissa  and 
eight,  8,  characters  per  inch  for  the  ordinate  axis. 
Therefore  S  =  10. /8. 

After  initializing  DATAX,  DATAY,  and  X,  the  input  data, 
Y,  is  scanned  to  determine  the  normalizing  factor.   If  this 
normalizing  factor  is  less  than  1.E-32,  an  error  statement 
is  printed  and  the  plotting  is  abortted. 

In  the   DO  LOOP   ending  with  statement  8,  each  line  of 
the  polar  plot  is  printed  after  a  call  is  made  to  PTPLOT  to 
establish  the  ploar  grid  information.   The  variable,  DIM,  is 
used  to  as  a  scaling  factor  for  the  polar  plot.  The  value  of 
1.0  will  cause  all  of  input  data  to  be  plotted,  however,  if 
only  the  values  less  than  one-half  of  the  normalizing  factor 
are  of  interest,  then  DIM  can  be  set  to  .5.  This  will 
enlarge  of  the  center  of  the  polar  plot. 

CALLED  BY:      MAIN 

CALLS  TO:       PTPLOT 

SART 
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c 


c 


c 


c 
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SUBROUTINE    POLPRT     INANE, Y>  1 

COMMON     1SYM.LINE  2 

DIMENSION    X(360l,     Y1360),    DATAX1360I,    DATAYI360),     LINEI130),     1SYNI  3 

1141  4 

DIMENSION    TITLAI2I.    T1TL2I2I  5 

DATA    TITLA/'PHI     •,'THET'/  6 

N    ■    360  7 

DIM   -1.0  8 

NST    ■     I  9 

KST    •     1  10 

C  It 

C              S    IS    SCALE    FACTOR    OF    PRINTER!  12 

C              ABSCISSA   CHAR.    PER    INCH   /    ORDINATE    CHAR.    PER    INCH  13 

C  I* 

S    ■    10. 0/8.0  15 

C      ZERO  OATAX  AND  DATAY  17 

C  18 

C  19 

DO  I  IAM.N  20 

D  -  IA-1  21 

DATA  XI  IA)  »  0.0  22 

DATA  Y( IAI  «  0.0  23 

1  XIIAi  -  D»3. 1*15927/180.0  2* 
C  25 
C  26 
C  FACTOR  IS  THE  NORMALIZING  DIVISOR  27 
C  28 

FACTOR  ■  Y( 1  I  29 

C  30 

DO  2  !'=-<?,  N  31 

2  IF  IFACTOR.LT.YI IAII  FACTOR-YIIA)  32 
C  33 
C  3* 

IF  INAME. EQ. II  TITLI'TITLAI  II  35 

IF  (NAME.EQ.2)  TITL1«TITLA(2)  36 

IF  I INAME. EO. 31. OR. (NAME.EO. 4). OR. (NAME.EO. 7  I .OR . I  NAME. E0.8 1 1  TITL  37 

12I1MTI  TLAUI  38 

IF  ( (NAME.EO. 51 .OR. INAME. EQ. 61. OR. INAME. EO. 91. OR. INAME. EO. 10) I  TIT  39 

IL2I  1 1  =T  ITLAI2)  40 

IF  (INAME. EO. 31. OR. (NAME.EO. 51. OR. INAME. EQ. 71. OR. INAME. EO. 9  I)  TITL  41 

12(2I>TI  TLAIll  42 

IF  ((NAME. E3. 41. OR. (NAME.EO. 6)  .OR. (NAME.EO. 8). OR. INAME. EQ. 101  I  TIT  43 

1L212I-I ITLA12)  44 

IF  (FACTOR.  G-T.  I.  E-32)  GO  TO  3  45 

IF  (NAME. LE. 21  WRITE  16.91  TITL1  46 

IF  INAME. GE. 31  WRITE  16,101  TITL2  47 

RETURN  48 


49 


C       NORMALIZE  DATA  TO  ONE  50 

C  *• 

C  52 

3  DO  4  IA-I.N  53 

4  Y(  IAI  •  YllAl/FACTOR  54 
C  55 
C  56 

IF  INAME. LE. 21  WRITE  (6,111  TITL1, FACTOR  «..,.-„_               \l 

IF  (  (NAME. GF. 31. AND. (NAME. LE. 61  I  WRITE  16,131  TITL2. FACTOR  58 

IF  (NAME.GE.7)  WRITE  16,12)  TITL2, FACTOR  59 

C       FILL  OATAX  AND  DATAY  ARRAY  FROM  X  AND  Y  ARRAY  60 

C  -                                                       61 

I  « 

DO   5    IAM.N  JJ 

DATA    XI  1  At     ■=    Yl  IA)»COSIXI  IAI)  64 

5  DATA  YIIA1-  Yl  I  A  I • S I N 1 X( I  A  I  I  65 
C  ** 
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C       SORT  DATA  BY  OROINATE  MAGNITUDE  68 


69 


CALL  SART  (OATAX, DATAY, Nl  JO 

C      DATAX  AND  DATAY  ARE  SORTED  BY  OESENDING  MAGNITUDE  ON  THE  OATAY  VAL  72 

C      SET  UP  FOR  PLOTTING  POLAR  GRIO  WITH  DATA  £3 

C  '* 

c  75 

DO  8  IYY-1 ,81  J* 

C  •' 

CALL  PTPLOT  I IYY.SI  £8 

C      LINE  IS  RETURNED  WITH  POLAR  GRID  INFORMATION  80 


81 


C              SET    UP    ">•    BIN    SIZE    UPPER   AND  LOWER  LIMITS                                                                             82 

C              ULL    IS    THE    LOWER    BIN    LIMIT  83 

UL     IS    THE    UPPER    BIN    LIMIT  84 

BIN    ■    0IM/80.0  86 

ULL    ■    0IM-I2MYY-1 )»BIN  87 

UL    -    ULL»2«B1N  88 


89 

90 


C      CYCLE  THROUGH  DATA  TO  F INO  WHICH  ONES  FALL  IN  »Y«  BINS                   91 

C  93 

IF  INST.GT.N)  GO  TO  7  ,                             '} 

C  95 

00  6  JJ-NST.N  '* 


IF     IMT4HJjl.lI.ULLI     GO    TO    7  97 

KST    ■    JJ  98 

AMAG  -  SQRTIOATAXI JJ I »OATAX  I  J J  I  »DAT AYIJJ I »OATAY IJJ I »  99 

C  100 

C               CHECK    THAT    MAGNITUDE     IS    NOT    OVER    DM  101 

C  102 

IF    (  AMAG.GT.OIM)    GO   TO   6  103 

C  10*. 

C                OK     IS     THE     FINAL     LINE     POSITION    FOR    THE    ■••  105 

C  106 

OK    •    OATAX|JJI*S*<>0.0/OIN»61.0  107 

IF    IOK.LT.10.0)     GO    TO    6  108 

K     ■     INHOK)  109 

K    -    iabsikI  110 

OK    «    AH  SI  OK)  111 

IF     IIOK-KI  .GT.0.5  I     K"K»1  Ii| 

IF     (.If  .  L  I  .  lO.O.O1*  .OK.GT.I  1  I  .01     GO    TO    6  113 

LINEIKI     «     ISYNI*)  IU 

6  CONTINUE  115 
C  116 

7  CONTINUE  1 17 
NST    '    KST»1  118 

C  119 

C                PRINT     OUT     ONE    LINE    OF     PLOT  120 

C  121 

WRITE     16,1-1     LINE  122 

8  CONTINUE  123 
C  12* 

RETURN  125 

C  126 

•9    FORMAT     M0X.1A4.'    COMPONENT    OF    THE    ELECTRIC    FIELD    IS  LESS'/IOX.                      127 

1     "THAN     l.E-6*i     THEREFORE     THIS    FIELD    WAS    NOT     •  /  10X  ,  •  PLOT  TEO.          EXEC  128 

2UTI0N    WILL    CONTINUE     AS    NORMAL. •//•  129 

10  FORMAT     UOX.'IHE     MAXIMUM    VALUE     OF    THE    8  1  STATIC    PATTERN    FOR    '/  130 

1  10X,  lA<,.«-<  ,1A<>,  •     I  INCIDENT-SCATTERED)     IS    LESS    THAN  •/                                              131 

2  IOX,    '     1.E-30.I          POLAR   PLOT    NOT    CALLE0.V//1  132 

11  FORMAT     I'l'.m,'     ELECTRIC    FIELD    ANTENNA    PATTERN    FOR  SPECIFIED    PIA               133 
1NE.« i 9X ."NORMAL  12 ING    FACTOR"     '.E10.5I  1  3A 

12  FORMAT     ('I8ISTATIC    SCATTERING    PATTERN    F OR • . 1  A* . • - • , I A* , • I  I NC 1  DENT-  135 
1SCAMERE0)     P0LARI2ATI0N.<  .9X  ,  •  NORMAL  I  2  I  NG    F  AC  TOR-  •  ,  E  10  .  5  »  136 

13  FORMAT     IMBACKSCATTERING    PATTERN    FOR  •  ,  I  A* ,  • -• ,  I AA , '  I  INC  IOENT- SCAT  T  137 
1ERE0I    POLARIZATION. '.9X. 'NORMALIZING    FAC TOR •■ t E 10 . 51  138 

I*    FORMAT     I1X.130A1 I  139 

ENO  1*0 
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PTPLOT 

PURPOSE:       to  establish  the  grid  information  for  the 
polar  plot. 

METHOD:        In  the  DO  LOOP  ending  at  statement  1  the 
alpha-numeric  characters  are  transferred  to  ISYN  in  order 
to  pass  via  COMMON  to  other  subroutines.  In  the  statements 
following  statement  2,  the  equations  for  the  plotted 
concentric  circles  are  established.  Below  statement  7  the 
grid  marks  on  the  090-270  axis  are  inserted. 

CALLED  BY:      POLPRT 

CALLS  TO:       LINECK 
NUMB 
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SUBROUTINE    PTPLOT     IIYY.Sl 

THIS    SUBROUTINE    SETS    UP    POLAR    GRID    INFORMATION 

COMMON     1SYM.LINE 

DIMENSION    L  I  NE  (  130),     ISYMI1AI,     I  SYNt  I  «•  I 

OAT  A     1SYN/IH»,1M.,1H     , 1H» , 1 H/, 1  HO, lHl , 1 H2 , 1H3, 1HA , 1 H5, 1 H6 , 1HB, 1H9/ 

INTtGER    Y.YY.W 

SET    UP    ISYM    FROM    I SYN    FOR    COMMON 

00    1    K'lil4 
ISYMIKI     -     ISYN1KI 
1    CONTINUE 


CLEAR    LINE    ANO    SET     TO    BLANK 


00  2     1=1,130 

2  LINEI  II    *    ISYMI31 

Y    -    *1- IYY 

IF     IY.EO.OI     GO    TO    7 

SET  UP  EQUATIONS  FOR  CONCENTRIC  CIRCLES 

YY  «  Y»Y 

1  '     I YY*2.5/2I*S 

X  =  6l.O»S0HT(2500.0-Z» 

CALL  LINECK  (X,YI 

IF  (Y.GT.32.0R.Y.LT.-32I  GO  TO  3 

X  -  61 .0»SU«T< 1600. O-ZI 

CALL  LINECK.  (X,YI 

3  IF  ( Y.GT.2<..OK.Y.LT.-2*l  GO  TO  * 
X  -  61.0»SQRT(900.0-ZI 

CALL  LINECK  IX, Yl 
*  IF  (Y.GT.  16.0R.Y.LT.-16)  GO  TO  5 
X  •  61 .0»SUKT (400.0-Z I 
CALL  IINCCK  IX, Yl 

5  IF  IY.GT.8.0R.Y.LT.-8I  GO  TO  6 
X  ■  61 .0»SQRT(lOO-2) 

CALL  LINECK  IX, Yl 

SET  UP  EUUATIONS  FOR  MULTIPLES  OF  30  DEGREES 

6  X  •  6l.OU.73205l»Y»S 
CALL  LINECK  (X.YI 

X  ■  61. 0»Y»S/l. 732051 


1 

2 

3 

A 

5 

6 

7 

8 

9 

10 

11 

12 

13 

1* 

\l 

17 
IB 
19 
20 
21 
ZZ 
23 
2* 
25 
26 
27 
28 
29 
30 
31 
32 
33 
3* 
35 
36 
37 
38 
39 
*0 

A*2 
A3 
A* 
♦  5 
A6 
AT 
A8 


CALL  L INECK  (X,Y) 

PUT  IN  POLAR  PLOT  NUMBER  LABELS 

CALL  NUMB  IYI 
W  •  I  Ad  SI  Y  I 

FILL  IN  POLAR  PLOT  AT  000,  090,  180,  AND  270 


10 


IF  (W.NE, 
L  INEI55I 
L  INE  157  1 
L INE 1591 
L  INE  I  63  I 
L INE 165  I 
L  INf  167  1 
IF  IW.NE. 
L INEI 561 
LINCI5B I 
LINEI60I 
L  INEI62I 

I  INI  If''-  I 

L INEI66 I 
IF  (W.NE. 
I  I  NE.  (  5  7  I 
I INt 159  1 
L INE (60  I 
I  INE  I  62  I 
L1NEI63I 
L  INEI65I 
IF  IM.NE. 
L  INEI58I 
LINEI60I 
LINEI62I 
L INEI6A I 
IF  (W.NE. 
LINE (59 | 
LINEI63I 
CUNT INUE 
RETURN 
END 


01  GO 
I  SYMI 
I  SYMI 
I  SYMI 
I  SYMI 
ISYMI 
1SYMI 

21  GO 
ISY1I 
ISYMI 
I  SYMI 
ISYMI 
ISYMI 
I  SYMI 

<.)  GO 
I  SYMI 
ISYMI 
ISYMI 
ISYMI 
I  SYM( 
ISYMI 

61     GO 
ISYMI 
ISYMI 
ISYMI 
ISYMI 

81  GO 
ISYMI 
ISYMI 


TO 

2 

2 

2 
2 

2 
2 

ru 

2 

2 

2 

2 

2 

2 

TO 

2  I 

21 

2 

21 

21 

21 

TO  11 

2, 

l\ 
TO  12 

ii 
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READ 


PURPOSE:       to  interpret  and  translate  the  input  data 
cards. 

METHOD:        The  program  utilizes  free  format  for  the  data 
cards,  that  is,  the  program  uses  character  recognition  to 
determine  which  parameters  are  being  read.    In  the  IF 
statements  containing  A(1),  A  (2),  A  (3)  ,  and  A  (4),  the  first 
four  characters  on  the  data  card  are  compared  to  the  first 
four  letters  of  the  key  words.  This  will  determine  the  type 
of  parameters  that  card  contains.   The  other  IF  statements 
determine  which  parameters  are  being  read. 

Subroutine  BLNK  is  called  to  remove  the  blank  spaces  on 
the  parameter  cards.   Subroutines  EQUAL  and  LEFT  are  called 
to  determine  the  position  of  the  equal  character  and  the 
left  paren,  respectively.  Subroutine  NUMBER  is  called  to 
convert  the  alpha-numeric  characters  to  numbers,  either 
fixed  or  floating  point.   This  numerical  value  is  assigned 
to  the  parameter  just  determined. 

A  detailed  explanation  of  the  data  cards  is  found  in 
appendix  II  titled  "USERS  MANUAL". 

CALLED  BY:      BLNK 
EQUAL 
LEFT 
NUMBER 
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SUBROUTINE    READ     MA,  IB, IBIS 
1 IFLAG,KFLAG,KGEN,LOAD,LZD,M 
2AFFP,AFF1,ABIP.AB1T,AM,BN,C 
3HII , PHSF ,PHSI , THAF . THAl.THl 
4,vniT,X,XNP,Y,YNP,2,ZLL0.7.N 
CIMMON    /A/    AIROI 
CCMPLEX     VULMl  I  ,ULDIl  I 
OIMENSIUN     I.M  I  I  ,      IBI  LI  .     XII 
DIMENSION    XNPI  1  I  ,     YNPI  I) ,     I 
DATA    AA, AB, AC, AD, AE.AF.AG.A 
1AM, AX/' A>  ,  «B' , "C ■  , •0* ,• E* , ' 
2'Q- .>R»  ,'S'  ,  •!•  ,  'U'  ,'W'  .'X- 
DATA    PLANK, COMMA, MINUS, PLEF 
•     •/  •/ 

5  7.2957  79 

4 

-1 

-1 

1 

■I 


V 

RAD    * 


C,  ICARO,  IGAIN, 1GRD,  I  NEAR,  INT, ISCAT.IWR, 
SG,NBAP,NB1P , NFFP.NGEN,NM,NP,ABAP, ABA  I, 
MM.ER2.EK3.ER4,FMC  ,HGI  , PHAF.PHAI  , PHIF.P 
f  ,  fill  I  ,THSF  .  THSI  , SlG2,SlG3,SIG4,TD2,T0  3 
P.STEPI 

.     Ylll,     2(11,     KGENI1I,     KFLAGI1I 
np I  1  I  ,    1201 II 

H,AI,iK,Al,AMA,AN,AO,AP,AQ,AR,AS.Ar,AU. 

F',•6•^,H•,■I,,•^,.,L•,lM•,,N1,,0•,,P•, 

/ 

T, POINT, RIGHT, SLANT/'     •,',•,•-•,•(',•.• 


INI 
IBISC 
IGAIN 
INEAR 
I  SCAT 
I  WR  -- 


I 


IF     I1FLAG.E0.6I     GO    TO    2 
IF     I MSG.NE.OI     GO     10    4 

1  READ    ( 5, 76,END=72)    A 

2  IF     I  (  A  I  I  I  .NE.ACI  .0R.IAI2I.NE.BLANKI.0R.C  A I  3>.NE.BL ANK )  .OR.  I  AI4I.NE 
1.3LANM )     GO    TO    3 

WRITE     16, 741    A 
GO    TO    1 

3  WRITE     I  6,751 

GO    TO    5  t 

4  READ    I5,76,EN0*72)    A 

5  IIARO    =     ICAROM 
WRITE      16.771     ICARO.A 

IF     I (MSG.NE.OI  .AND.  KAI1I.EQ.AEI  .AND. IAI  2)  .  EG.  AN)  .  AND.  I  A(  31  .EQ.ADI 
111     GO    TO    70 

IF     IIHSG.NE.OI.ANO.il  All  I .EC 
1.  AND.  I  A(4)  .EQ.API  I  I     GO    TO    6S 

IF     IIAI II .EQ.ACI .AND. IA(2I.( 
t.EO.BLANKI I    GO    TO    73 

IF    IHSG.GT.OI    GO    TO    4 

CALL    BLNK     IAI 

N    -    4 

INSULA) ION 


EQ.ASI.AND.tAI  2  I.FO.ATI  .AND.I Al  3I.EQ.A0I 

i9 

. EQ. BLANK  I. AN  0.1 AI3I.E0. BLANK)  .ANO.IAI 4 1 


IF     IIAI II 
1G0    TO    10 
KFLAGI20) 


NE.AI I .0R.IAI21.NE.ANI .OR. IAI 3) .NE.AS) . OR. I  A  I  4 ) .NE . AUI  I 


1 
2 
3 

5 

6 

7 

8 

9 

10 

II 

12 

13 

14 

15 

h 

18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
3* 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 


CALL    LEFT     INI 

6,  ,F.,!!*JN.,:NE:A'"  •°r-'*'n»h  .ne.aai.or.ia(n»2).ne.adi.or.<ain»3).ne 

1. A|)|  GO  TO  7 
KFLAGI4)  .  1 
CALL  EOUAL  IN) 
CALL  NUMBER  IN, N2, XI, IX) 
BM  .  XI 

IF  IAIN2I .EO.R IGHTI  GO  TO  4 
IF  I Al N2I .NE. SLANT )  GO  TO  71 
N  ■  N2»  1 
GO  TO  6 

7(  J^jj*!.*' jgf  j»D>.0R.t»IN»ll.NE.AII.0R.<AIN»i).NE.AE).0R.<»INO».NE 

KFIAGI6I  ■  I 

CALL  EQUAL  INI 

CALL  NUMBER  IN,N2,X1, 1X1 

ER2  -  XI 

IF    IAIN2I .EQ. RIGHT)    GO    TO    4 

IF     IAI  N2I  .NE.SL  ANT  I    GO    TO    71 

N    •    N2» 1 

GO    TO   6 

8.  '^     !!4'^,;^-*c,'OK-'*<N»H.NE.AOI.OR.IAIN»21.NE.ANI.OR.lA(N»3I.NE 
l.AOM     GO    TO    9 
KFLAGI5I     -    1 
CALL    EjUAL     IN) 
CALL    NUMBER     IN, N2, XI, IX) 
SIG2     *     XI 

IF    IAIN2I .E0.R1GHTI    GO    TO    4 
IE    IAIN2I.NE. SLANT)    GO    TO    71 
N    *    N2»  1 
GO    TO    6 

9.'L  !!  *IN>;NE.AL).OR.(AIN»l).N6.AO).OR.I»IN»2l  .NE.AS)  .OR.  IAIN*  3).  NE 

1  •  *  S  I  I     (..()     TO     71 
KFLAGI7)     -     | 
CALL    EQUAL     INI 
CALL    NUMBER    IN, N2, XI,  IX) 

T02    '    XI  « 

IF    IAIN2I.EQ. RIGHT)    GO    TO   4 
IF    IAIN2I.NE. SLANT)    GO    TO    Tl 
N    •   N2tl 
GO    TO   6 

WIRE  » 


49 
■50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 

72 
73 
74 
75 

?? 

78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 


94 
95 
96 


74 


10 


IF    I  I  At  D.NE.AWI  .OR.  I  A I  2 ) . NE . A  I ) .OR. 1 A< 3 ) .NE .AR) .OR. < A  I  *  I .NE .AE ) » 

go  ro   13 

CALL    LEFT     (Nl 


II 
1 


IF     l(AIN).NE.AR)  .  OR.  I  AiNM  I  .NE.AA).0R.IA<N»2).NE.AD).OR.<AIN»3I.NE 

.All)     GO    TO    12 

KFLAGI2)     =■     1 

CALL    EQUAL     (N) 

CALL    MUMPER     (N,N2,X1,IXI 

AM    =     XI 

IF    (AIN2I .FO.R IGHT)    GO    TO    * 

IF     IAIN2).  NE.  SLANT  1    GO    TO    71 

N    «    N2«  1 

GO    TO    1 1 

12  IF     I  I  A  IN)  .NE.ACI .OR. I  A  IN* 11 .NE.AOI.OR.I  A  I N»2 > .NE . AN) .OR. ( A  IN » 3) .NE 
1.A01I     GO    TO    71 

KFLAGI3)     »     1 

CALL    ECUAL     IN) 

CALL     NUMBER     IN, N2, XI,  IX) 

Q  h«*«        -         %     1 

IF    <A(N2I .EO.P IGHTI    GO    TO    * 
IF     IAIN2I.NE. SLANT)    GO    TO    71 
N    =    N2»l 
GO    TO    1 1 

EXTERNAL    MEDIUM 

13  IF     I(AI1I.N£.AE).0R.IAI2).NE.AXI .OR.I A(  3 ) .NE. AT ) .OR. ( A ( *  I .NE . AE II 
1G0    TO    1 7 

K.FLAGI8)    ■=    1 
CALL    LEFT     IN) 

1<.     IF     I  I  AINI.NE  .AC)  .OR.  I  A  IN*  1)  .NE.AOI.OR.I  A  I  N>2 1  .  NE  .  AN  I  .0R.IA(N»3).NE 
1  -Ad) I    GO    TO    15 
KFLAGI9)     ■=     1 
CALL    ECUAL     INI 
CALL    NUMBER     IN, N2, XI, IX) 
SIG3    «     XI 


IF    IAIN2I.E0.R  IGHT  I    GO    TO   * 
SLANTI 


GO    TO    71 


IF    (AIN2I.NE 
N    -    N2«l 
GO    TO    1* 

15     IF     <<AINI.N(.A3).QR.<AINM).Ne.AI).0R.(AIN»2I.NE.AE).0R.lAIN«3).NE 
I. ALII     GO    TO    16 
KFLAGI10I     ■     1 
CALL    EQUAL     INI 


97 
98 
99 
100 
101 
102 
103 
10* 
105 
106 
107 
108 
109 
110 
111 
112 
113 
11* 
115 
116 
117 
lid 
119 
120 
121 
122 
123 
12* 
125 
126 

\27 
128 
129 
130 
131 
132 
133 
13* 
135 
136 
137 
138 
139 
1*0 
1*1 
1*2 
1*3 
1** 


CALL     NUMBER     IN, N2, XI, IX) 

ER3    '    XI 

IF    <A(N2).EQ.RICHTI    GO    TO 

IF     IA(N2).NE. SLANT)    GO    TO 

N    •    N2» 1 

GO    TO    I* 


S 

71 


16 


.!ASIl*GO,TOEflU'OR"<*IN*"'NE,*0,'OR',*'N*2KNE-*S,*OR-'*(N*3,*NE 
"XFI AG( 111-1 

CALL     EUUAL     IN) 

CALL    NUMBER    (N,N2,X1,  IX) 

TD3    -     XI 

IF     ( AIN2) .EO.R IGHT )    GO    TO    * 

IF    (AIN2I .NE. SLANT)    GO    TO    71 

N    .    N2» 1 

GO    TO    1* 


LOAD 


.17 

18 

I 

19 
20 


GO    }u4{V*NE'*L''OR'<*'2,'NE-*0,-OR'<*,3,-NE'**'-OR-,*«*,-NE'*0" 


KFLAGIl*)     •     1 

GO    TO     19 

IF     I IAI 1) .NE.AII 

GO    TO    22 
KFLAGI2*)     •     1 

J     '     I 

CALL     LFFT 


OR 


<A<2).NE.AHA).OR.IAO).NE.API.OR.(AC*).NE.AE)I 


(Nl 


CALL    NUMBEP     IN, N2, XI, IX) 


5' 


TO    2 


IF    IIX.LE.O) 

LIUU  I     .     IX 

N    ■    N2» 1 

CALL    NUMBER     IN, N2, XI, IX) 

RMAG    «     XI 

N    •    N? •  1 

CALL    NUMBER    IN, N2, XI, IX) 

RUE1",    •    XI 

RREAL    -    RMAG*COS(ROEG/RAD| 

RIMAG    •    RMAG»SINIROEG/RAOl 

7LL0III     -    CMPLX(RREAL.R)MAG) 

LOAO    ■     1 

IF     IAIN2I.EQ.R IGHT I    GO    TO    * 

IF    <A(N2).NE. SLANT)    GO    TO    71 

1    -    1*1 

N    -    N2»l 

GO    TO    20 


1*5 
1*6 
1*7 
1*8 
1*9 
150 
151 
1^2 
153 
15* 
155 
156 
157 
158 
159 
160 
161 
162 
163 
16* 
165 
166 
167 
168 
169 
170 

hi 

17* 
1  75 
176 
177 
178 
179 
180 
181 
182 
(83 
18* 
165 
186 
187 
188 
189 
190 
191 
192 
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21  KFLAGI2*)     -    -1 
LOAD    "    -1 

GO    TO    * 

FREQUENCY 

22  IF     <l4ll).NE.AFI.OR.iAI2I.NE.ARI.OP.(A(31.NE.AE).OR.IA(*).NE.Aal) 
l&O    TO    23 

KFIAGI1)     -     1 

CALL    LETT     INI 

CALL    NUMBER     <N,N2,XI,IX) 

FMC    •     XI 

GO    TO    * 

PLOT 

23  IF     ((A(l).NE.API  .OR.  (412I.NE.AL  I. OR. 1*131. NE.*OI. OR. <*!*). NE.ATII 
1GO    TO    31 

KFLAG<22I     »    I 
CALL    LEFT     (Nl 

2*     IF     UAINI.NE.AFI.OR.IAtN»|I.NE.AA).OR.<AIN»2).NE.AR).OR.IA<N»3).NE 
l.AFII    GO    TO    25 
IGAIN    •     1 
NFFP    •     I 
GO    TO    27 

25  IF     l(4(M).NE.*B).0R.I41N'll.NE.AI  I  .OR.  I  A  I  N»  2 )  .NE  .  AS  I  .OR  .  I  A I N  ♦  3)  .  NE 
I. ATI  I     GO    TO    26 

IBISC    ■     I 
NBIP    ■     1 
GO    TO    27 

26  IF     IIAlNI.NF.AR).0R.IAINMI.NE.AAI.0R.IA(N.2).Nl-.ACI.0R.IA(NOI.Nf 
l.AKII     GO    TO    71 

I  SCAT     -     1 
NBAP    -     1 


27  DO   28    I -N.80 
K    •    1*1 

IF    (All  I. E8. SLANT  I    GO    TO    29 

28  CONTINUE 


GO   TO    7  1 
29   N    ■   K 

IF     (IAIN) .NE.ATI  .OR. IAIN*  I I .NE.AHI.OR.I AINf2l.NF.AE) .OR . ( A( N» 3  I . NE 


193 
19* 
195 
196 
197 
198 
199 
200 
201 
202 
203 
20* 
205 
206 
207 
208 
209 
210 
211 
212 
213 
21* 
215 
216 
217 
218 
219 
220 
221 
222 
223 
2'* 
225 
226 
227 
228 
22H 
2  30 
231 
232 
233 
23* 
235 
2  36 
237 
218 
239 
2*0 


30 


31 


I  .41)  I     GO    TO    10 
CAll     EQUAL     (N) 
C4LI     NUMBER     <N,N2,Xl , IX) 
IF     (NFF  P.EO. 1  I     AFF  T-X  I 
IF     INBI  P.  tO.  II     Adl  T*X  1 
IF     ltl84P.E0.il     ABAT'Xl 
IF     ( AIN2I .EQ. RIGHT )    GO    TO    * 
IF     (AIN2I.NE. SLANT)    GO    TO    71 
U    ■    N2«l 

GO    TO    2*  j 

IF     I  lAINI  .NE.API  .OR.I  A(N»1  >.NFJ.AH).OR.I  A(N»2I.NE.A1  ))     GO    TO    71 
CALL    EOUAL     INI 
CAll     NUMBER     IN.NJ.XI  ,  IX) 
If     (NFFP. EQ. II     AFFP-X1 
IF     INRIP.EU.l I     ARIP-X1 
IF     INB4P.Eg.ll     ABAP-Xl 
IF    (AIN2I .EU.R IGHTI    GO    TO    * 
IF     IAIN2) .NE. SLANT  I    GO    TO    71 
N    «    N2» I 
GO    10    2* 

OUTPUT 

IF     II Al  ll.NE.AO)  .OR.I A( 2) .NE . AU I .OR. I A(  3  I .NE. AT  I . OR. I  A  I  *  I .NE . API  I 
1G0    TO    ** 
KFLA&I22I     >    1 
CALL    LEFT    IN) 


32    IF    KAINI.NE.ABI  .OR.I  AINHI.NE.A1  )  .OR .  I  A  I  N»  21  .NE  .  AS  1  .OR .  <  A  <  N»  31 .  NE 
l.AI) I    GO    TO    33 


KFLAGI1BI     ■     1 

I 8ISC    •     1 

CALL    ECUAL     INI 

CALL    NUMBER     IN,N2.X1,|XI 

PHSI     ■     XI 

N    -    N2»  I 

C4LL     NUMBER     IN, N2, XI, IX) 

PHSF     .     XI 

N    -    N2»  I 

CALL    NUMBER     IN.N2.XI , IXI 

TMSI    ■    XI 

N    -    N2» I 

CALL  NUMBER  IN.N2.X1, IXI 

THSF  ■  XI 

IF  (»<N2I .EO. RIGHT)  GO  TO  * 

IF  (AIN2) .NE. SLANT )  GO  TO  71 

N  ■  N2»  I 

GO  TO  32 


2*1 
2>,2 
2*3 
2** 
2*5 
2*6 
2*7 
2*8 
2*9 
250 
251 
252 
253 
25* 
255 
256 
257 
2  58 
2  59 
260 
261 
262 
263 
26* 
265 
266 
267 
268 
269 
270 
271 
272 
273 
27* 
275 
2  76 

iJ1 
278 
279 
280 
281 
282 
283 
28* 
285 
286 
?87 
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N    -    N? 

37 

CONTINUE 

c 

c 

c 

3a 

GO    10    7  1 

N2    ■    N2»l 

INEAR     =     1NEAR-1 

39 

IF    (AIN2I  .EQ. RIGHT) 

GO 

TO    * 

IF    ( A( N21 .NE. SLANT  1 

GO 

TO    71 

N    ■    N2»  1 

GO    TC   32 

33l!AF)5*Go'TOE3SF,,OR,,*,Nfl,'NE,M,,OR-,*(N*2,-NE-*R,-OR-l*,N*3,-NE  "o 

KFLAGI16I     ■     1  "I 

IGAIN    =     I  292 

CALL    ECUAL     (N)  ?93 

PMAI    NU;?ER    <N.N2.X1,IX)  f^ 

N    *    N2M  |96 

CALL    NUMBER    IN.N2.  XI.  IXI  ™J 

N    »    N2H  "9 

CALL    NUMBER     IN. N2. XI, IX)  300 

N    =    N2»I  |£2 

THAF    ""xi"    •"•«»'«»»««•  30* 

IF    IAIN2I.EQ. RIGHT)    GO    TO   *  12? 

IF    IAIN2I.NE. SLANT)    GO    TO    71  f0,6. 

c             GO*tBU  308 

3*UABI)*GO,foE;oN,"OR-,*,N*l,-NE'AE,'OR'l*<N*i!,-NE-*A,-0,,-,*tN»3,-NE  31» 

KFLAGI19I    -    1  3|2 

INEAR  x  2  313 

CALL  EQUAL  INI  I* 

IF  (A(N).EQ.PLEFT)  GO  TO  35  If? 

INEAR  »  1  | J* 

1*1  317 

GO  TO  36  318 

C  319 

C  320 

C  321 

35  00  37  L«l,50  322 
I  ■  L  323 
N  =  N»l  US 

36  CALL  NUMBER  IN, N2, XI, IXI  HI 
XNPI1  I  -  XI  326 
N  -  N2»l  |27 
CALL  NUMBER  IN,N2fXl.IXI  l|S 
YNPIII  -  XI  1|9 
N  ■  N2»l  |30 
CALL  NUMBER  IN, N2, XI, IXI  |3l 
ZNPIII  .  XI  "2 
IF  I  INEAR. EO.U  GO  TO  39  HI 
INEAR  -  L»l  |3* 
IF  (AIN2I.EQ. RIGHT)  GO  TO  38  \H 


3  37 
338 
339 
3*0 
3*1 
3*2 
3*3 
3** 
3*5 
3*6 
3*7 
,  3*8 

*  3*9 

*0  IF  (  (  A(  Nl  .NE  .AB)  .OR.  (  AINU  I  .NE.AA)  .0R.(  A(N»2(  .NE.AU.OR.  (AIN»31  .NE       350 
1 .AK)I  GO  TO  *1  351 

KFLAGI17I  •  1  352 

ISCAT  =  1  353 

CALL  E3UAL  (N)  35* 

CALL  NUMBER  IN, N2, XI, IXI  355 

PH1I  -  XI  356 

N  «  N2f 1  357 

CALL  NUMBER  IN, N2, XI, IX)  35B 

PHIF  '    XI  359 

N  =  N2»l  360 

CALL  NUMBER  IN, N2, XI, IXI  361 

THII  •  XI  362 

N  •  N2» 1  363 

CALL  NUMBER  IN, N2, XI, IXI  36* 

THIF  -  XI  365 

IF  (AIN2) .EQ. RIGHT)  GO  TO  *  366 

IF  (A< N2) .NE. SLANT  I  GO  TO  71  367 

N  ■  N2«  1  368 

GO  TO  32  369 

3  70 

*1  IF  I  IAIN). NE.ACI.OR.IAINHI.NE.AU). OR. l»IN»2I.NE.ARI. OR. <»IN»3).NE      371 

I. ARM  GO  TO  *3  372 

KFLAGI 151-1  373 

I  up.  •  1  3  7* 

375 

3  76 

377 

00  *2  H-N. 80  378 

IF  (AIK I.EO.R  IGHT  1  GO  TO  *  3  79 

N  ■  KM  300 

IF  (AUI.EO.  SLANT)  GO  TO  32  381 

*2  CONTINUE  382 

383 
GO  TO  71  38* 


77 


C  385 

43  IF  (IAIN). NE. AS)  .OR. I  A IN»1 ) .NE . AT ) .OR. ( A (N»2) .NE. AE) .OR. ( A(N» 3) .NE  386 
I. API)  GO  TO  Tl  387 

CALL  EQUAL  (Nl  388 

CALL  NUMBER  IN, N2, XI, IX)  389 

SIEP  -  XI  390 

IF    (AIN2I.EQ. RIGHT)    GO    TO    4  391 

IF    (AIN2I.NE. SLANT)    GO    TO    71  392 

N    =    N2M  393 

GO    TO    32  39<i 

C  395 

C  FEED    POINT  396 

C  39/ 

44  IF    II  A(  D.NE.AF)  .OR.IA(2I.NE.AEI.OR.IA(3).NE.AE).OR.1AI4I.NE.ADI>  398 
1G0    TO    45  399 

KFLAGI13)     •     1  400 

GO    TO    46  401 

45  IF    I  I A( II .NE.AG) .OR.I AI2) .NE.AEI .OR.IAI 3>.NE.AN>.OR.(A(4I.NE.AE) )  402 
1G0    TO    49  403 

KFLAGI23)     *■     1  404 

46  NGEN  »  0  405 
CALL    LEFt     (Nl  406 

47  CALL  NUMBER  IN, N2, XI, 1X1  407 
NGEN  =  NGENM  408 
KGENINGENI  =  IX  409 
IF  IAIN2I .EQ.RIiHTI  GO  TO  *  410 
N  =  N2» I  411 
CALL  NUMBER  IN, N2, XI,  IX)  412 
VMAG  ■  XI  413 
N  =  N2» 1  4|4 
CALL  NUMBER  IN, N2, XI, IX)  415 
VOEG  ■  XI  416 
VREAL  ■  VMAG»COSIV0EG/RAO)  417 
VIMAG  =  VMAG'SINI VOEG/RAOI  418 
VOLTINGEN)  ■=  CMPLX  (  VRE  AL  i  V  I  MAG  )  419 
IF  (AIN2I.EQ. RIGHT)  GO  TO  4  420 
IF  (AIN2I.NE. SLANT)  GO  TO  71  421 
IF  l(A(N2l.  EQ.  SLANT)  .AND.  I  A  (N2  ♦  I )  .EQ.  BLANK)  I  GO  TO  48  422 
N  -  N2» I •  423 
GO    TO    47  424 

48  REAO  (5,761  A  425 
ICARO  '  ICARDM  426 
WRITE  (6,77)  ICARO, A  427 
N  ■  1  428 
CALL  BLNK  IAI  429 
GO    TO    47  430 

C  431 

C  432 


C              DESCRIPTION  433 

49  IF    KAI  11  .NE.AOI  .OR.(A(2I.NE.AEI.OR.IA(3).NE.AS)  .0R.IAI4).NE.ACI)  435 
I  GO    10    5  2  lu 

KFLAGI12I     ■    1  437 

J    ■    0  438 

CALL     LEFT     INI  439 

50  CALL  NUMBER  IN, N2, XI,  IX)  4411 
J  *  J»l  441 
NM  =  J  44} 
IAUI  '  IX  44  i 
N  =  N2»l  444 
CALL  NUMBER  IN, N2,  XI,  IX)  us 
IBIJI  =  IX  <,<,«, 
IF  (AIN2I.EQ.RIGHTI  GO  TO  4  447 
IF  (AIN2)  .NF. SLANT)  GO  TO  71  4  4  11 
IF  (IAIN2I.EQ. SLANT). ANO. I  A t N» 1 1 . EQ. BLANK!)  GO  TO  51  449 
N  -  N2»l  450 
GC    TO    50  451 

51  REAO  (5,761  A  45J 
ICARD  -  ICARDM  451 
CALL  BLNK  (A)  4,*,* 
HR1TE  (6,77)  ICARO, A  455 
N  «  1  456 
GO    TO    50  457 

C  458 

C                G6CMETRY  459 

52  IF     KAI  ll.NE.AG)  .OP. I  A(  2)  .NE.  AE  I  .OR.  I  A(  3  )  .NE  .  AOI  .OR.  (  A(  4  I  .NE  .  ANA)  )  46? 
1     GO    TO    55  462 

KFLAGII2)     -     1  463 

JJ    ■    0  464 

CALL    LEFT     IN)  465 

53  CALL  NUMBER  IN, N2, XI, IX)  466 
J  J  "  JJ»1  467 
NP  -  JJ  468 
XIJJ)  -  XI  469 
N  -  N2»l  470 
CALL  NUMBER  (N,N2,X1,IX)  Wl 
YIJJ)  •  XI  472 
N  ■  N2» I  473 
CALL  NUMBER  IN, N2, XI, IX)  474 
2IJJ)  «  XI  475 
IF  (AIN2I. EO. RIGHT  )  GO  TO  4  476 
IF  (AlN2l.NE. SLANT)  GO  TO  71  477 
IF    (IAIN2I.EQ.SLANTI. ANO. (AINU). EO.BLANKII     GO    TO    54  478 

.        N    •    N2»l  479 

GO   TO    53  480 
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5*    READ    (5,76)    A  481 

ICARD    «      ICARD*!  482 

WRITE     16,771      ICARD, A  483 

CALL    BLNK     (A)  4B4 

N    -     1  4d5 

GO    TO    53  4B6 

C  487 

C               INTERVAL    FOR    CALCULATION  488 

C  489 

55  IF  (I  A(  1)  .NE.Al)  ,OR.(  AI2I.NE.AN)  .OR.IAI  31  .NE.  AT  I  .OR.  I A  ( *  )  .NE.AE))  490 
ISO  TO  56  491 

KFLAGI21I  «  t  492 

CALL  LEFT  IN)  493 

CALL    NUMBER    IN,  N2,  XI,  1X1  494 

INT    «     IX  495 

IF    IA(N2).EQ.RIGHTI    GO    TO  *                                                                                                                 496 

GO    TO    71  497 

C  498 

C                GROUND  499 

C  500 

56  IF     I  (  A(  I)  .NE.AG)  .0R.(  A(2)  .NE.AR)  .OR.IAI  3l.Nt.AOI  .OR.  (  A  I  *  I  .NE  .  AU)  )  501 
1GU    TO    66  502 

KFLAGI25)    »    1  .503 

KFLAGI26)    ■    1  50* 

IGRD    *    2  505 

CALL    LEFT     IN)  506 

57  IF    I1AIN1 .NE.API .OR.IAfNMI .NE.AE) .0R.(A(N»2) .NE- .  AR  I  .OR.  (  A  (N»  3)  .NE  507 
' l.AFI)    GO    TO    58  508 

1GRD    *    1  509 

GO    TO    64  510 

58  IF     MAIN)  .NF.Ail  .0R.(A(N*1).NE.A0).0P..(  A(N»2).NF.AOI.OR.(AlN«3).NE  511 
l.ADI I    GO   TO   59  512 

ER4    ■     30.  513 

SIG4    -     .02  514 

GO    TO   64  515 

59  IF     I  I  Al  Nl  .ME  .API  .OR. I  AINU)  .NE.A01.0R.I  AIN*2I  .NE.AOI  -OR. I  A(N«3).NE  516 
l.ARI)     GO    TO    60  51  7 

ER4    «    4.  518 

SIG4    -     .001  519 

GO   TO   64  520 

60  IF  ((A(N)  .NE.ASI  .OR. (AINU)  . NE.AE). OR. (  AIN*2)  .NE.AA)  )  GO  TO  61  521 
ER4  •  80.  522 
SIG4  -  4.  523 
GO    10    64  524 

61  IF     I (Al Nl .NE.AH) .0R.IA(N»1)  .NE. AE I .OR. ( A (N* 2  I .NE . A  I I .OR . ( A ( N » 3  I . NE  525 
i  .  Ar.i  i    r.n    rr>  fc?  *  j*. 


LAG))    GO    TO  62  526 

28 


CALL    EOUAL    (Nl  527 

CALL    NUMBER    IN, N2, XI, IX)  52 


HGT    -    XI  529 

IF    IA(N2I .EO.R IGHTI    GO    TO  4                                                                                                                          530 

IF    IAIN2I.NE. SLANT)    GO    TO  71                                                                                                               531 

N    «    N?* 1  532 

GO    TO    57  533 

62  IF     I  IAIN)  .Nt .AC)  .OR.IAIN* 1)  .NE.AOI. OR. ( A  I N» 2 ) . NE . AN ) .OR . I  A  I N* 31 . NE  534 
l.ADI)     GO    TO    63  535 

CALL    EQUAL    (N)  536 

CALL    NUMBER     IN, HI, XI, IX)  537 

SIG4    <=     XI  5  38 

IF     (A(N2).EQ.RIGHTI    GO    TO  4                                                                                                                                    539 

IF    I Al N2l .NE. SLANT)    GO    TO  71                                                                                                               540 

N    ■    N2*l  541 

GO    TO    5  7  5.2 

63  IF     (I AINi .NE.A3I .0R.(A(N*1I .NE.AI1.0R.I  A(N*2I .  NE  .  AE  I -OR  .  (  A  (  N*  31  .  N  t  543 
I. AD)     GO    TO    71  544 

CALL    EQUAL     IN)  545 

CALL    NUMBER     (N,N2,X1,IXI  546 

ER4    "    XI  547 

IF  (A(N2) .EQ. RIGHT)  GO  TO  4                                              548 

IF  (AIN2) .NE. SLANT)  GO  TO  71                                                    549 

N  '    N2U  550 

GO  TO  57  551 

C  552 

C  553 

C  554 

64  00  65  K=N,80  555 
IF  lAIKI.EO. RIGHT)  GO  TO  4  556 
N  «  KU  557 
IF  (A(K).EQ.SLANTI  GO  TO  57  558 

65  CONTINUE  559 
C  560 
C  561 
C  562 

GO   TO    7  1  563 

C  564 

C  565 

66  IF     I ( Al 1) .NE .AS1 .OR. ( A(2I -NE. ATI .OR.IAI 3). NE.AOI .OR. (A(4I. NE.API I  566 
1G0    TO    67  567 

IFLAG    ■    2  568 

RETURN  569 

C  510 

67  IF    (I Al  II .NE.AC) .0R.IAI2) .NE.AH) .OR.(A(  3 ) .NE . AA ) . OR . ( A ( 4  I .NE . AN) I  571 
IGO    TO    68  572 

IFLAG    -    3  573 

RETURN  574 


68     IF     MAI  II. NE.AE). OK. (  A(2I.NE.AN)  .OR.IAI  3). NE.AOI  I     GOTO     71  5  76 


W 


IFLAG    •     1  5,B 

RETURN  570 

69  IFLAG  •  5  S90 
RETURN  ijHI 

70  lFLAl,  ■  4  51I2 
RETURN  ia\ 

71  MSG  »  I  5H<, 
KFLAG130)     ■     ICARD  i,ub 

72  If  UFLAG. NE. 51  MRITE  (6,781  jjjj'j 
IFLAG  ■  5  5B8 
RETURN  58v 

C  590 

73  IFLAG  ■  6  iqi 
ICARD  -  ICARO-l  5<j2 
RETURN  59J 

C  594 

C  595 

74  FORMAT     (5X.80A1 I  507 

75  FORMAT     (////5X,'oATA    CARDS'//!  IZL 

76  FORMAT     I80A)  I  .,.,., 

II   fBShat    If'^lStrtU   CARD/STOP   CARO   MISSING 600 

eno  79 


RITE 


PURPOSE:       to  generate  a  list  of  branch  currents  from 
the  input  loop  currents. 

METHOD:        The  generation  of  branch  currents  is 
accomplished  in  the  DO  LOOP  ending  at  statement  2.  The 
branch  currents  are  stored  in  CJ  (I)  by  the   latter  part  of 
the  DO  LOOP  ending  at  statement  3.  If  the  branch  currents 
are  reguested  for  output  (IWR  positive) ,  the   DO  LOOP  ending 
at  statement  5  accomplishes  this. 

CALLED  BY:      GANT1 
GFFLD 

CALLS  TO:       NONE 
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c 

c 

00  3  K=l,NM 

KA  ■  14  IK) 

KB  •  I6IK) 

CJA  -  1 .0, .01 

CJB  *  1  .0,.0) 

NDK  »  NOIK) 

c 

c 

SUBROUTINE  KITE  I  I  A , I B ,  INM,  IU11,  II , I  2 i I  3 • HO. ND.NM.C Ji CGi  1GR0)  1 

COMPLEX  CJIll ,CGI1 I.CJA.CJB  2 

DIMENSION  lAlil,  16(11,  11(1).  1211).  13(1).  ND(INN,4),  NDIll  3 

A MAX  -  .0  * 

5 

6 

7 

8 

9 

10 

11 

12 

13 

1* 

00  2    II  -=l,NDK  15 

1  •  MO  I K , I  1  I  16 
Fl  '  I.  17 
IF  (KB. FQ. 121 II)  GO  TO  1  18 
IF  (KB. EQ. Ill  II  I  Fl  — 1.  19 
CJA  =  CJA»FI»CJt 1 »  20 
GO    TO    2                                                         .                                                                                                             21 

1  IF  (KA.  to.  131  II  )  Fl  — 1.  22 
CJB    *    CJBfFI*CJI II  23 

2  CONTINUE  2* 
C  25 
C  26 

CGIKI    =    CJA  27 

KK    =    K.NM  28 

CGIKK)     »    CJB  29 

ACJ    '     CABS(CJA)  30 

BCJ    =    CABSICJBI  31 

IF     IACJ.GT.AMAX I     ANAX»ACJ  32 

IF     (BCJ.GT.AMAX)     AMAX-BCJ  33 

3  CONTINUE  3* 
C  35 
C                                                                                                                                                                                                     36 

IF     I  IWR.GT.OI    GO    TO    *  37 

RETURN  38 

4  IF  IAMAX.LE.O.)  AMAX-1.  39 
MR  I TE  (6.81  *0 
NMG  »  NM  41 
IF     UGRD.GT.O)     NMG    ■    NM/2                                                                                                                                             42 

C  43 

00    5    K=l .NMG  ** 

CJA   ■    CGlKI  45 

KK    •    K»NM  46 

CJB    -    CGIKK)  47 

CCJA    •    CABS(CJA)  48 


CCJB    -    CABSICJBI  49 

ACJ    ■=    CCJA/AMAX  50 

BCJ    '    CCJB/AMAX  51 

PA    «    .0  52 

PB    *    .0  53 

IF     IACJ.GT.O.)     PA     -    57.29578*ATAN2(AIMAGICJAI .REAL(CJA) I  54 

IF    IBCJ.GT.O.)    PB    -    57.295/8»ATAN2(AIMAG(CJB>,REAL(CJB) I  55 

5  WRITE     (6,71     K,IA(K)  ,C JA.CCJA.AC J.PA.IB(K)  .CJB.CCJB.BCJ.PB  56 

C  58 

WRITE     (6,6)  59 

RETURN  60 

C  61 

C  62 

6  FORMAT     I IHOI  63 
/    FORMAT    (2X,I2,2(2X,I2,2X,EI  1.5.1X.E11.5,  IX,  El  1 .5  ,  1X.EH  •  5,  IX,  F6. 1 )  64 

II  65 

8  FORMAT  I/2I46X, 'NORMALIZED'  ,5X  )/•  SEC', 21*  NODE' , 4X , 'RE AL • ,6X , • I M»  66 

IG INAPT'  , 3X, 'MAiNI  TUOE' , 3X , < MAGN I TUOE '  , 3X, 'PHASE' J )  67 

ENO  68 


81 


SART 


PURPOSE:        to  sort  data  for  polar  plot. 

METHOD:        This  subroutine  sorts  the  values  of  the 

points  to  be  plotted  by  the  polar  plot  package  starting  with 

the  greatest  positive  value  of  y  to  the  greatest  negative 

value.   In  the  DO  LOOP  ending  at  statement  1,  the  value  of 

(x  ,y  )  is  interchanged  with  the  value  of  (x  ,y  )  if  y   is 
i   i  J   J       J 

greater  than  y  . 
i 

CALLED  BY:      POLPRT 


CALLS  TO:        NONE 


SUHBOUUNf      ',*»!       IHAI AX.PAI AY.NI  1 

DIMENSION    DATAXI500),  DATAYI500)                                                                                                                           J 

C  3 

C       THIS  ROUTINE  SORTS  OATA  IN  DATAY  BY  MAGN1TU0E  * 

C  5 

NN    .    N-l  (, 

C  7 

00    2     !• 1,NN  8 

NM    -     I.I  9 

C  10 

00    I     J=NM.N  |l 

If     IDATAYUI  .GE.OATAY(JI)    CO    TO    1  12 

STOft  ■  DAT  AY  I  I  I  13 

OATA  Yl II  •  OATAYIJI  |« 

OATA  Yl J)  ■  STOR  15 

STOR  -  OATAM  I  )  1  fc 

OATA  XII)  •  OATAXIJI  17 

OATA  XI Jl  •  SIOA  18 

1  CONTINUE  19 
c  20 

2  CONTINUE  21 

c  22 

RETURN  3} 

ENO  24 
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SGANT 


PURPOSE:       to   calculate  the  mutual  impedance  between 
filamentary  monopoles. 


2   •                                                  f    2 

A 

I' 

B 
1    , 

SOUI 
DIP( 

%*                         ^     c 

P21*                   \ 

I  1 

*CE                                  RECEIVING 
)LE                                       DIPOLE 

Z  =  P11  +  P12  +  P21  +  P22 


METHOD:        In  the  induced  emf  formulation,  the  mutual 
impedance  of  coupled  dipoles  is 


Z  =  - 


/ 


I  (t)  E  (t)  dt 

2        1 


where   I    (t)    denotes  the  current  distribution    (normalized   to 
2 

unit  terminal  current)    on   dipole   2,    and   E    (t)    is   the   field 

1 

of  dipole  1  when  it  transmits  with  unit  terminal  current. 

Distance  along  the  axis  of  dipole  2  is  denoted  by  the 

coordinate  t.  E   may  be  expressed  as  the  sum  of  the  fields 

1 

from  each  of  the  monopoles  comprising  dipole  1.  Furthermore, 

the  integral  is  the  sum  of  the  integrations  over  each  of  the 

monopoles  comprising  dipole  2.   Thus,  the  dipole-dipole 

mutual  impedance  may  be  expressed  as  the  sum  of  four 

monopole-monopole  impedances. 

It  may  be  convenient  to  draw  the  above  figure  in  terms 
of  monopoles  with  the  current  distribution  shown  as  dotted 
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lines.   (The  monopole  letters  remain  the  same.) 


P11 


P12 


B 


P21  / 


P22 


D 


The  surface  impedance  is  calculated  just  above  statement 

2.  B01  denotes  J   /J   where  J   and  J   are  the  Bessel 

0     1        0      1 

functions  of  order  zero  and  one  with  complex  argument,  ZARG. 

It  is  assumed  that  all  the  wire  segments  have  the  same 

radius,  conductivity  and  surface  impedance. 


In  the  DO  LOOP  ending  with  statement  3,  SGANT  calculates 
the  segment  lengths  D ( J)  .  DHIN  and  DMAX  denote  the  lengths 
of  the  shortest  and  longest  segments.  If  the  wire  radius  or 
the  segment  lengths  are  clearly  beyond  the  range  of 
thin-wire  theory,  N  is  set  to  zero  at  statement  4  followed 
by  RETURN  to  the  main  program  to  abort  the  calculation. 

At  statement  5,  the  program  selects  a  segment  K,  and  a 
few  statements  below  this  it  selects  another  segment  L.  K  is 
a  segment  of  test  diple  I,  and  L  is  a  segment  of  expansion 
mode  J.  The  mutual  impedance  between  segments  K  and  L  is 
obtained  by  calling  subroutine  GGS  or  GGMM.  In  statement  18, 
this  impedance  is  lumped  into  C(MMM).  The  mutual  impedance 

Z    between  dipoles  I  and  J  is  the  sum  of  four 
ij 

segment-segment  impedances. 

The  variables  IFLAG  and  JFLAG   are  used  if  a  ground 
plane  is  present  for  the  calculation  of  the  mutual  impedance 
elements.   If  IFLAG  is  equal  to  JFLAG,  the  mutual  impedance 
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terms  will  not  have  the  effects  of  a  ground  plane  since  both 
monopoles  lie  on  the  same  side  of  the  ground  interface.  If 
the  monopoles  are  on  the  opposite  sides  of  the  interface 
(IFLAG  not  equal  to  JFLAG) ,  the  reflection  coefficient 
correction  must  be  applied  to  the  mutual  impedance  elements. 
This  same  technique  is  applied  in  subroutines  GNFLD  and 
GFFLD. 

In  SGANT#  segment  K  has  endpoints  KA  and  KB,  and  segment 

L  has  endpoints  LA  and  LB.  It  is  convenient  to  think  of  KA 

and  KB  as  points  1  and  2  on  segment  K,  and  LA  and  LB  as 

points  1  and  2  on  L.   The  four  segment-segment  impedances 

can  be  defined  as  P(IS,JS).   The  first  subscript  IS  refers 

to  the  terminal  point  on  segment  K,  and  the  second  subscript 

JS  refers  to  the  terminal  point  on  L.  Thus  IS=1  or  2  if 

dipole  I  has  its  terminal  point  12(1)  at  KA  (point  1)  or  KB 

(point  2) ,  respectively.  Similarly,  JS=1  or  2  if  mode  J  has 

its  terminal  point  12  (J)  at  LA  or  LB.  The  impedances 

P(IS,JS)  are  defined  with  the  following  reference  directions 

for  current  flow:  from  point  1  toward  point  2  on  each 

segment.  If  dipole  I  has  this  same  reference  direction  on 

segment  K,  FI=1;  otherwise  FI=-1.  Similarly  FJ=1  or  -1  in 

accordance  with  the  reference  direction  for  mode  J  on 

segment  L.  In  statement  18,  P(IS,JS)  is  multiplied  by  FI  and 

FJ  before  its  contribution  is  added  to  Z 

ij 

Subroutine  GGMM  calculates  the  impedances  Q(KK,LL)  which 
are  like  the  P(IS,JS)  but  have  different  conventions  for 
reference  directions  and  subscript  meaning.  The 
transformation  from  the  Q  impedances  to  the  P  impedances  is 
accomplished  in  the  DO  LOOP  ending  with  statement  13. 

If  the  wire  has  finite  conductivity,  the  appropriate 
modification  is  applied  to  the  impedance  matrix  just  above 
statement  15.   The  terms  arising  from  the  dielectric  shell 
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on  an  insulated  segment  are  obtained  from  subroutine  DSHELL 
just  above  statement  16.  Finally,  the  lumped  loads,  ZLD,  are 
added  to  the  diagonal  elements  of  the  impedance  matrix  in 
the  DO  LOOP  ending  at  statement  23. 

K  is  a  segment  of  test  dipole  I,  and  L  is  a  segment  of 
expansion  mode  J.  When  the  segment  numbers  K  and  L  are 
egual,  SGANT  calls  GGMM  to  obtain  the  mutual  impedance 
between  two  filamentary  electric  monopoles.  These  monopoles 
are  parallel  and  have  the  same  length.  Monopole  K  is 
positioned  on  the  axis  of  the  wire  segment,  and  monopole  L 
is  on  the  surface  of  the  same  wire  segment.  Thus,  the 
displacement  is  equal  to  the  wire  radius.  The  two  monopoles 
are  side-by-side  with  no  stagger. 

When  segments  K  and  L  intersect,  SGANT  again  calls  GGMM 
for  the  mutual  impedance  between  the  two  filamentary 
monopoles.  Monopole  K  is  situated  on  the  axis  of  wire 
segment  K,  and  monopole  L  is  on  the  surface  of  wire  segment 
L.  The  axes  of  segments  K  and  L  define  a  plane  Prand 
monopole  K  lies  in  this  plane.  Monopole  L  is  parallel  with 
plane  P  and  is  displaced  from  it  by  a  distance  egual  to  the 
wire  radius. 

CALLED  BY:      MAIN 


CALLS  TO:       CBES 

DSHELL 

GGMM 

GGS 
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SUBROUTINE    SGANT     II  A  ,  I  B  ,  I  NM  ,  IN*T  ,  I  SC  .  I  1  ,  I  2  .  I  3,  JA  ,  JB  ,  MD.  N,  NO  .NM.NP  ,  *  I 

lM,BM,C,CGO,CMM,D,EP2,EP3.ETA,FHi,GAM,SGO,X,Y,2,il.O,2S,ERR,  IGRDI  2 

CCHPLEXERR  3 

COMPLEX    ZG^H^S.EGO.GD.CGOS.SGOS.SGDT.BOl  * 

CUMPLEX    Plf ,Pl2IP2l,P22,Qll,012,Q2l.Q22,EP2,EP.ETA,GAM,EP3  5 

COMPLEX     EPSiLA.CKEA.BEIA, ZARG  6 

COMPLEX    P(2,2I,0(2,2),CGOU I . SGDI 1 J ,CI I I.ZLOIll  7 

DIMENSION    XII).     Vlll,     HI).     Oil),     IAU),      IBID.     MOIINM,*)  8 

OIMENSION     lllll,     12(1).     nil),     JAI1),     18(11.    NDtlT,     ISCUI  9 

OATA    E0,TP,U0/8.854E-12,6.283l8,1.2566E-6/  10 

EP    ■    EP3  11 

ICC    -     (N«N»NI/2  12 

13 

00    1    l«It ICC  1* 

i  cm  ■   i.o,.o)  15 

16 

IS    *    1.0, .01  17 

IF    (CMM.LE.O. )    GO    TO    2  18 

OMEGA    =     TP«FH2  19 

EPS1LA    ■    CMPIXIE0.-CMMM.E6/0KEGA)  20 

CWEA    ■     I .0, I. >«OMEGA»EPS)lA  21 

BETA    •    OMEGA»SORTIU0)»CSQRMEPSILA-EPI  22 

2APG   -    BETA»AM  23 

CALL    CBFS    I2ARG.B01 )  2* 

IS    •    BfTA»BOl/CKEA  25 

2  IH  •  lS/{ TP«AM«GAMI  26 
OMIN  -  I.t30  27 
DMAX    -     .0  28 

?9 

00    3    J- I ,NM  30 

K    ■     IA( j|  31 

L    ■    1BI J)  32 

OIJ)    -    SQRTI(X<K>-X(U>**2MYIK>-YILM»*2»IZ(K)-2ILII**2)  33 

IF     IDIJI.LT.OMIN)     DMIN-DIJ)  34 

IF     101 JI.GT.DMAX)     DMAX-O(J)  35 

EGO    =    CEXPIGAM»01 Jl  )  36 

CGDIJI    »    IEG0M./EG0I/2.  37 

3  SGDUl    •    IEG0-1./EG0I/2.  38 

39 

IF     (0MIN.LT.2.«AMI     GO     TO    *  40 

IF     (CABS(GAM«AM> .GT. 0.061     GO    TO    4  41 

IF     (CAhSIGAM*0MAX).GT.3.)     GO    TO    4  42 

IF    IAM.GT.0.)     SO    TO    5  43 

4  CUNTINUE  44 
N-0  45 
WRITE  (6.241  AM.DMAX.DMIN  46 
WRITE    (6,25)  47 

48 


49 

50 
.AND. IK.GT.NN/2) I      IFLAG-l  51 

53 

54 
55 
56 
57 
58 
59 
60 
.AND.IL.GT.NI1/2I )    JFLAG-1  61 

6Z 
63 
64 
65 
66 
67 
68 

DO    19     11  =  1  .NDK  69 

I     •    HOIK,  III  70 

MM    -     II-[)»N-()»I-I 1/2  71 

Fl     "     1.  " 

IF    (KB. E0.I2I I) I    GO   TO   6  73 

IF     (KB.  EO. Ill  II)     Fl  — 1.  .  74 

IS    •    1  75 

GO    TO    7  76 

IF    IKA.E0.13I  1)1    Fl  — 1.  7T 

I S    ■    2  78 

79 
DO    19    JJ=l,NDL  80 

J    «    MOI L.JJI  81 

MMM    •     MMtJ  82 

IF    I  I  .GT.  J)    GO    TO    19  83 

FJ    •    I.  84 

IF     (LB.  EQ.I2U)  I    GO    TO    8  85 

IF    (LB.  EO.  II IJ)  )    FJ— I.  86 

JS    -    1  87 

GO    TO    9  88 

IF    (LA.E0.I3I  J)  )    FJ— 1.  89 

JS    •    2  ,0 

IF     (NIL.NE.OI     GO    TO    18  91 

NIL     ■     1  H 

IF     (K.EQ.LI     GO    TO     14  93 

IND    •     ILA-KAI«ILB-KAI») LA-KB I  * ( L B-KB)  94 

NGRO    •      IGRO  93 

IF     I1FLAG.EQ.JFLAGI     IGRO— I  96 


DO  19 

K-l  ,NM 

IFLAG 

*  0 

IF  (1 IGRD.GT.O) 

NDK  ■ 

ND(K) 

KA  - 

14  IK  1 

KB  ■ 

I3IK) 

OK  * 

Dl  K) 

CGOS 

-  CGO(K) 

SGOS 

•  SGOIKI 

DO  19 

I  •  1  .  NM 

JFLAG 

=  0 

IF  II IGRD.GT.O) 

NUL  « 

NOIL) 

LA  ' 

IAILI 

LB  - 

IBID 

OL  - 

D(  L) 

SCOT 

•  SGD(L) 

NIL  - 

0 
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c 


SF&MFNTS  K  AND  I  SHARE 

10 

KG 

-  0 

JM 

■  KB 

JC 

•  KA 

KF 

-  1 

ISO  =  ( *  B-LA )•( K8-LB) 

IF 

1  1N0.NE.0I  GO  TO  11 

JC 

■  KB 

KF 

■  -1 

JM 

=  KA 

KG 

'     3 

11 

LG 

-  3 

JP 

•  LA 

IF 

■  -I 

IF 

(LB.  EQ.  JC)  GO  TO  12 

JP 

"  LB 

LF 

'  1 

LG 

*    0 

IF     ( INO.EO.O)    GO    TO    10  97 

St&MENTS    K    ANO    L     SHARE    NO    POINTS  98 

CALL    GGS     (X(KA)  , Y ( KA I , 2 ( KA I .X ( KB  I ■ Y ( KB )  ,  2  I  KBI  ,  X(  I A  I  ,  Y<  L  A  I  .  2  I  L  A  I  ,Xl  99 

ILBI.KILBI  ,  2UBI  ,  AM . OK ,CGOS . S&OS . OL . SGUT ,  INT,ETA,GAN,P(l,ll,P(l,2l,  100 

2PI2, 1 ) ,P(2,2I .ERR, IGROI  101 

IGfcO    «     N&RO  102 

GO    TO    IB  03 

SHARE    ONE     POINT     (THEY     INTERSECT!  0* 

05 
06 
07 
08 
09 
10 
11 
12 
13 
I* 
15 
16 
1 1  7 
13 
19 
20 
21 

12  SGN    *    KF»LF  122 
CPS1     «     MX(JP)-X(  JCI  (•IX(JM1-X(  JC|)»(Y(  JP|-Y(JCIl»(Y(JM|-V(  JCIIMZ  123 

l(JP)-2(JCI)»(2<JM)-2(JCIlF/(DK*DLi  12* 

CALL    GGMM    <  .O.OK,  .O.DL.AM.CGOS.SGOS.SGOT.CPSI  .ETA.GAN.QI  1, 1  I  ,Q<  1,2  125 

II ,0(2,1 ),Q(2,2I )  126 

IF  <  K.PO.&I  .J)     SGN— SGN  127 

C  128 

DO    13    KK=1  ,2  129 

KP    •     IA8SIKK-KG*  30 

C  31 

DO    13    LLM.2  32 

LP    '     UBSILL-LGI  33 

PIKP.LP)    •    SGN*0(KK,LLI  3* 

13  CONTINUE  35 
C  36 

IGRD-NGRO  ,  37 

GO    TO    18  '  58 

C  K-l       (SELF    REACTION    OF    SEGMENT    Kl  139 

I*    01 1     ■(  .0,  .01  *0 

Q12    -    I .O..0I  1*1 

IF    ICMM.LE.O.I    GO    TO    15  1*2 

GO   ■    GAM»DK  43 

2G   -    2H/(SG0S»»2»  ** 


Oil    ■    2&«ISGDS»CGOS-GDI/2.  \  *J 

012    ■    2&»IGD»CGDS-SGDS)/2.  }** 

15  ISCK  -  ISCIKI  1*7 
Pll  '  I .0..0I  1*8 
P12  '  (  .0..01  J* 
IF  ( I SCK.EO.OI  GO  TO  16  fO 
IF  IBM. LE. AMI  GO  TO  16  51 
CALL    DSHELl     IAM,BN,OK.CGOS.SGDS,EP2.EP,ETA,GAM,Pll,Pl2l  52 

16  Qll  -  Pll»QU  ll 
CALL"GGMM*(.0,3K,.0(DK,»H,CGDStSGOS,SGDSil.  tETA,GAH,PU,P12fP2ltP2 

12 1  '■ 

Oil    -    PI  l»Qt 1  57 

012   «    P12»012  58 

Plltl  I    -    Oil  *n 

P(lt2l     «    012  5? 

PI2.1I     •    012  *1 

P(2,2I     =    Oil  'J? 

IF     IKA.NE.LA)     GO     TO    17  bi 

GO    TO    18  ** 

17  Pll .11     -    -012  *5 

pi  1.2)  ■  -on  55 

p(2,ti  *  -on  5; 

PI2.2I     *    -012  5* 

18  CIMMMI  =  CIMMMI »FI»FJ*P( IS. JSI  69 

19  CONTINUE  J? 

DO    23    l«l,N  H 

MM    -     (  l-ll»N-(  1*1  —  1  1/2  J* 


1  J    •    MMH 


75 


JJA    ■     JAI II  1% 
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Jl     ■     JJA 

112    -    121 II 

111    •     11(11 

IF     (  I  I,'. I  U.IBI  Jl  l»  J1-J1»NM                                                                                                                                       "9 

JJB    •     JBI  I  I  'I 

J2    »    JJB  l\ 

IF    ( II2.E0.IB1 J2II  J2-J2»NM                                                                                                                             83 

CIIJI    •    C(  IJI»2LOUll»ilOU2l  %i 

JJJ  •  JJA  J5 

C  86 

00  22  K-1,2  J' 

NOJ  •  NDIJJJI  88 

00  21  JJ-l.NOJ  2? 

J  •  M0( JJJ. JJ I  'I 

IF  (  J.fcO.l  I  GO  TO  21  "2 
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IF    IIJUI.Nt.IUI    GO    TO   21  193 

I  J   -    MN*J  19* 

Fl    =    I.  195 

IF    IK. £0.21    GO    TO    20  196 

IF     (111  JI.NE.lll)     Fl  — 1.  197 

CIIJ)     -    CUJI»FI»ZLDIJ1I  198 

GO    TO    21  199 

20  IF    II3UI.NE.I3II  II    Fl  — 1.  200 
CIIJI    *    CUJ)*FI»2LDIJ2)  201 

21  CONTINUE  202 

203 

22  JJJ   ■    JJB  20* 

205 

23  CONTINUE  206 

207 
RETURN  208 

209 
2*    FORMAT     C3X,'*N   •    ■  ,  E 1 0. 3 , 3X , • DMAX    ■    •  ,E  1  0.  3 ,  3X  ,  •  DM  IN   -    '.EIO.SI  210 

25    FORMAT     (•     WARNING    *•••••••♦••••«••••»*•••••♦♦••••»••*•«•*•••»••»•/  211 

1,'    THIS    PROBLEM    EXCEED   LIMIT    OF    THIN   MIRE    CONDITION,    THE    RESULTS  212 

2    ARE    NOT   CORRECT' I  213 

ENO  21* 
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SORT 


PURPOSE:       to  define  the  set  of  dipole  modes. 

METHOD:        In  the  DO  LOOP  ending  at  statement  3,  the  set 
of  dipoles  is  defined  by  filling  the  vectors  11(1)  and  13(1) 
(the  endpoints  of  dipole  I)  ;  12  (I)  (  the  terminal  point  of 
dipole  I);  and  the  vectors  JA  (I)  and  JB(I)  (the  monopoles 
comprising  dipole  I)  with  the  node  numbers  and  segment 
numbers,  respectively.  The  DO  LOOP  ending  at  statement  8 
determines  MD(J,K)  (the  list  of  dipoles  sharing  segment  J) 
and  ND  (K)  (the  number  of  dipoles  sharing  segment  J). 

CALLED  BY:      MAIN 
CALLS  TO:       NONE 
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SORT  tlA,IB,ll,I2,l3iJA,JB,ND,ND,NM,NP,N.MAX,MIN,ICJ,IN 

.  13(1),  JA<1) 
.  NOUI.  MDIIh 


SUBROUTINE 
IMI 
DIMENSION  JSPI20) 
01  MENS  I  ON  II  III,  1211 
OIMENSION  14(1),  18(1 
I  -  0 


00  3  K=1,NP 
NJK  «  0 

DO  1  J-l.NM 

INO  -  I  1AI Jl-K)*(  I8IJI-K) 

IF  I INO.NE.OI  GO  TO  1 

NJK  =  NJK»1 

JSP(NJK)  ■  J 

CONTINUE 

MOO  =  NJK-1 

IE  (MOD.LE.OI  GO  TO  3 

00  2  IMD-l.MOO 

1  -  I»l 

IF  ll.GT.ICJI  GO  TO  2 

I  PD  «  IMDU 

JAI  •  JSP( IMOI 

JAI I  I  ■  JAI 

JBI  ■  JSPIIPOI 

JBII1  '  JBI 

11(11  -  IA1JAI) 


I),  JBI1) 

INM.il 


c 
c 

c 

c 

2 

IE     1IA(  JAD.EO.KI 

12(1)     ■     K 

13(1)     •     lA(JBI) 

IF     ( IAI JBI I.EO.KI 

CONTINUE 

IK  D-IBI  JAI) 
1 M I l-IBI JBI) 

3 

CONTINUE 

N    •    1 

00    4    J-l.NM 
NOUI    •    6 

c 

c 

4 

00    4    K-1,4 

MOIJ.K)     •     0 

111    •    N 

IF    (N.GT.ICJ1     II  1 

-    ICJ 

1 

2 

3 

5 
6 
7 
B 
9 
10 
11 
12 
13 
1* 
15 
16 
IT 
18 
19 
20 
21 
22 
23 
2* 
25 
26 
27 
28 
29 
30 
31 
32 
33 
3* 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 


DO  8  1-1,111 
J  -  JAI II 

00  7  L-1,2 
NOUI  •  ND(J)»l 
K  -  1 

M  -  0 

5  MJK  -  MOIJ.K) 

1  E  (MJK  .NE.O)  GO  -TO  6 
M  «  1 

MOIJ.K)  •  I 

6  K  ■  K»l 
IF  1K.GT.4) 
IF  IM.E8.0l 

7  J  -  JBI  I) 

8  CONTINUE 


GO  TO  7 
GO  TO  5 


MIN 

MAX 


100 
0 


00  9  J-l.NM 
NOJ  •  NOlJI 
IF  (NOJ.GT.MAXI 
IF  (NDJ.LT.MlNl 

RETURN 
ENO 


M4X-NDJ 
M1N-NDJ 


49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 

i? 

73 
74 
75 
76 
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SQROT 


PURPOSE:       to  solve  the  set  of  simultaneous  equations  to 
determine  the  currents  on  the  thin  wire  structure. 

METHOD:        This  subroutine  considers  the  matrix  equation 
ZI  =  V  which  represents  a  system  of  simultaneous  linear 
equations.  NEQ  denotes  the  number  of  simultaneous  equations 
and  the  size  of  the  matrix  Z. 

On  entry  to  SQROT,  S  is  the  excitation  column  V.  On 
exit,  the  solution  I  is  stored  in  S.   Z(I,J)  denotes  the 
symmetric  square  matrix.   Also  on  entry,  the  upper-riqht 
trianqular  position  of  Z(I,J)  is  stored  by  rows  in  C(K)  with 

K  =  (I  -  1)  *  NEQ  -  (I  *  I)  /  2  +  J  . 


If  112  =  1,  SQROT  will  transform  the  symmetric  matrix 
into  the  auxiliary  matrix  (implicit  inverse) ,  store  the 
result  in  C  (K)  and  use  the  auxiliary  matrix  to  solve  the 
simultaneous  equations.  If  112  =  2,  this  indicates  that  C(K) 
already  contains  the  auxiliary  matrix. 

The  transfromation  from  the  symmetric  matrix  to  the 
auxiliary  matrix  is  accomplished  in  the  DO  LOOP  endinq  at 
statement  5.   The  solution  of  the  simultaneous  equations  is 
accomplished  in  the  remainder  of  the  proqram. 

CALLED  BY:      GFFLD 


CALLS  TO:       NONE 
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SUBROUTINE    SQROT    ( C . S 1 1 WR , I  12 iNEQI  I 

COMPLEX    C(  II ,S(1 I tS4  2 

N    -    NEO  '3 

IF    (  112. EQ. 21    GO    TO   6  4 

CI  1  1    -    CSQRUCI  III  5 

C  6 

DO    I    K=2,N  7 

1  CIKI  '  C(KI/C( 1  I  8 
C  9 
C  10 

DO    5    l»2,N  11 

I  MO    -     I -I  12 
IPO   •    I»l  13 

10  •    I  I -1 >  *N-I 1*1-1 1/2  1* 

II  •    10*1  15 
C  16 

00    2    1-1,1  HO  [7 

11  -    (l-il»N-(L«L-l)/2M  18 

2  CUM  ■  Cllll-CILI  l-CILI)  19 
C  20 

Cf  1 1  >    ■    CSORMCIll  II  21 

IF    I IPO.CT.NI    G0    TO   5  22 

C  23 

00   4    J-  IPO.N  24 

IJ    ■    ID»J  25 

C  26 

00  3  M. i. | MO  27 

MD  »  IM-I I-N-1M-M-MI/2  28 

HI  -  MO»l  29 

MJ  -  HO»J  30 

3  CIIJI  •  CI  IJ1-CIHJ)*C(HI  I  31 
C 

4  CIIJI  •  CIIJI /C (II) 
C  3* 

5  CONTINUE  35 
C  36 

6  Sill  •  SUI/CU)  37 
C  38 

00  8  1-2, N  39 

IMO  •  l-i  40 

C  41 

42 

/2»I  43 

..-I  44 

C  45 

1  I  -  (!-l)»N-U»l-!l/2U  46 
6  Sll)  •  SI  II/CI1I)  47 

C  48 


C 


00  7  L»l. I  HO 

LI  -  U-Il»N-(L«L-U/ 

Sll)  -  Si  I l-CILI )»S(l 


3i 


NN  -  (  <N»1 l»N)/2  J' 

SIN)  ■  SINI/CINN)  £° 

NMO  «  N-l  I' 

C  ' 

00  10  1-1, NMO  ll 

K  =  N-l  54 

KPO  =  K*l  ll 

KD  -  (K-l )»N-(K»K-KI/2  56 

C 

DO  9  L-KPO.N  58 

9  SIKI  -  S(K)-C1KL)»S(L)  *? 
C 

KK  -  KO»K  *| 

10  SIKI  -  S(KI/C(KK)  *£ 

c 

IF  (  IWR.LE.O)  GO  TO  13  J5 

CNOR  -  .0  |° 

00  11  I «1.N  Jg 

SA  -  CABSlSI II)  ,  *' 

11  IF  IS*. GT. CNOR)  CNOR-SA  J° 

IF  (CNOR.LE.O.)  CNOR-1.  J,| 

00  12  I-l.N  7* 

SS  -  SI  II  I? 

SA  •  CABSISS)  IS 

SNOR  -  SA/CNOR  >'  JJ 


?F  "sa?GT.0.1  PH  ■  57.29578«ATAN2IAIHAGISS),RE»USS>»  79 

12  WRITE  16.14)  1 ,SNOR,SA,PH,SS  |° 

KR1TE  (6,151  *{ 

13  RETURN  "* 

14  FORHAT  (IX,1I5,IFI0.3,IFH.7,1FIO.O,2F15.6)  JJ 

15  FORHAT  I1H0I  St 
END  ■* 


V 


9 . «BLNK 

..?•*•••  'EQUAL 
.READ -If:.. 


LEFT 
* 'NUMBER 


• 


SORT 


A 
I 

N 


.  .CBES 
DSHELL 


SGAHT.  U  J...     ^ 

GANT1 © 

*  GNFLD GNF 


*  * • .  .  .  •  • GFF 

GFFLD.  •  t I  *  ^ 


....  LINECK 

.  «polplot»  • :; . 

•  •  •  •   WFTMR 

*  "  '  »SART 


NUMB 
POLPRT*  : : 


®.  .  •  »EXPJ 
*  •  •  .  . 
•  •  «GGMM 

.  • GDISS 

(§)•  |  •  •  •  •  *RITE 

•  •  • 

•SQROT 

CALLING  SEQUENCE  OF  THE  SUBROUTINES 
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SYMBOL  DICTIONARY 


A         characters  of  the  input  data  cards 

ABAP      backscattering  phi  plane  angle  for  plotting 

ABAT      backscattering  theta  plane  angle  for  plotting 

ABIP      bistatic  scattering  phi  plane  angle  for 
plotting 

ABIT      bistatic  scattering  theta  plane  angle  for 
plotting 

ACSP      absorption  cross  section  for  phi  polarization 

ACST      absorption  cross  section  for  theta 
polarization 

AFFP      far-zone  phi  plane  angle  for  plotting 

AFFT      far-zone  theta  plane  angle  for  plotting 

AM        radius  of  the  thin  wire  of  the  structure 

BM        outer  radius  of  the  dielectric  shell  of  the 
insulation  of  the  wire 

C  elements  of  the  open-circuit  impedance  matrix 

CG  branch  currents  for  the  structure 

CGD  cosh  yd   for  a  given  segment 

CJ  loop  currents  for  the  structure 

CMM  conductivity  of  the  wire 

D  length  of  a  given  segment 

ECSP  extinction  cross  section  for  phi  polarization 

ECST      extinction  cross  section  for  theta 
polarization 

EFF       radiation  efficiency 

EP        loop  currents  induced  by  a  phi  polarized  wave 

EPP       phi-polarized  far-zone  feild  of  the  dipole 
mode 

EPPS      scattered  electric  field  in  the  phi  direction 
due  to  a  phi  polarized  wave 

EPTS      scattered  electric  field  in  the  theta 
direction  due  to  a  phi  polarized  wave 

EP2  complex  permittivity  of  insulation 

EP3  complex  permittivity  of  ambient  medium 

EP4  complex  permittivity  of  ground 

ERF  EP4/EP3 

ER2  relative  dielectric  constant  of  insulation 
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ER3       relative  dielectric  constant  of  the  ambient 
medium 

ER4       relative  dielectric  constant  of  the  ground 

ET        loop  current  induced  by  a  theta  polarized 
wave 

ETA       intrinic  impedance  of  ambient  medium 

ETPS      scattered  electric  field  in  the  phi  direction 
due  to  a  theta  polarized  wave 

ETT       theta  polarized  far-zone  field  of  the  dipole 
mode 

ETTS      scattered  electric  field  on  the  theta 
direction  due  to  theta  polarized  wave 

EX  near-zone  electric  field  in  x  direction 

EY  near-zone  electric  field  in  the  y  direction 

EZ  near-zone  electric  field  in  z  direction 

EO  8.854E-12 

FHZ  frequency  in  hertz 

FMC  frequency  in  megahertz 

GAM       intrinic  progration  constant  of  the  ambient 
medium 

GG        time-average  power  input 

GPP       power  gain  associated  with  the  phi  polarized 
component 

GTT       power  gain  associated  with  the  theta 
polarized  component 

height  of  the  structure  above  ground  plane 

first  node  of  a  given  segment 

second  node  of  a  given  segment 

indicator  for  bistatic  scatter  calculations 

indicator  for  the  data  cards 

dimension  corresponding  to  the  number  of 
simultaneous  linear  equations 

indicator  for  proqram  termination 

indicator  for  antenna  qain  calculations 

indicator  for  presence  of  the  ground  plane 

indicator  for  the  type  of  far-zone 
calculations 

indicator  for  near-zone  calculations 

dimension  corresponding  to  the  number  of 
mon poles 

INT       number  of  integration  steps 
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HGT 

IA 

IB 

IBISC 

ICA8D 

ICJ 

IFLAG 

IGAIN 

IGRD 

INC 

INEAR 

INM 

ISC 

ISCAT 

IWR 

11 

112 

12 

13 

Jl 

JB 

KFLAG 

KGEN 

LOAD 

LZD 

MAX 

MD 

MIN 

MSG 

N 

ND 

NGEN 

NM 

NPL 

OMEGA 

PH 

SCSP 

SCST 

SGD 

SIG2 

SIG3 

SIG4 

SPPM 

SPTM 

STPM 

STTM 

indicator  for  the  insualtion 

indicator  for  backscatter  calculations 

indicator  for  current  distribution  output 

endpoint  node  of  a  given  dipole 

indicator  for  auxiliary  matrix 

terminal  node  number  of  a  given  dipole 

endpoint  node  number  of  a  given  dipole 

first  segment  number  of  a  given  dipoile 

second  segment  number  of  a  given  dipole 

print  indicator 

list  of  generator/feed  locations 

indicator  for  structure  load 

list  of  impedance/load  locations 

maximum  of  the  number  of  segments  connected 
to  any  one  given  node 

list  of  dipoles  sharing  a  given  segment 

minimum  of  the  number  of  segments  that 
connected  to  any  one  given  node 

indicator  for  error  printout 

number  of  simultaneous  linear  eguations 

total  number  of  dipoles  sharing  a  given 
segment 

indicator  for  antenna  calculations 

number  of  segments 

indicator  for  polar  plot 

angular  frequency 

phi  angle  for  far-zone  calculations 

scattering  cross  section  for  phi  polarization 

scattering  cross  section  for  theta 
polarization 

sinh  yd   of  a  given  segment 

conductivity  of  insulation 

conductivity  of  the  ambient  medium 

conductivity  of  ground 

echo  area  phi  incident-phi  scattered  wave 

echo  area  phi  incident-theta  scattered  wave 

echo  area  theta  incident-phi  scattered  wave 

echo  area  theta  incident-theta  scattered  wave 
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TD2  loss  tangent  of  the  insulation 

TD3  loss  tangent  of  ambient  medium 

TH  theta  angle  for  far-zone  calculations 

TP  2n  (6.28318) 

00  1.2566E-6 

VG  antenna  complex  driving  voltages 

VOLT  list  of  VG«s 

X  x-coordinate  of  each  node 

XNP  list  of  XP's 

XP  x-coordinate  for  near-zone  calculations 

y  y-coordinate  of  each  node 

YNP  list  of  YP's 

YP  y-coordinate  for  near-zone  calculations 

Y11  complex  power  input 

Z  z-coordinate  of  each  node 

ZLD  complex  load  at  a  given  node 

ZLLD  list  of  ZLD's 

ZNP  list  of  ZP's 

ZP  z-coordinate  for  near-zone  calculations 

ZS  surface  impedance  of  the  wire 

Z 1 1  antenna  input  impedance 
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APPENDIX  B 
USEE'S  MANUAL 
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USER'S  MANUAL 

The  Antennas-Scatterers  Analysis  Program  (ASAP)  for  thin 
wire  structures  in  a  homogenous  conducting  medium  performs  a 
freguency  domain  analysis  of  antennas  and  scatters.  The 
program  is  applicable  in  the  presence  of  either  a  perfect  or 
a  finite  ground.   This  appendix  will  describe  and  explain 
the  data  cards  necessary  to  execute  the  compute  program. 
Although  the  program  was  written  for  the  IBM  360  computer 
system,  it  can  be  executed  on  another  system  with  minor 
modifications. 

The  program  utilizes  piecewise  sinusoidal  expansion  for 
the  current  distribution  with  Kirchhoff  Current  Law  enforced 
everywhere  on  the  structure.  If  the  structure  contains  end 
points,  the  currents  at  these  points  are  assumed  to  vanish. 

I.  Program  Limits 

The  thin  wire  assumptions  are  guestionable  and  the 
accuracy  and  convergence  deteriorate  if  the  radius  of  wire 
utilized  for  the  structure  exceeds  0.01   of  a  wavelength,  if 
the  longest  segment  is  greater  than  one-fourth  of  a 
wavelength,  if  the  length  ratio  of  the  longest  and  shortest 
segments  exceeds  100,  or  if  the  total  wire  length  is  less 
than  30  times  the  wire  diameter.   If  a  wire  is  bent  sharply 
to  form  a  small  acute  angle  (less  than  30  degrees) ,  the  thin 
wire  model  is  guestionable.  It  is  assumed  that  the  wire 
conductivity  greatly  exceeds  the  conductivity  of  the  ambient 
medium.  For  insulated  wires,  the  dielectric  layer  is  assumed 
to  be  electrically  thin. 

II.  Minimum  Data 

The  minimum  data  necessary  to  execute  the  program  is: 
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a.  description  of  structure 

b.  radius  of  wire  used  for  the  structure 

The  program  will  default  to  the  other  parameters  necessary. 
The  default  parameters  are: 

a.  wire  for  the  structure  is  copper 

b.  frequency  of  operation  is  300  mhz 

c.  homogneous  medium  is  free  space 

A  more  detailed  explanation  of  the  defaults  will  be 
discussed  when  the  data  card  for  the  parameter  is  described. 

III.  Outputs 

In  antenna  problems,  the  output  includes  structure 
currents,  impedance(s)  of  f eed (s) ,  gain,  polar  radiation 
plots,  and  near  field  calculations.  In  bistatic  scattering 
problems,  the  output  includes  structure  currents,  complex 
elements  of  the  polarization  scattering  matrix,  polar 
reradiation  pattern  plots,  and  echo  areas  produced  by  a 
plane  wave.   For  backscattering  problems  the  output  includes 
absorption,  scattering  and  extinction  cross  sections  in 
addition  to  the  outputs  of  bistatic  scattering.  Most  of  the 
outputs  are  suppressed  and  must  be  requested.  Since  the 
program  can  produce  a  large  volume  of  output,  care  should  be 
exercised  until  the  user  is  familiar  with  the  outputs. 

IV.  Data  Cards 

The  Analysis  Program  utilizes  free  format  for  the  data 
cards,  that  is,  the  program  utilizes  character  recognition 
to  determine  which  parameters  are  being  read.  Data  placement 
(location)  on  the  input  card  is  not  critial.  Blank 
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characters,  on  all  input  cards  but  the  COMMENT  data  card, 
are  ignored  and  may  be  used  at  the  discretion  of  the  user. 
Since  character  recognition  is  used,  only  the  first  four 
characters  of  the  key  words  must  be  present  and  correct. 

The  format  for  the   COMMENT  CABD  utilizes  standard 
FORTRAN  format  (i.e.  'C  in  column  1  followed  by  at  least 
four  blanks) .   The  COMMENT  CARD  is  the  only  type  of  input 
card  that  position  in  the  data  block  is  critial.   This 
(these)  card  (s)  must  be  placed  at  the  beginning  of  a  data 
block.  A  data  block  is  a  series  of  related  data  cards. 
Several  data  blocks  may  be  used  to  define  an  analysis 
problem.  This  will  become  clear  when  the  termination  cards 
(END,  STOP,  or  CHANGE)  are  discussed.   There  is  no  limit  to 
the  number  of  comment  cards  that  may  be  used.  As  a  check  for 
the  user,  all  input  data  cards  will  appear  on  the  output  as 
they  appear  in  the  input  deck. 

The  format  of  other  data  can  be  of  one  of  two  forms: 

a.  type  of  card  (option  1/option  2/ ) 

b.  parameter  (value)  . 

The  type  of  format  to  use  will  be  apparent  as  the  individual 
data  cards  are  discussed. 

The  numerical  values  for  the  parameters  may  be  stated  in 
any  one  of  the  following  forms.  The  program  will  translate 
the  number  to  the  proper  form  for  the  specified  parameter, 
either  fixed  or  floating  point.   All  of  the  following 
examples  have  the  same  value. 

0.0001  or  .0001  or  100. U  or  100U  or  .1M  or  0.1M  or  .0000001K 

-6  -3  3 

U  =  10  M  =  10  K  =  10 
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1-  .WIItE  This  card  is  used  to  define  the  parameters 
associated  with  the  wire  utilized  by  the  thin  wire 
structure.  Two  options  are  available  and  are  defined  as: 

RADIOS=value  of  the  radius  of  the  wire  in  meters 

CONDUCTIVTTY=value  in  megamhos  per  meters  . 

The  wire  data  card  must  appear  in  the  first  data  block  to 
define  wire  radius.   The  default  value  of  the  conductivity 
is  50  megamhos/meter  (copper) . 


IIRE  (   RADIUS=.001/   CONDU=28.5) 
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2.  INSOLATION  This  card  is  utilized  to  define  the  parameters 
associated  with  the  insulation  of  the  wire  used  for  the 
structure  to  be  analyzed.  If  this  card  is  omitted,  the 
program  assumes  that  the  structure  is  uninsulated.  Four 
options  are  available  and  are  defined  as: 

RADIOS=value  of  outer  radius  in  meters 

CONDOCTIVITY=value  in  micromhos  per  meter 

DIELECTRIC=value  of  relative  dielectric  constant 

LOSS  TANGENT=value  . 

The  conductivity  and  either  the  relative  dielectric  constant 
or  the  loss  tangent  (but  not  all  three)  options  may  be 
stated. 


Insolation  (  radios=.015/  cond=7./diel=5) 
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3-  EXTERIOR  MEDIUM  This  card  is  utilized  to  describe  the 
homogeneous  medium  surrounding  the  structure.  If  the  medium 
is  free  space,  this  card  may  be  omitted.   Three  options  are 
available  and  are  defined  as: 

DIELECTRIC=value  of  relative  dielectric  constant 
CONDUCTIVITY=value  in  micromhos  per  meter 
LOSS.  TANG2NT=value  . 

As  with  INSULATION   card  state  either  conductivity  or  loss 
tangent. 


/EXTE(LOSS=.45) 


106 


4.  DESCRIPTION  This  card  is  utilized  to  describe  the  shape 
of  the  wire  structure  to  the  program.  The  user  must  divide 
the  wire  structure  into  segments  of  the  appropriate  length 
and  number  each  node  starting  at  one.  A  node  is  a  point 
where  a  segment  begins  or  ends.  A  maximium  of  four  segments 
can  meet  at  any  given  node.   An  isolated  wire  must  contain 
at  least  two  segments  and  three  nodes.   The  structure  is 
described  by  stating  the  node  numbers  that  each  segment 
connects.  The  description  of  a  square  loop  might  appear  as: 

DESCRIPTION (1-2/2- 3/3-4/4-1)  • 

The  description  of  a  dipole  and  reflector  might  appear  as: 

DESCRIPTION (1-2/2-3/3-4/4-5/6-7/7-8/8-9/9-10) . 

If  the  description  will  not  fit  on  one  data  card  continue  on 
the  next  card  as  if  the  previous  card  were  longer.  The 
dipole  example  might  appear  as: 

DESCRIPTION ( 1-2/2-3/3-4/4-5/ 

6-7/7-8/8-9/9-10) . 

Note  that  the  last  character  on  the  card  to  be  continued  is 
a  slant  (/) .  As  many  cards  as  necessary  may  be  used.  The 
maximum  number  of  nodes  permitted  is  fifty.  If  ground  plane 
is  present,  the  maximum  number  is  twenty-five.  If  a  ground 
plane  is  present  and  the  structure  touches  the  ground  plane, 
the  lowest  node  numbers  MUST  be  used  for  the  touching  nodes. 
That  is,  if  the  structure  touches  the  ground  plane  at  two 
points,  node  numbers  1  and  2  MUST  be  assigned  to  these 
nodes. 


'DESCRf  1-2/2-3/3-4/4-1) 
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5.  GEOMETRY  This  card  is  used  to  state  the  physical  location 
in  rectangular  coordinates  of  each  node  of  the   DESCRIPTION 
CARD  .  The  rectangular  grid  is  in  units  of  meters.  If  node  1 
is  located  at  x1,y1,z1  and  node  2  at  x2,y2,z2  and  node  3  at 
x3,y3,z3#etc. ,  the   GEOMETRY  CARD   might  appear  as: 

GEOMETRY  (x1,  y1, z1/x2, y2, z2/x3, y 3, z3/ ) 

As  with  the   DESCRIPTION  CARD,  continuation  cards  are 
permitted. 


'GEO  M(.  1,0,.  1/-.  1,0.1/-.  1,0-.1/.1,0,-.1) 
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6.  FEED  For  antenna  analysis  the  feed  point  (s)  and 
voltage  (s)  must  be  stated.  In  the  foremention  dipole  and 
reflector  example  if  the  feeds  were  at  node  2  with  a  voltage 
source  of  .5  at  an  angle  of  -90  degrees  and  at  node  4  with  a 
voltage  source  of  .5  at  an  angle  of  +90  degrees  the   FEED 
CARD   might  appear  as: 

FEED (2 ,.5,-9 0/4,. 5 ,+9  0) 

The  order  of  the  information  for  each  voltage  source  is  node 
number,  magnitude,  and  phase  angle.  This  order  is  repeated 
until  all  sources  are  stated.  If  the  source  information  will 
not  fit  on  one  card,  use  another  card  similiar  to  the 
initial  one;  that  is,  repeat  the  word  "FEED11.   If  only  one 
voltage  source  is  applied  to  the  structure,  only  the  node 
number  must  be  stated.  In  the  dipole  example,  if  the  drive 
is  at  node  3,  the   FEED  CARD   might  appear  as: 

FEED(3) 

A  default  source  of  one  voltage  at  zero  degree  phase  is 
assumed.  Voltage  sources  should  only  be  stated  for  nodes 
with  only  two  segments. 


TEED  (2,.  5, -9  0/4,  .5, +9  0) 
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7.  LOAD  This  card  is  used  to  describe  the  loads  to  be  placed 
at  various  locations  on  the  structure.  The  format  for  this 
card  is  similiar  to  that  of  the   FEED  CARD,  that  is,  the 
word  "LOAD"  is  used  in  the  place  of  "FEED".   The  order  of 
the  information  on  the  card  is  the  same.   Since  this  card  is 
frequency  dependent,  it  must  be  changed  if  the  frequency  of 
operation  is  changed.  No  default  parameters  are  available. 
The  structure  is  assumed  unloaded  unless  this  card  is  used. 
Once  the  structure  is  loaded,  it  will  remain  loaded  for  the 
remainder  of  the  data  block  series.  To  unload  the  structure 
the  following  card  may  be  used: 

L0AD(-1) 


'LOAD  (1,1 20, -4 5/3, 120,  +  4 5) 
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8.  OUTPUT  This  card  is  used  to  request  output  data.  Most  of 
the  output  is  in  tabular  form.   More  than  one   OUTPUT  CARD 
is  permitted  per  data  block,  but  not  for  the  same  type  of 
output.   If  only  the  antenna  input  impedance,  antenna 
efficiency,  or  time-average  power  input  is  of  interest,  no 
OUTPUT  CARD  is  necessary.  These  parameters  are  automatically 
printed  if  a   FEED  CARD   or  GENERATOR  CARD   is  utilized.  One 
or  more  of  the  following  options  may  be  used  to  request  the 
various  outputs  available. 

FAR  FIELD=phi  initial,  phi  final,  theta  initial,  theta 
final 

This  option  gives  the  components  of  the  electric  field 
intensity  in  the  far  field  as  phi  and  theta  varies  between 
limits  specified  in  one  degree  divisions. 

BACKSCATERING=phi  initial,  phi  final,  theta  initial, 
theta  final 

This  option  gives  the  absorption,  scattering,  and  extinction 
cross  sections,  and  the  complex  elements  of  the  polarization 
scattering  matrix  for  an  incident  plane  wave  illuminating 
the  structure  from  the  spherical  direction  of  phi,  theta  as 
both  vary  between  limits  specified  in  one  degree  divisions. 

BISTATIC=phi  inital,  phi  final,  theta  initial,  theta 
final 

This  option  gives  echo  area  and  the  complex  elements  of  the 
polarization  scattering  matrix  for  an  incident  plane  wave 
illuminating  the  structure  from  the  spherical  direction  phi, 
theta  final  of  the  backseat tering  output  option,  reradiated 
in  the  phi,  theta  direction  as  both  vary  between  limits 
specified  in  one  degree  divisions.  A  bistatic  output  request 
must  be  accompanied  with  a  backscattering  request  in  the 
same  data  block. 

STEP=value  in  degrees 
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This  option  will  cause  any  of  the  above  output  options  to  be 
stepped  at  a  different  interval  size.  That  is,  if  one  of  the 
above  options  is  to  be  stepped  at  ten  degrees  intervals,  use 
this  option.  This  option  overrides  the  one  degree  stepping. 

CURRENT 

This  option  gives  the  currents  on  the  structure  which  are 
produced  by  the  feed/generator  voltages  and/or  the  incident 
plane  wave  of  the  backscattering  request. 

NEAR=x1,y1,z1 

or 


NEAR=(x1,y1,z1/x2,y2,z2/x3,y3,z3/etc. 


) 


This  option  gives  the  value  of  electric  field  components  in 
the  near  field  for  the  antenna  at  the  point  or  points 
specified. 


/OUTPUT  (FARF=45,50,25,50) 


112 


9-  PLOT  This  card  will  produce  normalized  polar  plots  in  the 
specifed  plane  for  the  stated  option.   The  plane  is 

specified  by  stating  either  "PHI= "  or  "THETA= '• .  The 

PLOT  CARD  overrides  the  limits  of  the   OUTPUT  CARD  for  the 
same  option.  If  only  a  normalized  pattern  is  of  interest, 
only  a   PLOT  CARD  is  necessary.  If  a  table  of  values  and  a 
normalized  pattern  is  desired,  both  a   PLOT  CARD  and  OUTPUT 
CARD  must  be  used.   Only  one   PLOT  CARD   is  permitted  per 
data  block.   The  following  pattern  plots  are  available: 

FAR  FIELD/plane 

This  option  will  plot  the  far  field  intensity  for  each 
component  of  the  electric  field. 

BACKSCATTERING/plane 

This  option  will  plot  the  normalized  magnitude  of  each  of 
the  elements  of  the  polarization  scattering  matrix. 

BISTATIC/plane 

This  option  will  plot  the  normalized  magnitude  of  each  of 
the  elements  of  the  polarization  scattering  matrix  produced 
by  the  incident  plane  wave  stated  by  final  limits  of  the 
backscattering  option  of  the  output  reguest. 


'PLOT  (FARF/THET=90) 
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10.  GROUND  This  card  is  used  to  describe  the  ground 
parameters  if  a  ground  plane  is  present.  If  no  ground  plane 
is  present,  the  structure  is  assumed  to  be  in  free  space  or 
the  homogeneous  medium  of  the   EXTERIOR  MEDIUM  data  card. 
Seven  options  are  available  and  are  defined  as: 

PERFECT 

This  option    will   analyze  the   structure   over   a   perfect    ground 
plane. 

GOOD 

This  option  will  analyze  the  structure  over  a  good  ground 
plane  where  the  conductivity  of  the  ground  is  .02  mhos/meter 
and  the  relative  dielectric  constant  is  30. 

POOR 

This  option  will  analyze  the  structure  over  a  poor  ground 
plane  where  the  conductivity  of  the  ground  is  .001 
mhos/meter  and  the  relative  dielectric  constant  is  4. 

SEA 

This  option  will  analyze  the  structure  over  salt  water  where 
the  conductivty  of  the  water  is  4.  mhos/meter  and  the 
relative  dielectric  constant  is  80. 

HEIGHT=value  in  meters 

This  option  will  analyze  the  structure  with  origin  of  the 
GEOMETRY   card  this  height  above  the  ground  plane.  The 
lowest  point  of  the  structure  must  not  lie  below  the  ground 
plane.  It  may  lie  on  the  ground  plane. 

CONDUCTIVITY=  value  in  mhos/meter 

This  option  is  used  to  state  the  value  of  conductivity  of 
the  ground  plane  if  the  default  values  mentioned  above  are 
not  utilized. 
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DIELECTBIC=  value 

This  option  is  used  to  state  the  relative  dielectric 
constant  of  the  ground  plane  if  the  default  values  mentioned 
above  are  not  utilized. 


/GROUND  (HEIG=10/COND=.002/DIEL=10) 
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11.  INTERVAL  FOR  CALCULATION  This  card  is  used  to  state  the 
number  of  intervals  to  be  used  for  calculating  the  elements 
of  the  impedance  matrix  with  Simpson' s-rule  integration.  A 
large  value  for  the  number  improves  the  accuracy  at  the 
expense  of  greater  execution  time.  For  most  problems  a 
suitable  combination  of  speed  and  accuracy  is  obtained  with 
a  value  of  four,  the  default  value.  If  the  rigiorous 
closed-form  impedance  expressions  in  terms  of  the 
exponential  integrals  is  desired,  set  this  value  to  zero. 

INTERVAL=value 


'lNTE(6) 
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12-  GENERATOR  This  card  is  similiar  to  the  FEED  CARD  in  use, 
except  that  the  segment  numbers  are  stated  instead  of  the 
node  numbers.  This  is  useful  if  three  or  four  segments  meet 
at  a  node.  The  positive  terminal  of  the  generator  is 
connected  to  the  specified  segment  such  that  current  is 
forced  in  the  the  positive  direction.  The  positive  direction 
of  current  flow  is  from  the  first  stated  node  number  of  that 
segment  toward  the  second  stated  as  ordered  on  the 
DESCRIPTION  CARD. 


/GENE  (2,  .5, -9  0/4  ,.5, +9  0) 
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13.  IMPEDANCE  This  card  is  similiar  to  the  LOAD  CARD  in  use, 
except  that  the  segment  numbers  are  stated  instead  of  the 
node  numbers.  As  with  the  GENERATOR  CARD,  this  is  used  if 
three  or  four  segments  are  connected  to  a  node.   The 
impedance  will  be  connected  to  the  positive  terminal  of  the 
specified  segment.  The  format  of  this  card  is  the  same  as 
the  LOAD  CARD. 


'IHPE  (1,1  20, -4 5/3,  120, +4 5) 
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14.  CHANGE  This  card  at  the  end  of  the  data  block  signals 
the  program  that  the  following  data  cards  are  changes  to  the 
previously  read  data,  for  the  next  run.  If  a  "CHANGE  CARD" 
is  used,  the  outputs  must  be  reguested  again  in  the  next 
data  block. 


15.  END  This  card  signals  the  program  that  this  is  the  end 
of  a  data  block  series  and  to  reinitialize  data  for  the  next 
problem.   An  "END  CARD"  cannot  be  used  with  a  "CHANGE  CARD". 


16.  STOP  This  card  signals  the  program  that  all  of  the  data 
cards  have  been  read  and  to  terminate  itself  when  execution 
is  completed.   This  card  must  be  used  as  the  last  card  in 
place  of  the  "END  CARD"  of  the  last  data  block  series.   A 
"STOP  CARD"  cannot  be  used  with  an  "  END  CARD"  in  the  same 
data  block. 
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C      AN  EXAMPLE  PROBLEM 

C 

C      V  ANTENNA 

C 

WIRE(RADIUS  =  1M) 

GEOM (0,-. 18, ♦. 18/0, -.09,*. 09/0, 0,0/0, 0.09,. 09/0,. 18,. 18) 

DESC(l-2/2-3/3-4M-5) 

FEED(3) 

OUTPUT(FARF  =  <t5,50,65,80/STEP  =  5) 

CHANGE 

OUTPUT(BIST=45,45,45,*5/BACK=0,0,10,12) 

OUTPUT(CURRENT) 

CHANGE 

C 

C  CHANGE    STRUCTURE    SHAPE    TO    DIPOLE 

C 

GEOM(0,-. 25,0/0, -.125, 0/0, 0,0/0, .125, 0/0, .25,0) 

PLOT(FARF/PHI=90) 

GROUND{HEIGHT=.2  5/GOOD) 

STOP 


THE  ABOVE  DATA  DECK  WILL  PRODUCE  THE  OUTPUT  ON  THE 
FOLLOWING  PAGES. 
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